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1952 


Uber spezielle Tragfliigelsysteme. 
Von K. Nickel. 


1, Einleitung. Uber die folgenden Untersuchungen wurde in gedrangter Form bereits anlaB- 
lich der wissenschaftlichen Jahrestagung der Gesellschaft fiir angewandte Mathematik und 
Mechanik vom 28. bis 31. Marz 1951 in Freiburg i. Br. vorgetragen (Vortragstitel: Integral- 
gleichungen erster Art in der Strémungslehre und ihre Auflésungen). Sie beschaftigen sich mit 
Tragfliigeln in reibungsfreier, inkompressibler und stationdrer Potentialstrémung. Unter sogleich 
noch naher anzugebenden Voraussetzungen wird darin nach den auf die Tragfliigel ausgeiibten 
Drucken, Kraften und Momenten gefragt. Die Anzahl der Tragfliigel sei n, um den Umfang der 
Betrachtungen nicht zu sehr zu erhdhen, soll der Grenziibergang n > oo (unendlich viele Trag- 
fliigel: Tragfliigelgitter oder Systeme von solchen), ither den in Freiburg ebenfalls berichtet 
wurde, hier nicht vollzogen werden. 

Da bis heute die dreidimensionale Tragfliigeltheorie noch sehr groBe mathematische Schwierig- 
keiten bereitet, soll auch hier eine Beschrankung auf die beiden (im Falle n = 1 vielbearbeiteten) 
Sonderfalle der ,,tragenden Linie“‘ und der ,,Profilumstrémung“ stattfinden, also die Speziali- 
sierung auf verschwindende Fliigeltiefe bzw. unendliche Fliigelspannweite. Durch diese Grenz- 
iibergange kommt man in beiden Fallen auf zweidimensionale Strémungen, deren Behandlung 
durch die Hilfsmittel der Potentialtheorie und konformen Abbildung sehr erleichtert wird. Die 
nachstehend angewandte Methode ist die besonders bequeme Berechnung des Strémungsverlaufs 
aus einer Belegung der Tragfliigelkontur mit Wirbeln, wie sie der Prandtlschen Theorie der 
tragenden Linie und der Birnbaumschen Theorie der diinnen Profile zugrundeliegt. 

Die Beschrankung auf diinne Profile bei der Berechnung der Profilumstrémung wurde aus 
zwei Griinden getroffen. Einmal fiihrt die Birnbaumsche Methode formal auf dieselben Formeln, 
wie die Prandilsche Theorie der tragenden Linie, so daB also bei einer gemeinsamen Behandlung 
die mathematische Arbeit der Lésung dieser Gleichungen nur ein einziges Mal geleistet zu werden 
braucht. Zum anderen ist es zwar, etwa mit konformer Abbildung, jederzeit méglich, mit be- 
liebiger Genauigkeit die Strémung um ein System von Fligelprofilen zu berechnen. Soll aber etwa 
die Wirkung einer kleinen Abanderung dieser Profilkonturen untersucht werden, so ist die ganze 
miihsame Berechnung von neuem auszufithren, ohne dafs Klarheit iiber das zu erwartende Resul- 
tat besteht. Im Gegensatz hierzu sind zwar die Ergebnisse der Theorie der diinnen Profile von 
vornherein mit einem gewissen Fehler behaftet, der von den getroffenen Naherungsannahmen 
herriithrt, jedoch lassen sich allgemeingiiltige Formeln z. B. fiir Auftriebsdichte, Gesamtauftrieb 
und Langsmoment angeben, aus denen sofort der Zusammenhang zwischen Profilgestalt und 
Profileigenschaften abgelesen werden kann. 

Die Anzahl der zulassigen Tragfliigelsysteme muf jedoch noch eingeschrankt werden, da sich 
nicht simtliche méglichen Systeme mit den nachfolgend entwickelten Formeln behandeln lassen. 
Betrachtet man namlich eine Anzahl beliebiger Profilkonturen C mit einer Wirbelbelegung y (z) in 
der komplexen z-Ebene, so wird von diesen Wirbeln an der Stelle € die Stérungsgeschwindigkeit 
u + iv erzeugt, die aus der Formel 


u(t) iv @) = 55 >| oS ae (*) 


berechnet werden kann (es ist iiber alle Konturen C zu summieren). In einer Theorie kleiner 
Stérungsgeschwindigkeiten, wie sie die Prandtlsche und Birnbaumsche Theorie darstellt, inter- 
essiert jedoch hiervon nur die Vertikalgeschwindigkeit v, die durch den Imaginarteil des rechts- 
stehenden Integrals gegeben wird: 


= — Imlay [PEL 
c 
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In den weiteren Rechnungen zu den eingangs gestellten Fragen ergibt sich die Aufgabe, aus 
Gleichung (*) bei gegebener linker Seite die Wirbelverteilung y(z) zu berechnen. Es zeigt sich, 
daB es zwar moglich ist, fiir Gleichung (*) eine allgemeine Auflésungsformel anzugeben, fiir die 
Gleichung (*) ist bis heute jedoch im allgemeinen Fall keine solche Formel bekannt. Nun stim- 
men die beiden Gleichungen iiberein, falls u in der ersten Gleichung verschwindet. Dieser Fall 
tritt fiir beliebige Wirbelbelegungen y(z) jedoch nur ein, wenn alle Konturen C Gerade sind und 
alle auf einer Geraden liegen. Aus diesen Griinden wird daher im folgenden die (sehr starke) 
Einschrankung getroffen, daB® aus der Zahl aller Tragfliigelsysteme nur diejenigen speziellen An- 
ordnungen herausgegriffen werden, bei denen die Tragfliigel (als tragende Linien) nebeneinander 
oder (als Profile) hintereinander liegen. 

Entsprechend den oben angefiihrten zwei Grenzfallen der Tragfligeltheorie teilen sich die 
nachfolgenden Ausfithrungen in zwei Hauptteile: Ziff. 2 befaBt sich mit der Theorie der tragen- 
den Linie, Ziff. 3 mit der Theorie der diinnen Profile je fiir die geschilderten speziellen Tragfliigel- 
systeme, lingere mathematische Betrachtungen sind in den mathematischen Anhang Ziff. 4 ver- 
wiesen, worin auch einige niitzliche Formeln zusammengestellt sind. Zur Einfiihrung in die je- 
weiligen Theorien beginnen Ziff. 2 und 3 je mit dem bereits bekannten Fall n = 1. Nach dem Auf- 
stellen der Grundgleichungen fiir beliebiges n folgt dann in Ziff. 2 der Nachweis, daB es zu beliebig 
vorgegebenen induzierten Abwindgeschwindigkeiten langs der tragenden Linien genau eine zu- 
gehorige Auftriebsverteilung gibt, als Beispiel wird dann der Fall zweier benachbarter, gleich 
eroBer Fliigel mit geringstem induzierten Widerstand durchgerechnet. In Ziff. 3 werden ge- 
schlossene Formeln fiir Auftriebsdichte, Gesamtauftrieb und Langsmoment beliebiger Profil- 
formen angegeben und allgemeine Folgerungen daraus gezogen, die den Zusammenhang zwischen 
Profilform und Profileigenschaft erhellen. Als Beispiel wird noch ein zahlenmafig gegebener 
Doppelfliigel behandelt. 

2. Tragende Linien nebeneinander. a) Kin Tragfliigel erstrecke sich in Spannweitenrichtung 
langs der y-Achse von a bis b (vgl. Abb. 1; die benutzte Bezeichnungsweise folgt dem in der Luft- 
fahrtliteratur iiblichen Koordinatensystem, nach dem die x-Achse in Anstrémrichtung, die 
y-Achse in Tragfliigelrichtung liegt). Die értliche Zirkulation sei [(y), die von dem Tragfliigel her- 
vorgerufene Abwindgeschwindigkeit werde mit w(y) bezeichnet. Unter den Vernachlassigungen — 
der Theorie der tragenden Linie hangen dann nach L. Prandtl! Abwindgeschwindigkeit und Zirku- 
lation nach folgender Gleichung zusammen 

b 


wR ee dI'(n) dy (De 


Gy i 
Dabei ist an der singularen Stelle 7 = y das Integral als Cauchyscher Hauptwert aufzufassen, 
was durch das Symbol ¢ gekennzeichnet werden soll. 
Bei bekannter Abwindgeschwindigkeit w(y) stellt (1) eine lineare Fredholmsche Integral- 
gleichung erster Art fiir dI’/dy dar. Ihre Lésung wurde zuerst von A. Betz? angegeben. Sie 
lautet 


b 
dI'(n) 4 1 w(y) 
a r ——— 7 (b — —a) dy}. 
Gone Goes me eee (2) 
Dabei ist die Gesamtzirkulation I eine beliebige Konstante mit der Bedeutung 
b 


nS (eae geal iae 


Da die Zirkulation [’(y) an den Tragfliigelenden verschwinden soll, da also ["(a) = Eb) 0 
sein soll, findet man [’ = 0. Damit ergibt sich durch Integration von (2) 


n A ; 
= EO) 4 1 ie 
a | di See | TCEEn an (¢ 7 — V(b — y)(y—a) dy) dt. (3) 


a a 


* L. Prandtl u. A. Betz, Vier Abhandlungen zur Hydrodynamik und Aerodynamik. Neudruck aus den 
Verhandlungen des III. Internationalen Mathematiker-Kongresses zu Heidelberg und aus den Nachrichten 
der ig gh der Wissenschaften zu Géttingen. Géttingen 1927. 

. Betz, Beitrage zur Tragfliigeltheorie, mit besonderer Beriicksichti i i 
Teves Die ae ee efliig 5 ricksichtigung des einfachen rechteckigen 
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Hierin 1aBt sich die Integrationsreihenfolge vertauschen und das innere Integral geschlossen be- 
rechnen. Man findet so die bekannte Betzsche Umkehrformel 


b 
Pa) == | wo)los5 


a 


b 
abo yn) tym Veo—y) (y a) (@—n)(q—a) | 
— dy. A) 


b) Betrachtet man n Tragfliigel nebeneinander lings der y-Achse statt des einen der Abb. 1, 
von denen sich der y-te von a, bis b, erstrecken mége (vgl. Abb. 2), so ist die Gleichung (1) zu er- 
setzen durch das entsprechende Gleichungssystem 


by 
ess dIy(qv) dy» # 
WwW, (¥,) s— ease dn» Yu—Nr (u Ik Dai ANDeO n) . (5) 


Dabei und im folgenden bezieht sich der Index immer auf den p-ten Fliigel, so daB z. B. w,(y,,) 
die Abwindgeschwindigkeit an der Stelle y, ist, die demnach zwischen a, und b, liegt. 


Abb.2. 


Entsprechend zu oben kann man auch hierin die Abwindgeschwindigkeiten w,(y,,) als be- 
kannt voraussetzen und nach den dadurch bestimmten Zirkulationen I’, (7,) fragen. Bei dieser 
Auffassung ist (5) ein Integralgleichungssystem. Seine Lésung ist unter sehr allgemeinen Voraus- 
setzungen + 


rr) 2 (lye. 


= x 
dn, 4 n 
| |/— i @.—my(ox— m0) 
“= 6 
n by ( ) / = ( ) 
| Prat + SIG FOV — Hb, —yN(aa— I) yu] OLB) 
a u—te Vo gaa 
ft: 
Dabei ist 


nite 
—1 Ze n— } dr, (ne) z 
P,,1 (1) =! A "ov Z | ae oa 
aa : 


ein Polynom vom Grade n—1. Infolge der Bedingung ; 
ese (Oe or, (ar) = La (ba) .0 (7) 


verschwindet analog zu oben in P,,_; das héchste Glied, die Summe der Gesamtzirkulationen, 
so daB dieses Polynom also tatsadchlich héchstens vom n—2-ten Grade ist. Durch Integration 
itber (6) erhalt man wie oben 


i (1) 
ay \/— ire. — ty) (ax —t) (8) 


: f auld fp 
ue Wul¥u ee 
x Pao (t,) ne, y) ¢ RL, Sj LTO JaNe Yu) dy,, | dt,. 
ou 

1 K. Nickel, Math. Z. 54 (1951), S. 81. Es soll an dieser Stelle ausdriicklich darauf hingewiesen werden, 
daB (6) das System (5) auch bei beliebigen komplexen Integrationsgrenzen a,, b, und Funktionen Wu 
auflést, ferner auch fiir beliebige Integrationswege zwischen a, und by. An der angegebenen Stelle ist dies 
noch nicht bemerkt worden. 

26% 
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Im Gegensatz zu dem obigen Falle n = 1 ist jetzt eine weitere Vereinfachung nicht mehr még- | 
lich, da die Vertauschung der Integrationsreihenfolge auf hyperelliptische Integrale fithrt, die 


sich im allgemeinen nicht geschlossen auswerten lassen. 


c) In den Koeffizienten des Polynoms P,,; stehen n willkiirliche Konstante zur Verfiigung, — 


um zu erreichen, da® die Bedingung (7) erfiillt wird, da also simtliche Zirkulationen I, an den 
Fliigelenden verschwinden. Es sol] nun zunachst gezeigt werden, daB diese Konstanten tatsachlich 
in jedem Falle ausreichen, aber auch alle gebraucht werden, daB also durch die Bedingung (7) 
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die Einzelzirkulationen I,(7,) aus Gleichung (8) eindeutig bestimmt sind. Dabei braucht man © 
sich nicht darum zu kiimmern, daB es ja schon bekannt ist, da der hiéchste Koeffizient von P,_; _ 


durch seine spezielle Form von selbst verschwindet. Fiir den nachfolgenden Beweis hat man also 
iiber n willkiirliche Koeffizienten von P,_ zu verfiigen. 

Nach (8) ist stets von selbst I',(a,) = 0, d.h. die Zirkulationen I, verschwinden an den 
linken Fliigelenden. Also entsprechen den 2n Bedingungsgleichungen (7) nur noch die n-Glei- 


chungen 
b, b, ; n 
fa, f tet) ay, A Say 
ean ae a 


b 


2 \—1e—m@—m \2 
4 x=1 “u 


et) Ie Ro if 
Dies ist ein lineares, unhomogenes Gleichungssystem fiir die n unbekannten Koeffizienten in Pan 


Ein solches hat bekanntlich genau dann eine Lésung, wenn die zugehérige Determinante nicht 
verschwindet. Diese lautet bis auf eine Konstante 


l “ Cael 

| Yo 

| : = any (1; OSL 2st) 
pea nares Necrerd | 
9 I 

und ist, wie in Ziff. 4a, gezeigt wird, immer > 0. Also gilt das Ergebnis: Schreibt man vor, daB 


D 


die Zirkulationen J’, an den Fliigelenden verschwinden sollen, so gibt es zu beliebigen Abwind-_ 
geschwindigkeiten w,(y,) immer genau einen zugehérigen ,,Satz‘° Zirkulationsverteilungen — 


E 
: : . wu(¥u) ; n | 
: | i¢ Yu" |-# Dj) a) (Aree) dy, dy, 


I’, (7,), der das Gleichungssystem (5) befriedigt und dieser wird durch die Formel (8) geliefert. 


Die unbekannten Koeffizienten des Polynoms P,; sind dabei durch die Bedingung (7) ein- 
deutig bestimmt. 

d) Berechnet man den induzierten Widerstand dieser n Tragfliigel, so kommt man nach 
Prandtl} bis auf eine unwesentliche Konstante auf den folgenden Ausdruck 


W= > [ Pio) eH) ay (9) 


y=] 


Die dritte Grundaufgabe der Tragfliigeltheorie stellt nun die Aufgabe, die J’, so zu bestimmen, 
da8 der Widerstand W zu einem Minimum wird. Als Nebenbedingung wird darin die Forderung 
gestellt, daB der Gesamtauftrieb 

n by 

A= >’ | r,(y,) dy, (10) 
yv=1. 


a 


dabei konstant gehalten werden soll. 
Dieses Problem wurde von M. Munk? gelist. Die giinstigste Zirkulationsverteilung ergibt sich 
nach dieser Arbeit aus der Forderung, daB die Abwindgeschwindigkeiten w,(¥,) konstant und 


LAS a OBS O24: 
* M. Munk, Isoperimetrische Aufgaben aus der Theorie des Fluges, Diss. Géttingen 1919.1 
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gleich groB sein sollen. Setzt man danach w, = w, =--+ = w, = konst. in (8) ein, so ergibt sich 
mit der Formel ! 


n 


by 2 Cae 
wu ¥u) " n 
Dag ee V — ie, —y,) (a, —y,) dy, = (— "Qa s (11) 


= 


wobei Q, ein normiertes Polynom vom Grade n ist, 


it, (12) 


I’, (n,) = € | yy Es) 
~ ee IL (bx—t)(ax—b) 
wo R,, ein normiertes Polynom vom Grade n bedeutet. Die Konstante C ergibt sich dabei aus der 
Bedingung (10), die n Koeffizienten von R, sind wieder wie oben so zu bestimmen, daB die Zirku- 
Jationen an den Fliigelenden verschwinden. Fiir n = 1 at sich die Integration in (12) ausfiihren 
und man kommt auf die bekannte ? elliptische Auftriebsverteilung. Fiir n > 1 handelt es sich 
in (12) um hyperelliptische Integrale und eine Zuriickfithrung auf elementare Funktionen ist 
nicht mehr méglich, jedoch wurde naherungsweise von H. Schlichting ® fiir verschiedene Trag- 
fliigelsysteme die Abminderung des induzierten Widerstands in Abhangigkeit vom gegenseitigen 
Abstand der Tragfliigel untersucht. Fiir n = ov, alle Tragfliigel gleich breit und mit gleichen 
Seitenabstanden, hat A. Betz‘ die dritte Grundaufgabe gelést. In diesem Extremfall lassen sich 
Zirkulationsverteilung und Widerstand wieder durch elementare Funktionen darstellen. 
In dem Ausdruck (9) fiir den induzierten Widerstand kann man wegen w,,(y,,) = konst. noch 
eine Integration sparen. Integriert man ndmlich (9) partiell und beachtet dabei (7), so findet 
man mit Konstanten C, und C, 


b, b, 
> dl, > y Rn(yv 
WV = C, [> dy, dy, = C, | i - y (y») dy, . (13) 
ee ha, |/— Hf bx—y0) (ax 9») 


e) Im Falle n = 2 ergeben sich fiir Zirkulation (12) und induzierten Widerstand (13) elliptische 
Integrale. Diese sind ausfiihrlich tabuliert *), und erleichtern damit die Berechnung von /’, und 
W. Da der symmetrische Fall (Tragfliigel mit Spalt) ein gewisses praktisches Interesse hat, soll 
er hier als Beispiel gebracht werden. Es bereitet keine Schwierigkeiten, die noch offenen Kon- 
stanten in (12) und (13) zu bestimmen und die auftretenden Integrale anschlieBend auf die 
Legendresche Normalfoim zu bringen. Da die Rechnungen jedoch etwas umfangreich sind, soll 
hier nur das Ergebnis angegeben werden: 

Gibt man den beiden Tragfliigeln die Breite 1 und bezeichnet die Spaltbreite zwischen ihnen 
mit s (vgl. Abb. 3), so ergeben sich die beiden Auftriebsverteilungen 


ji+s \1-Fs ( ey ( Vile ) 
al ra)Bt So et eam RE Sd p 
(2428-499) K(2 UES) a -harel2 ae 


P(y)=— A +3) 


t 
wn 
Se 


2-+s 
mit 


1 ace 3 
y=tt@ts 14g $s ante. 0Sp Sale. 


Dabei ist A der konstantgehaltene Gesamtauftrieb beider Fliigel, E(k, p) und F(k, p) sind die 
elliptischen Legendreschen Normalintegrale, E(k) und K(k) die vollstindigen elliptischen Inte- 
grale. In den beiden Grenzfallen s = 0 und s = o ergeben sich, wie zu erwarten, elliptische Auf- 
triebsverteilungen, im ersten Fall iiber den ganzen Tragfliigel hinweg, im zweiten Fall itber jeden 


1 K. Nickel, a. a. O. S.90, Formel (17). 

2 L. Prandtl, a.a.O. S. 32. . 

3 H. Schlichting, Mitteilungen der Deutschen Akademie der Luftfahrtforschung, Heft 2. 

4 4. Betz, Der Energieverlust an spaltférmigen Unterbrechungen von Tragfliigeln (maschinenschrift- 
lich). 
; A. M. Legendres Tafeln der elliptischen Normalintegrale erster und zweiter Gattung. Herausgegeben 
von F. Emde, Stuttgart 1931, 
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er beiden Einzelfliigel. Der Zusammenhang zwischen induziertem Widerstand W und s ist 


W= W, rileees ==) on 
1l+s 
2x (VEE) 
Bes) : ( per 
+2649) K(2 sag | = Co RE eos 
Als Grenzwerte ergeben sich W/W, = 1 fiir s = 0 und W/W, = 2 fiir s = oo. Dazwischen steigt 
die Kurve monoton mit monotoner Ableitung. In Abb. 3 sind einige Auftriebsverteilungen (ge- 


strichelt die elliptischen zum Vergleich) und der Verlauf von W/W, gezeichnet. Man sieht, daf 
die Auftriebsverteilungen schon fiir sehr kleine 


ae 


Wo 


a ed Spaltbreiten nur noch wenig von den elliptischen 
oe re Verteilungen abweichen. Damit erklart sich das 
$=0 sehr starke Ansteigen des in- 
duzierten Widerstands bei klei- ~*~ é ; 
nen Spaltbreiten. Z. B. ist bei “3 
7 einer Spaltbreite s = 1/100 , 


(5cm bei 10m Spannweite) der Dbbas: 


induzierte Widerstand schon 
um 41% héher als ohne Spalt. Dazu muB allerdings 
bemerkt werden, dain diesem Fall bei den itblichen 


Fliigeltiefe: ein Meter = 20fache Spaltbreite) die 
Voraussetzungen der tragenden Linie nicht mehr 


! | 
0 OS 40 15 6,0 Go 
Abb. 3. 


erfiillt sind, so daB eigentlich dreidimensional gerechnet werden miiSte. Immerhin gibt aber die 
Theorie der tragenden Linie doch schon einen gewissen Einblick in die zu erwartenden Gesetz- 
maBigkeiten. Weiter ist noch darauf zu achten, daf es sich bei den berechneten Werten von W/W, 
um Geringstwerte handelt, die nur fiir die dargestellten Auftriebsverteilungen giiltig sind, daB 
also alle anderen Auftriebsverteilungen mit Spalt noch héhere Widerstandswerte besitzen. Immer- 
hin wird man jedoch nach den experimentellen Ergebnissen im Falle eines Fliigels vermuten, 
da sich ebenso wie dort die tatsachlichen Widerstandswerte von den Geringstwerten nicht stark 
unterscheiden werden. Nachfolgende Tabelle gibt einige Werte von W/W,. Eine solche Tabelle 
wurde zuerst von K. Pohlhausen berechnet und von L. Prandtl mitgeteilt. 


Ss 0,01 | 0,02 | 0,03 | 0,04 | 0,05 | 0,06 0,07 | 0,08 | 0,09 | 0,10 | 0,15 | 0,20 | 0,25 


W/W, | 1,412 | 1,470 | 1,510 | 1,540 | 1,566 | 1,587 1,606 | 1,623 | 1,638 | 1,652 1,706 | 1,745 | 1,776 


; 0,30 | 0,35 | 0,40] 0,45 | 0,50 | 0,60 | 0,70 | 0,80 | 0,90] 1,00] 1,50] 2,00] 2,50 


W/W, | 1,800 | 1,820 | 1,837 | 1,851 | 1,863 | 1,884 1,899 | 1,912 1,922 | 1,931 1,957 | 1,971 | 1,979 


: 3,00 | 4,00 | 5,00 | 10,00 | co 
W/W, | 1,984 | 1,990 | 1,993 | 1,998 | 2,000 


1 Siehe S, 52 der in FuBnote 1 von S. 364 zitierten ,,Vier Abhandlungen“. 
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f) Der in e) behandelte Tragfliigel mit Spalt ist aerodynamisch dasselbe Gebilde, wie es zwei 
nebeneinanderfliegende Tragfliigel sind. Nur stellt man bei einem Blick auf die Abb. 3 fest, daB 
die in e) gefundenen Zirkulationsverteilungen fiir die benachbarten Tragfliigel nicht zulassig sein kén- 
nen, da (bis auf den Fall s = oc) das Rollmoment um die Mitte jedes der beiden Einzelfliigel nicht 
verschwindet, so daB die Tragfliigel sich in Drehung versetzen wiirden, wie das in Abb. 4 skizziert 
ist. Will man also fiir zwei freie, symmetrische, im Abstand s befindliche Tragfliigel diejenigen 
Zirkulationsverteilungen bestimmen, die ein Minimum des induzierten Widerstandes ergeben, 
so mufs man aufer konstantem Gesamtauftrieb noch konstantes (und zwar verschwindendes) 
Rollmoment vorschreiben. Diese Aufgabe wird zwar bei M. Munk nicht behandelt, sie 1aBt sich 
aber in Erweiterung eines fiir den Einzelfliigel ausgesprochenen Ergebnisses des Verfassers 1 
leicht lésen. Man erhalt als Bedingung fiir die gesuchten giinstigsten Auftriebsverteilungen, daB 
die zugehérigen Abwindgeschwindigkeiten w,(y,,) linear in y,, sein miissen, und zwar ergibt sich 


w;(y¥;) =4—By,, wWe(yo) =O + B ¥2 


(fiir mehr als zwei Tragfliigel entsprechend). Dabei sind die Konstanten « und f aus den beiden 
Forderungen vorgegebenen Gesamtauftriebs und verschwindenden Rollmoments zu bestimmen. 

Mit Hilfe der Gleichung (8) kann man auch hier sofort bis auf die noch offenen Konstanten 
die Lésung dieser Aufgabe hinschreiben. Nur bereitet jetzt die zahlenmaige Auswertung der 
Integrale erheblich mehr Schwierigkeiten als fiir das Beispiel e). Deshalb soll hier von einer Durch- 
rechnung abgesehen werden. Man kann sich jedoch iiber die Verhaltnisse wenigstens qualitativ 
Klarheit verschaffen: Durch die zusitzliche Bedingung, da8B die Rollmomente fiir die beiden 
Einzelfliigel verschwinden miissen, werden sich fiir s + 0 die Zirkulationsverteilungen noch mehr 
den Ellipsen annahern. Weiter wird die neue W/W,-Kurve iiber der bisherigen liegen, da der 
induzierte Widerstand sich nur vergréfern kann, wenn, wie hier, zusdtzlich noch weitere Neben- 
bedingungen einzuhalten sind. Da fiir s = oo die elliptischen Teilauftriebsverteilungen auch die 
neue Forderung von selbst erfiillen, findet man hier dieselbe Auftriebsverteilung und also auch 
denselben Wert W/W, = 2. Der mutmaBliche Verlauf der Kurve W/W, ist in Abb.5 skizziert, 
jedoch ist dabei der Punkt s = 0 exakt bestimmt. Die Forderung verschwindenden Rollmoments 
fiir jeden Einzelfliigel lautet hier namlich 


1 
1 
[ Poi(\y! =) dy =0 
—1 
und fiihrt auf die Abwindbedingung u(y) = «4, +f, | y | mit zwei noch offenen Konstanten a, 


und f,. Die hierzu gehorige Auftriebsverteilung 1aft sich nun geschlossen angeben und ist schon 
an anderer Stelle? berechnet worden. Die Konstantenbestimmung liefert 


wi=z4ly 4+ (32 ®)|y 


4 


und 
i . ENG peses 
ity) = atk (4) yl—y +(37—8)y*log A =F) 
vgl. Abb.5. Dabei ist A der vorgeschriebene Gesamtauftrieb beider Fliigel. Mit diesen Funktionen 
findet man schlieBlich 


nn 4 (24 —16m + 32%) =1,2537..., 
wobei hier W, der Wert der elliptischen Auftriebsverteilung gemaf e) ist. 

g) Die vorstehend abgeleiteten Formeln fiir /’, und W lassen sich natirlich auch auf andere 
Probleme anwenden, so z. B. auf die zweite Grundaufgabe, die auf die Prandtlsche Integro- 
differentialgleichung fithrt. Bei solchen Anwendungen wird es jedoch meist ecinfacher sein, bei 
gegebenen Abwindgeschwindigkeiten w, (y,,) die Zirkulationen J”,(y,) nicht exakt nach der For- 
mel] (8) zu berechnen, sondern J’, in eine Reihe zu entwickeln, von der man dann, je nach der 
gewiinschten Genauigkeit, mehr oder weniger Glieder beriicksichtigt. Am zweckmaBigsten diirfte 


1 K. Nickel, Z. angew. Math. Mech. 31 (1951), Sri 
2 K. Nickel, a. a. O. (vorhergehende Fufnote) 5. 77, 
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wohl fiir die meisten Falle der folgende Ansatz sein: 


rn) = (VV — BG) (a — ad Ootamtrath +o tnt) @=12 ms (4) 


wobei also die Beiwerte y,, der Reihe von der Tragfliigelnummer unabhangig sind. Der Ansatz (14) | 


entspricht vollkommen der im Falle n = 1 iiblichen Setzung. ‘Man kommt so auf sehr einfache 
Integrationsformeln, bei denen nur begreiflicherweise die Integrationskonstanten, die von den 
Tragfligelgrenzen a, und b, abhangen, fiir gréBere n recht komplizierte Ausdriicke werden (man 


kommt auch hierbei auf hyperelliptische Integrale). Sieht man jedoch von ihrer genauen Be-. | 


stimmung ab, so findet man aus den vom Verfasser } angegebenen Formeln die sehr einfachen 
Ausdriicke 


n 


| — IT (by —)(a—n,) dn, = (—1)" Pose (y,) (b=0,1,..-) (15) 


Jk 
ny = 


b, 
dM eee 
0 af eile lv A=1 

a, 


(P,,., ein normiertes Polynom vom Grade n + k) und 


1 n by yntk d (+) (k 7 0, Me ONG 5) 
poiG : a =(-1"" (16) 
7 Yu—"s n 
g 2, 16: =») (a2—ty) 0 (k= —n,. , —}) 
(Q; ein normiertes Polynom vom Grade k). Setzt man (14) in (5) ein, so findet man mit (15) und (16) 
by 
ula dln) any a wae 
WwW, (¥,) ee an a, g dn» Yu—N AS Riek Cy) (17) 


wo wieder R,_11, ein normiertes Polynom vom Grade n —1 -+ k ist. 

Im jeweils vorliegenden Einzelfall lassen sich die Koeffizienten von R,_ 14 ;, oft aus bekannten 
Nebenbedingungen, wie z. B. Symmetrieeigenschaften usw., ohne groBe Rechnung bestimmen. 
Jedoch mu man sich immer auf erheblich gréBere Rechenarbeit gefaBt machen, als im Falle 
n— 1. Denn fur tk — 0, also 


n, 


Enea eaiy | — I,m) (sm), 


erhalt man aus (17) fiir die Abwindgeschwindigkeiten w,(y,) schon ein ,,Nullpolynom“ (mn —1)- 
ten Grades, auf dessen Koeffizienten man keinen Einflu8 besitzt. Man muB also im Ansatz (14) 
erst einmal eine Anzah] Glieder mitnehmen, nur um dieses Nullpolynom wieder zu kompensieren 
und benétigt dariiber hinaus noch weitere Glieder, um die vorgegebenen Abwindgeschwindig- 
keiten w,(y_) geniigend genau annahern zu kénnen. 

3. Diinne Profile hintereinander. a) Ein Tragfliigel erstrecke sich in Tiefenrichtung in Rich- 
tung der x-Achse von a bis b und werde von links mit der Geschwindigkeit u,, angestrémt (vgl. 
Abb.6). Unter den Vernachlassigungen der Theorie der diinnen Profile 2? bestimmt sich dann die 
Starke der Wirbelbelegung y der x-Achsé aus der Gleichung 


b 
df(x) ss y(é) 
ig Bae. CF earie (18) 


Dabei ist f(x) die bekannte Ordinate der ,,Profilskelettlinie“ (vgl. Abb. 6). Die Integralgleichung 
(18) geht in die Gleichung (1) itber, wenn man df(x)/dx durch w(y) und y (£) durch — = u,, dI"(n)/dy 


ersetzt. Also kann man auch hier die Auflésungsformel (2) gebrauchen. Nur ist jetzt darauf zu 
achten, da y(&) wegen der Joukowskischen AbfluBbedingung so zu bestimmen ist, da keine 
Umstrémung der Hinterkante eintritt, daB also (b) = 0 ist. Bestimmt man J’ in (2) aus dieser 
Forderung, so findet man folgende Auflésungsfoimel der Gleichung (18): . 


b 
Sele bee of) a ne ‘ 
7 (5) te 1/3 ¢ dx E—x =e (19) 


1Siehe FuBnote 1 von S. 365, > W. Birnbaum, Z. angew. Math. Mech. 3 (1923), S. 290. 
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Wegen der Kutta-Joukowskischen Formel: Auftrieb = 9-u,-y (0 = Luftdichte) ist mit (19) bis 
auf einen konstanten Faktor auch schon die Auftriebsdichte bestimmt. Gesamtauftrich und 
Langsmoment, etwa um den Punkt x = 0, ergeben sich aus (19) sofort durch Integration. Ver- 
tauscht man die Integrationsreihenfolge und fiihrt die innere Integration aus, so findet man den 


Auftrieb 


b b 
: ay 
A=ou,| (d= —29ut, | LO) /2—* as (20) 
und das Langsmoment (um den Nullpunkt) 
b b 
Meu, | y (E)EdE =2out, fae es *| =F. (21) 


Aus den Formeln (19) bis (21) lassen sich sofort fiir belie bige Profilmittellinien f(x) die wichtig- 
sten Eigenschaften ablesen. 


Sta) be 
Ee eee an on oe 
= oh es ees 
= Tt 


Protilkontur. 
leo 4 = ay d v. 
v 


We 
Abb. 6. Abb. 7. 


b) Entsprechend dem Schritt in Ziff.2b, von einem zu n Tragfliigeln kann man auch hier zu 
n hintereinanderliegenden Profilen iibergehen. Fiihrt man die entsprechende Bezeichnungsweise 
wie dort ein (vgl. Abb.7), so ersetzen sich die Gleichungen (18) bis (21) durch 


by 
df (%) 1 - (eS 
Xu ane DX BAU (Gl, Ze) (22) 


dxy gy — Xu 


mit det Auflésung 


n n by, n 
Loeusil, Polen df (Xp) 1 \ a, — Xp 3 
pe) | Se > 6 Hee | KES One 
oy : 


und dem Gesamtauftrieb 


n n by / n 
a: ay - 2 df. (x)u |/ a, — Xu 
A — OS [> () dé, ea 20 es zy | pp | To ines, dx, (24) 
v= = = 

ay 


ay 


sowie dem Moment um den Nullpunkt 
1 > i step Xp 
7 (b, ——— ay) <= Xy (ea 18 byes te dx, . 
g=1 eles 


* by n Pu 
df u (Xp 
My= eu. >) | 70(6)E d= Zou, Df SR 
p=l a, fed hy : 5 (25) 


c) Aus diesen Formeln (23) bis (25) sollen nun einige allgemeingiiltige Aussagen gefolgert 
werden (giiltig im Rahmen der Theorie der diinnen Profile), die z.T. in Sonderfallen schon be- 
kannt sind 1, Ihr Beweis folgt in d). 

Dabei ist unter dem Anstellwinkel « immer der Anstellwinkel aller Profile des Systems ver- 
standen (vgl. Abb.7); es soll also nicht etwa ein Einzelprofil aus dem Verband gelést und die 
Wirkung einer Drehung dieses speziellen Profils untersucht werden. Zu Vereinfachung fiir einige 
der folgenden Aussagen soll ferner noch die Vereinbarung getroffen werden, dafi die Hinterenden 
der einzelnen Profilskelettlinien alle auf einer Geraden liegen sollen. Diese Voraussetzung ist 
durchaus zulassig, denn in der Theorie der diinnen Profile kommt es, wie ein Blick auf die 


1 Man vergleiche fiir n = 2 etwa R. Grammel, Die hydrodynamischen Grundlagen des Fluges. S. 76 ff. 
Braunschweig 1917, ferner H. Séhngen, Luftfabrt-Forsch. 17 (1940), 5. 17. Siehe auch K. Nickel, Z. angew. 
Math. Mechan. 31 (1951), 5. 297. 
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Gleichungen (22) und (23) lehrt, nur auf die Neigung f, der Skelettlinie an, nicht dagegen auf 
ihre absolute Hohe f,. Man kann also innerhalb der getroffenen Voraussetzungen (geringe Héhen- 
differenzen der Skelettlinie) die einzelnen Fliigel beliecbig nach oben oder unten verschieben, kann 
also die obige Forderung f,(b,) = 0 immer erfiillen. 

1. Auftriebsverteilung ! c,(x) = 2y(x)/u., Gesamtauftrichb A und Moment M, sind linear 
im Anstellwinkel «, es gilt also 


dca( xv 
Ca (2) = Ca (%,) + A , 


3. Es gilt allgemein 


4. Es gibt genau einen Profil-Neutralpunkt N (aerodynamic center), um den das Moment M 
anstellwinkelunabhangig wird, d.h.: bezogen auf N ist dc,,/da = 0. Man erhalt diesen Neutral- 
punkt aus dem Schwerpunkt des homogen mit Masse belegten Fliigelsystems, verschoben um 
1/, der Summe der Fliigeltiefen in Anstrémrichtung. 

5. Fir f,(x,) = 0 ist ebenfalls c,, => 0. Dabei tritt das Gleichheitszeichen bei c,, nur fiir 
f.(%,) = 0 ein. Oder: Legt man in Anblasrichtung durch den hintersten Punkt des Fligelsystems 
eine Gerade, so kann diese bei Gesamtauftrieb Null nie ausschlieBlich auf einer Seite des Systems 
liegen. 

6. Wenn sich die Tangente der Fliigelkonturen vom vordersten bis zum hintersten Punkt 
immer gleichsinnig weiterdreht, so ist c,,, + 0 (c,, bezogen auf den Neutralpunkt). Oder: Druck- 
punktfestigkeit ohne Wendepunkt im Profil ist nicht méglich. 

7. Sonderfall der Streckenprofile mit demselben Anstellwinkel: Fiir diese ist c,(x,) = ¢., = 
Cm, = 0. Der Gesamtauftrieb ist bei allen Anstellwinkeln derselbe, als wenn die einzelnen Fliigel 
zu einem einzigen Fliigel zusammengeriickt waren; der Druckpunkt liegt im Neutralpunkt. 

d) Der Beweis der ausgesprochenen Behauptungen 1 bis 7 kann wie folgt gefiihrt werden: 

1, Die Neigung der Skelettlinie f(x) mit dem Anstellwinkel a wird durch 


f(x) —o (26) 
angenahert. (Drehung des Koordinatensystems um den Winkel a, fiir die betrachteten kleinen 
Winkel « kann sin a durch « und cos « durch 1 ersetzt werden.) Aus der Linearitat der Ausdriicke 
(23), (24) und (25) folgt damit sofort die Behauptung. 

2. Setzt man (26) in (23) ein, so findet man 


b 
i D 
dea(ér) 2 dy(&,) 4 BY ik SS I 1 ray, — Xp 
do We) Tdtcow We ie 2 E, — xy ag! oem dx, 
= a 
lb 


2 ey eee 
Ea ies 


nach Formel (32) [vgl. Ziff. 4d)]. 


' Die Verwendung von ca(x) zur Kennzeichnung der ortlichen Auftriebsdichte und von ca als Beiwert 
des Gesamtauftriebs wird kaum zu MiBverstandnissen AnlaB geben. In der Theorie der tragenden Linie 


wird ebenfalls seit langem schon ca als Beiwert des Gesamtauftriebs neben ca(y) als drtlichem Auftriebs- 
beiwert benutzt, 
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3. Setzt man (26) in (24) ein, so kommt mit (29) aus Ziff. 4d) 


dca = 
i See Foe dx, = 27. 
y (by ==) A=1 bs Sade 


v=1 


4. Die Abszisse des Neutralpunktes sei «y. Fiir das Moment um diesen Punkt findet man 
dann aus (25) und (26) 


n 


dMy_ d = 
Gt ran ote | Eile rae de 


v 


zie fs +90 = 
=m0u2 en (ene ApS 6-0 32> meee 


Soll nun dMy/d« = 0 sein, so findet man hieraus, da® das nur fiir 


> (63 —a? : 
Eoi—a) 


Peas 
ya has — 7 D>, —4,) 


> (by — ay) Via 
méglich ist. Hierin stellt das erste Glied den Schwerpunkt des Tragfliigelsystems dar, wenn man 


sich die Strecken (a,, b,) gleichmafig mit Masse belegt denkt, der zweite Summand ist ein Viertel 


der Summe der Fliigeltiefen. 
5. Durch ae Integration aus (24) erhalt man 


ys df, (%y) nee iy Ree _ fe a yy ee a,a— % g 
A——2@u, oie xy | 1 0 A exe es = Si (x Vea T= . 


falls f,(b,) = 0 ist, me eingangs verabredet wurde sh Der eee 


cay) Tie 
dx, a2 veer. 


im Integranden des letzten Integrals ist langs der Fligel itberall > 0, wie in Ziff.4b) bewiesen 
werden wird. Also ist das Vorzeichen von f,(x,) fiir das Vorzeichen von A und damit von c,, 
verantwortlich. Insbesondere kann offensichtlich c,, = 0 unter der Bedingung f, > 0 nur im 


Falle f, = 0 ane 
6. Nach (24), (25) ee 4. ist das Drehmoment um den Neutralpunkt 


mm en. [vt ) (&, — xy) dé, = 
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Dieser Ausdruck soll durch partielle Integration umgeformt werden. Dazu setzt man zur Ab- 


kiirzung 


b, x 
ia : x (62 —a2) A / 
imeine. = 1 a,—é 
U(x )= >) | ot | 14-2 —— + a) ee 
res A 
ia 2d (bo — a) eat 
e o=l1 
dy ay, 


Weiter soll im folgenden zugelassen werden, daB die Profilskelettlinien in endlich vielen Punkten 
geknickt sind. Seien dies die Stellen ce (v= 12; 5. ms ape 1, 2, eh) eo daB also 
fit (Cu,) +Fp— (C,,) ist (unter f{. und f” sei die obere bzw. untere Ableitung an der betrachteten 


Stelle verstanden). Partielle Integration von My liefert dann 


n b, 2 Le 
My = 2gu2,) >! { 7 U(x) dx, + >) U (a) fs (a) —fe (b—a)] + | 
y=1 a, y=2 


id 


ee ucala oo (1 


v= vr=1 


wenn man beachtet, daB U(a,) = U(b,_,) (v = 2, 3,...m) ist und daB U(a,) = U(b,) = 0 gilt. 
Dreht sich nun die Tangente der Skelettlinie gleichsinnig beim Fortschreiten von a, nach b,, 

etwa im mathematisch positiven Sinne (Profile 
| nach oben gewélbt), so ist iiberall 


| d? v 7 ‘ 
| dee £0, file) <fiO-a) 


2 
dx? 


und Se (C,,) =< Su (C,,) j 


Wie in Ziff.4c) gezeigt wird, ist nun immer 
U(y,) <9, wobei das Gleichheitszeichen nur 
\ fiir y, =a, und y, =65, eintritt. Nach Glei- 
chung (28) wird dann auch My => 0 und damit 
auch c¢,,, > 0. Dabei kann das Gleichheitszeichen 
offenbar nur fiir f = 0 eintreten. 

7. Fir die Streckenprofile gilt f,(x,) = 0, also 
wird (26) zu —a. Setzt man (26) in (23), (24) 
und (25) ein, so kommen mit (29) aus Ziff. 4d) 
sofort die Behauptungen. 

e) Als Beispiel soll ein Doppelfliigel betrachtet 
werden, der sich von x = 0 bis x = 3 und von 
S (a) =x (3x07-2o*-I2 +678) < x—=4 bis x = 5 erstreckt (vgl. Abb.8). Die Profil- 
Cg 48 V2(3-2)(4- 26-2) (z+ 7) Ss mittellinie sei innerhalb dieser beiden Strecken 
LI oh Ieee gegeben durch 
nao ee, f(x) = (3 x3 — 20 x2 — 147 x + 678) . 

Mit den Formeln (32) bis (35) des mathematischen 
Anhangs ergeben sich dann die in Abb.8 dar- 


gestellten Verteilungen fiir ¢,,(x) und de,(x)/dax. Die Profilmittellinie ist so gewahlt, daB c,,(x) 
gleichzeitig an den beiden Eintrittskanten x = 0 und x — 4 verschwindet (stoBfreier Eintritt) 


! 
| 
| 
| 
| 
| 
{ 


Abb. 8. 


4, Mathematischer Anhang. a) Behauptung: Fiir folgende Determinante gilt 


ne | b,, ng} dn, | 


| | |/— 8 6:—m) @—m | 


la 
| 


=!) (50) ea Ld eee 


v 
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Beweis: Bekanntlich gilt fiir die Vandermondesche Determinante 1 


ea lage IT me 1 Oe. 
Ee 
Also ist auch 


IT (Hv — Hu) 


o—1 T uy v= 


vie | ree 
I hi be ea ey girs >0, 
i i ees aaa, (a, — np) fe |/—He—m (a, —7) 

= — =a 


und damit wird 


by | e bn 0 
ae | ; ; “A HT (in — nu) 
iD) 2 — dy, he ee : j = dy,-+-dn,>0. 
|/ = FE bimm) (am) | . y Ht | — 1 (bs) (a2) 


b) Behauptung: 


ae eae Se 
d QA — Xp 
/ dsp ae 


Beweis: Offenbar geniigt es, nachzuweisen, daB gilt 


d 


nh 
dx, IT 


Wie man leicht nachrechnet, gilt immer fiir %,> x, und alle y = 1,2,...n 


O.-—— Xa. 


= 0% 


bj xy 


a, — Xy a), — Xy 


bj—1 — Ky 


bye = fie, | 


| — %,| >| a 
1 1 


| bn — Xp | | bn — x» | 


Daraus erhalt man sofort durch Produktbildung 


und dies ist nur eine andere Schreibweise fiir die obige Behauptung. 


c) Man setzt 


eu - $$ 


Lb 
p—l 5 
—- 
xy = > | rene on Sal are Toe &, 


I 


y=) Gg a 
: v u 
mit 
' a (Oe @*) a2 
=i 
SIX ware oe =e 4 y (bo — dp) 
>» (bo — ae) ee 


(xy, entspricht im Falle n = 1 der hintere Neutralpunkt). 
Behauptung: Es gilt immer U(x,) <0, wobei das Gleichheitszeichen nur fiir x, = a, und 


x, = 6, eintritt. 


1 P. B. Fischer, Determinanten, S.115—116. Berlin u. Leipzig 1932. 
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Beweis: Der Integrand ist am linken Ende des Integrationsbereichs (bei a,) negativ, wechselt 
bei xy, das Vorzeichen und ist bei b, positiv. Es soll gezeigt werden, daB U(b,) = 0 ist. Damit 


gilt dann automatisch U(x,) < 0 fiir x, zwischen a, und b,. Da nach Setzung noch U(a,) = 0 
ist, ist damit die Behauptung bewiesen. 


Es gilt nach (29) und 2 
aye— Ey a en 
U(b,) = as | Sf) —e=e eee dé, 8 [let Lee 
ee : 
=e > 4 a] + ac <p Deal =o. 


d) Fir numerische Rechnungen kénnen manchmal die nachfolgenden Formeln von Nutzen 
sein, deren Beweis jedoch nicht gebracht werden soll. Ihre Herleitung ist vom Verfasser 1 ge- 
bracht oder lehnt sich eng an die dortigen Gedankenginge an: 


=o [mess Hest) dé, = D1.— (29) 


a Sf ES) § fe age : 1g Ale ee 
“= ) {V pe — fT" | ae nies resort: Co 
1 n Y ; n nee: x1 1 n 3 
£9 O—9) 3-9) +2 SO —ey, 
- - - 


Sie 4, (I/- Sipss cain | (32) 


De Bie: HEH) e451 Dm), (33) 
2487) Ups) a, == 3—1s, >'6,—a) F | 
ie; (34) 
#5[S@—a|—2 Srey, 


lot «# as Clas, Fea 
+f se(\/—He== | dé, =F fp — = a SH ie: 


y=1 


F4%|S@—«) | 552 U—9 7] >,—0)] 5 (G3) 


(Eingegangen am 7. April 1952.) 
Anschrift des Verfassers: Dr. Karl Nickel, Instituto Aerotecnico, Cordoba, Argentinien, 


* K. Nickel, Math. Z. 54 (1951), S. 89 ff, 
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Kritische Drehzahlen unter Torsion und Druck. 


Von H. Ziegler. 


1, Einleitung. Die kritischen Drehzahlen einer Welle sind von ihrer Beanspruchung abhangig. 
So ist z. B. bei der fliegenden Welle (Abb. 1), die am einen Ende in einem langen Lager drehbar, 
am anderen durch eine Scheibe besetzt und auf Druck D sowie auf Torsion W beansprucht ist, 
die kritische Winkelgeschwindigkeit @, eine Funktion von D und W. Bei der Untersuchung 
dieses Zusammenhangs ist R.Grammel!, der erstmals das Torsionsmoment beriicksichtigt hat, 
von der Annahme ausgegangen, dafs der Momentvektor des duBeren Kraftepaares bei der Ver- 
biegung der elastischen Linie axial gerichtet bleibe. Diese Voraussetzung ist bei einer in ihrer 
Ebene beaufschlagten Turbinenscheibe naiherungsweise erfillt, da sich hier die vom einzelnen 
Gasstrahl ausgeiibte Kraft bei der Deformation der Welle kaum verlagert. Ein solcher axialer 
Momentvektor ist nichtkonservativ?, und damit hangt zusammen, da sich unter gewissen 
Randbedingungen (z. B. im Fall von Abb.1) fiir W + 0 jede Drehzahl als kritisch herausstellt * 4. 


Axiale Momentvektoren kommen indessen auBerst selten vor. Bei der 
viel haufigeren schiefen oder axialen Beaufschlagung einer Turbinenscheibe | 
tritt neben einer Druckkraft ein Kraftepaar auf, dessen Momentvektor sich a! 
bei der Deformation schief stellt, da sich die ihn erzeugenden Einzelkrafte 
(Abb. 2), auch wenn sie nach Betrag und Richtung als konstant voraus- 
gesetzt werden diirfen, mit der Scheibe verlagern und damit keine festen 
Wirkungslinien besitzen. 


W 


Bei der Behandlung des analogen Knickproblems (das mit wm = 0 aus 
der Drehzahlaufgabe erhalten wird) hat der Verfasser® gezeigt, da im 
wesentlichen zwischen drei Fallen zu unterscheiden ist, die als pseudo-, 
quasi- und semitangential bezeichnet worden und konservativ sind, 
und W. T. Koiter® hat seither einen Fall behandelt, der als Sonderfall des 
pseudotangentialen aufgefaBt werden kann. Das pseudotangentiale 
Moment, das durch zwei, an einem starren Querbalken angreifende, kon- 
stante Krdafte realisiert wird, scheidet bei Turbinenscheiben aus; ebenso das 
Koitersche, in Kreuzgelenken auftretende Moment. Dagegen kommt bei 
doppelt beaufschlagten Scheiben (Abb. 2b) das durch ein einziges Krattepaar 
erzeugte quasitangentiale Moment und bei mehrfach bis gleichmafig 
beaufschlagten Scheiben (Abb. 2a) das semitangentiale Moment vor, W 
das einer gleichmafigen Verteilung von mehr als zwei gleichgroSen Umfangs- D 
kraften entspricht. Ist die elastische Linie durch x (z), y (z), ihre Neigung 
durch x'(z), y’(z), die Abszisse der Scheibe durch 2, und das an dieser an- 4), y. Fliegende Welle 
greifende Moment vor der Deformation durch W gegeben, so hat das quasi- Bat pine Der 
tangentiale Moment an der deformierten Welle unter der Voraussetzung, ee 
daB die beiden Einzelkrafte der y-Achse parallel seien, in erster Naherung 
die Komponenten® 


RX ay 


= W(x,,0,1). | (1.1) 


Im semitangentialen Fall dagegen wird es — unabhangig von der Orientierung der Koordinaten- 
achsen — durch 


1 a Ya 
w—w (Le, 2y,,1) (1.2) 
dargestellt. 


1 R.Grammel, Der EinfluB der Wellentorsion auf die kritische Drehzahi, Stodola-Festschrift S. 180, 
Zurich und Leipzig 1929; vgl. auch C. B. Biezeno und R.Grammel, Technische Dynamik, 5.796, Berlin 1939, 

2 H. Ziegler, Z. angew. Math. Mechan. 31 (1951), S. 266. 

3 H. Ziegler, Z. angew. Math. Physik 2 (1951), Sy 2607 

4 4, Troesch, Ing. Arch. 20 (1952), S. 258. 

5 H. Ziegler, Z. angew. Math. Physik 3 (1952), S. 96. 

6 W. T. Koiter, Proc. Koninklijke Akad. d. Wetenschappen, erscheint demnachst., 
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Nachdem sich! gezeigt hat, da® schon bei der Knickaufgabe die kritische Belastung von der 
Art des auBeren Momentes abhangt, sind auch beim Drehzahlproblem unter einem quasi- oder 
semitangentialen Moment (d.h. bei zwei- oder mehrfach, nicht in ihren Ebenen beaufschlagten 
Scheiben) verschiedene und zudem andere kritische Winkelgeschwindigkeiten zu erwarten als 
unter einem axialen. Zweck dieser Arbeit ist ihre Berechnung, und zwar unter Beschrankung 
auf einige besonders einfach gelagerte zylindrische Wellen mit zwei gleichen Biegesteifigkeiten, 
auf eine einzige Scheibe und auf Belastungen, die im Vergleich zur statischen Knicklast klein 
sind, 

Die Gesamtheit der bei der Untersuchung kritischer Drehzahlen auftretenden Fragen labt 
sich nur auf kinetischem Wege befriedigend beantworten. Der Rotor soll daher als Schwinger 
aufgefaBt und auf ein ortsfestes Koordinatensystem bezogen werden. Diese von der iiblichen 
abweichende Behandlungsweise wird gestatten, nebenbei einige Fragen zu klaren, die aufer- 
halb des engeren Problemkreises liegen. 


pe y 


L- 


>< 
Abb. 2. GleichmaBige Beaufschlagung Abb, 3, Elastische Linie im Abb. 4. Kraftespiel an der 
durch a mehr als zwei Krifte (semitangen- kritischen Zustand. Scheibe. 


tiales Moment), b zwei Krafte (quasi- 
tangentiales Moment). 


2. Die Bewegungsdifferentialgleichungen. Die in Abb.3 durch ihre elastische Linie dargestellte 
Welle sei als masselos, torsions- und drucksteif vorausgesetzt. Sie besitze zwei gleichgroBe 
Biegesteifigkeiten und sei beliebig, aber um ihre Achse reibungsfrei drehbar gelagert. In S sei 
eine starre Scheibe aufgesetzt, und zwar zunachst ohne jeden Zentrierungsfehler, so daB sie als 
vege Kreisel mit zur elastischen Linie tangentialer Figurenachse aufgefaft werden 
<ann. 

Die z-Achse des ortsfesten Koordinatensystems falle mit der geraden Stabachse zusammen. 
Da man sich auf kleine Deformationen beschranken darf, kann man die axialen Verschiebungen 
der Punkte P der elastischen Linie vernachlassigen und die Verschiebungen x(z), y(z) normal 
zur Achse samt den Neigungen x’(z), y’(z) als klein betrachten. Der tangentiale Einheitsvektor » 
in P hat dann in erster Naherung die Komponenten v, = x’, v1, = y’, v, = 1. 

Fiihrt man als zweites Bezugsystem dasjenige Hauptachsensystem der Scheibe ein, dessen 
¢-Achse den tangentialen Einheitsvektor ), enthalt und dessen €-Achse zur (x, z)-Ebene parallel 
ist, so ist seine 7-Achse in erster Naherung parallel zur (y, z)-Ebene. Mi®Bt man die Drehung der 
Scheibe im Hauptachsensystem durch den auf die §-Achse bezogenen Winkel q, so ist ihre Lage 
durch die Koordinaten x, y,, v;x, U;y, Y bestimmt, wahrend ihr Bewegungszustand durch die 
Translationsgeschwindigkeiten x,, y,; sowie die Winkelgeschwindigkeit u (u,, uy, u,) beschrieben 
werden kann. Da die zeitliche Anderung von ), durch ),; =U X 0, gegeben ist, gelten die Be- 
ziehungen 


Di, = Uy — yy Uy, Diy = — Uy vy, UL, (2.1) 


und diesen entnimmt man insbesondere, daB u, und uy von erster Ordnung klein sind. 


1 Siehe FuBnote S.5 auf S. 377. 
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Fir einen beliebigen Vektor a, dessen erste und zweite Kom 


ponente klein sind, gelten! in 
erster Naherung die Transformationen 


a; = a, — V4, 4, , G, =a, — yy @; , a, = a,. (2.2) 
Durch Anwendung derselben auf die Winkelgeschwindigkeit 1 erhalt man mit Riicksicht auf (2.1) 


ug = — Vy, Uy, = yx» ue =uUu,=. (2.3) 
Fiihrt man mit m die Masse, mit C das axiale und mit A das aquatoriale Tragheitsmoment der 


Scheibe fitr den Schwerpunkt S ein, so ist ihr Drall beziiglich S im Hauptachsensystem durch 
® (Au;, Au,, Cu,) und nach (2.2), (2.3) im ortsfesten System durch 


D, = D; + »,, D, =—Av,+Cor,,, | 
D, =D, +4, D, = <Arv,, + Cory, (2.4) 
D, =D, = Cy 


gegeben. 
Abb. 4 zeigt die in S angreifenden Krafte, namlich die aus der Axialkraft D und dem Mo- 
ment Y(W,, W,, W.) bestehende Belastung sowie die von der Welle herriihrende Dyname 


ES UP... Pests) , Mm (M,, My, M,). Formuliert man damit den Impuls- sowie den Drallsatz fiir 
die Scheibe, so erhalt man neben der Beziehung 


1D, = aa 1D) (2.5) 
die Differentialgleichungen 
Mx 4 = a P, r) } 
my, : ao P, 2 | 
— Avy, + C§(P 04. + @ x) = W,— M,, (2.6) 
A;, + C(p Viy + %1y) = Ws My, 
Cop = lee iin 2b 


3. Die unbelastete Welle. Ist die Scheibe unbelastet, so wird in (2.5) D=0 und in (2.6) 
W,. = W, = W, = 0. Bei stationarem Betrieb? ist ferner M, = 0, nach der letzten Gleichung (2.6) 
mithin g = 0,g =@ und gy = ot, so daB sich das System (2.6) auf 


Eee —Ai,,+Coi,,=—M,, | 
my,=— Py, Ai,, +Cwr,,=—M, | 


(3.1) 


reduziert. 

Die rechten Seiten von (3.1) lassen sich mit Hilfe der EinfluBzahlen* durch die Lagekoordi- 
naten der Scheibe ausdriicken. Ist namlich a die Durchbiegung an der Stelle z, unter einer 
hier wirkenden Einheitskraft, f die Neigung unter einem Einheitsmoment, y diejenige unter 
einer Einheitskraft und gleichzeitig die Durchbiegung unter einem Kinheitsmoment, so gelten 
die Beziehungen 


Ba = Wek cay, y1 =a Py —y M,, | (3.2) 
1,=y P, +p M,, yy 27) ye Mean 
Ihre Umkehrungen schreiben sich mit der positiven* Determinate 
N= apy (3.3) 
in der Form 
1 Eig = 
P, =F (B 1 — 9 12) » M,= Fly OX Vyy)» | (3.4) 
1 ] ; 
PpalP Yi — ¥ Vy) M, =F ( y %y +0 VY.) > 


1 Siehe die in FuBnote 3 von S. 377 zitierte tse oy: 

25C) BS Bt u. R. Grammel, Technische Dynamik, S. 777. ; oe 

Se TR ol be Kreisel, seine Theorie und seine Anwendungen, Bd. II, $.19. Berlin, Gottingen, 
Heidelberg 1950. 

ca Op B. Biezeno u. R. Grammel, Technische Dynamik, 5S. 89. 
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und damit 1aBt sich die Formanderungsenergie 
it 
U act Bact hye Y tects Moy, — 2, Vy) 
der Welle beilaufig auf die Gestalt 


eh Sl (xi + yf) —2 7 (%1 12 FI Vyy) +o (v2 + ty] (3.5) 


bringen. 
Setzt man (3.4) in (3.1) ein, so erhalt man das homogene Differentialgleichungssystem 


my yan 
: (3.6) 


—Av,,+Cor,, bi Avy =9, 


és . a 
Avia OD te Fy +5. = 0. 


Dieses System kann als Ausdruck des Impuls- und Drallsatzes fiir die nichtrotierende Scheibe 
unter den verallgemeinerten Kraften 


Qu. =—Co Viy ? Q.., = Co Vix (3.7) 


gedeutet werden; denn wenn man (3.5), (3.7) sowie die Bewegungsenergie 


. . Ay: : 
T= = (x? + yi) iat vio) (3.8) 
in die Lagrangeschen Gleichungen 
d aT oT aU 
dt 0° 0q ra 0q me Q, » (q = X1> Ya Vix» V1y) 


einsetzt, erhalt man in der Tat (3.6). Da U positiv definit und die Leistung des durch die ver- 
allgemeinerten Krafte (3.7) dargestellten Kreiselmomentes 


L= Qu ie M12 Si Q., Vy = 0 
ist, reprasentiert (3.6) die Bewegung eines konservativen, wenn auch gyroskopischen Systems 
mit positiv definiter potentieller Energie. Ein solches ist! aber stets stabil, und daraus folgt 
ohne Ausrechnung, da das allgemeinste Integral von (3.6) beschrankt und damit ungefahrlich ist. 
Um das Auftreten von kritischen Drehzahlen zu erklaren, miissen die Zentrierungsfehler 
der Scheibe beriicksichtigt werden. Dies geschieht iiblicherweise dadurch, da® man sich den 


Schwerpunkt der Scheibe in ihrer Mittelebene verschoben denkt, so daB er nicht mehr mit dem | 


DurchstoBpunkt der elastischen Linie (Abb. 3) zusammenfallt. Man kann statt dessen auch 
einige kleine Zusatzmassen einfithren, und wenn man diese auBerhalb der Scheibenmittelebene 
annimmt, beriicksichtigt man neben der Exzentrizitat der Scheibe auch ihren Aufkeilfehler, 
d. h. einen kleinen Richtungsunterschied zwischen ), und der Scheibenachse. Fir die Rechnung 
geniigt die Annahme einer einzigen Zusatzmasse 4 < m mit den Koordinaten 


§& =e, cosy, 1] =e, sil Y'; Ome; 


im Hauptachsensystem; von den endlichen GréBen e,, e, verkérpert dabei die erste fiir sich’ 


allein eine statische und die zweite in Verbindung mit der ersten eine dynamische Unwucht. 
Im ortsfesten System hesitzt die Zusatzmasse yu in erster Naherung die Koordinaten 


x, +e, cos, ¥, +e, sin@, a ae 


und ihre Riickwirkung auf die Scheibe besteht — vom Gewicht abgesehen — in den Tragheits- 
kraften 


T, = — (x, — e, sin py —e,cospg’), T,=—p(y, +e, cospy —e,singy*) (3.9) | 


' H. Ziegler, Ing.-Arch. 20 (1952), S. 49. 
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samt ihren, auf die Parallelen zu x, y, = durch S bezogenen Momenten 


Wr x = [ey (%, — ey sin PP — e, cos p Gg) , 
Wry = — pee (Vy TF €; Cos pY — e, sin yg) , | (3.10) 
W,, = we, (x, sin yp —¥, cos p — e,f) . 


Erginzt man die rechten Seiten von (2.6) durch die Terme (3.9), (3.10), und beschrankt 
man sich wieder auf die unbelastete Welle, so geht die fiinfte Beziehung in 


(C + wel) & = pre; (x, sin p — jy, cos gy) — M, (3.11) 


tiber. Bei gut ausgewuchteten Wellen diirfen die mit tu behafteten Glieder ohne Bedenken 


_unterdriickt werden, so daB man im Falle stationaren Betriebes wieder y = wt erhalt. Damit 


gewinnt man in der gleichen Naherung aus den vier ersten Relationen (2.6) statt (3.6) das System 


Se » 
mx, 4 A %1 | Ye = ey & cos wt, 


= B y : 
my, +9, ey =pe,w*sinat, 


fs cet: x ; (3.12) 
A vy5 + Cor, Rey tg Cay = = It &y €y'@" cos Mt, 
Ai : y x : 
M2 + Corny, —F% +7 Mix = Mey eo sin wt, j 


und dieses geht, wenn man die Lagekoordinaten der Scheibe mit 
Cee oe eae Vie Uy — (3.13) 
komplex zusammenfaBt, in 


mr se aah = [1 &; @* exp (tt), 
(3.14) 


= —tipe, e, w* exp (tt) | 


ce aes, : y l a 
Av—iCowv icieages, 
iiber. 

Da die linken Seiten in (3.6) und (3.12) iibereinstimmen, ist die allgemeinste Liésung des 
reduzierten Systems (3.14) beschrankt. Ein partikulares Integral des inhomogenen Systems 
kann in der Form 


r = aexp (1wt), v = bexp (iwt) 
angesetzt werden, und die Integrationskonstanten a, b bestimmen sich aus (3.14) mit der Ab- 
kiirzung 
FOR) Aaa || (C= A) |, (3.15) 
A ‘lege t A? ; 
zu 


__ ey wo? | o Let ee. 
TERT E 1 (C — A)w U | ; 
wewetly —e,(B— mo] 
De Ft) é Ley t mo*)|. 
Das Integral ist nur fiir f #0 brauchbar, dann aber auch beschrinkt. Fir f > 0 wachsen in- 


dessen a und b im allgemeinen iiber alle Grenzen, so da die kritischen Winkelgeschwindigkeiten 
durch f = 0 gegeben sind. Dafiir kann man nach (3.15) und (3.3) auch 


g (w?) [a m — B (C — A)] Am(C— A)=0 (3.16) 


schreiben und erhalt, wenn man die Kreiselwirkung zunichst vernachlassigt, also C= A = 0 
setzt, die bekannte Beziehung 


at 


peas (3.17) 


mo 


fiir die kritischen Winkelgeschwindigkeiten w) und —w. 
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Da nicht nur (3.16), sondern auch die in den folgenden Abschnitten herzuleitenden Be- 
stimmungsgleichungen fiir die kritischen Winkelgeschwindigkeiten in @ quadratisch sind, wollen 
wir uns grundsitzlich auf positive Winkelgeschwindigkeiten beschranken. 


4, Der Einflu8 der Kreiselwirkung. Davon abgeschen, daB die Gleichung (3.16) beide Trag- 
heitsmomente enthalt und noch nicht mit C = 2 A fiir diinne Scheiben spezialisiert ist, stellt 
sie die Grammelsche Beziehung! fiir die kritischen Winkelgeschwindigkeiten @, des Gleich- 
laufs dar. Sie fiihrt auf 


1 1 tym—p(C— A) 4/4 [am—f (C— AP + 4m (C— A) (4.1) 


und liefert fiir abgeplattete Rotoren (C > A), da hier das negative Wurzelvorzeichen ein 
imaginares «, ergeben wiirde, nur eine kritische Winkelgeschwindigkeit @,, die stets gréBer 
ist als Wp. 
Setzt man, durch die Grammelsche Darstellung geleitet, 
C—A=mF, (4.2) 
so erhalt man aus (4.1) und (3.17) nach kurzer Zwischenrechnung fiir das Verhaltnis der kri- 
tischen Winkelgeschwindigkeiten mit und ohne Kreiselwirkung 


Ok 1 
Wo V1 —p (hk?) 


(4.3) 


wobei 


1 1 B 2 y? 
pi) =La+4mzyt(i+ful_te (4.4) 
und fiir C > A allein das obere Wurzelzeichen giiltig ist. 


Abb. 5 zeigt den Anstieg von (4.3) mit zunehmendem k*, wobei der Auswertung von (4.4) 
speziell eine fliegende Welle mit Endscheibe zugrunde gelegt und der ,,Tragheitsradius k 
auf die Lange | der Welle bezogen ist. Von der Grammelschen Abb. 13 unterscheidet sich Abb. 5 
grundsatzlich nur dadurch, daf statt k als Abszisse k?/Il? verwendet ist; im iibrigen gibt sie, 
da (4.2) noch beide Tragheitsmomente enthalt, das Verhaltnis w,/q, fiir beliebige Abplattungen. 

In der Praxis kommen auch Rotoren vor, die als verlangerte Kreisel (C < A) zu gelten | 
haben, und fiir diese miissen beide Wurzelvorzeichen zugelassen werden. Da der Radikand in 
(4.1), der mit (3.3) auf die Gestalt 


= [a m +B (C — A)? — 7% m(C— A) 


gebracht werden kann, fiir C << A positiv, nach (4.1) aber kleiner als [a m — B (C — A)}?/4 ist, 
gibt es, wie bereits O. Féppl+ gezeigt hat, bei verlingerten Rotoren zwei kritische Winkel- 
geschwindigkeiten des Gleichlaufs. Sie lassen sich auch wieder mit Hilfe von (4.3) und (4.4) 
ermitteln, wobei jetzt aber beide Wurzelvorzeichen zu beriicksichtigen sind. Da (4.4) fiir k? < 0 
eine positive und eine negative Wurzel besitzt, ist die eine der beiden kritischen Winkelge- 
schwindigkeiten @, kleiner, die andere gréBer als a. 

In Abb. 5 ist (4.3) auch fiir negative k? ausgewertet; sie liefert also die kritischen Winkel- 
geschwindigkeiten des Gleichlaufs auch fiir verlangerte Rotoren. Man iiberzeugt sich an Hand 
von (4.3) und (4.4) leicht davon, daB w, = 0 eine Asymptote der tiefer liegenden Kurve, k = 0 
eine solche der héher liegenden ist, und daf beide Kurven die gemeinsame Asymptote w;/a@, = 
(« / A)!’ besitzen, die fiir die fliegende Welle speziell in w/w, = 2 itbergeht. 

Die kritischen Winkelgeschwindigkeiten des Gegenlaufs, die der Beziehung 


og = lm +B(C +A tY/ Flam +h (C + AP—Am(C + A) (4.5) 


geniigen, werden vom System (3.12) nicht geliefert, und daraus folgt, daB sie sich weder mit 
statischen noch mit dynamischen Unwuchten erklaren lassen. 

1 R. Grammel, Der Kreisel, S. 16. 

® O. Féppl, Schweiz. Bau-Ztg. 65 (1947), S. 259. 


schrieben, in denen die Stérungs- 
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Gleich- und Gegenlauf werden iiblicherweise als Gleichgewichtsformen des aus der Welle 
und der ideal ausgewuchteten Scheibe bestehenden Systems fiir einen mit der Winkelgeschwindig- 
keit ® bzw. —q rotierenden Beobachter definiert. Als solche stellen sie Lésungen des redu- 
zierten Systems (3.14) von der Form r = a exp (twt), v = b exp (i wt) baw. r = aexp (—iot), 
v = b exp(—iwot) dar und haben somit fiir den ruhenden Beobachter als spezielle Kigen- 


schwingungen zu gelten. Mit der ersten ist die von einer Unwucht herrithrende, mit exp (i wt) 
proportionale Stérung in Resonanz, 


mit der zweiten dagegen nicht. Schon 
A. Stodola, der erstmals den Gegen- 
lauf theoretisch und experimentell 
untersuchte!, hat aber auf die Mig- 
lichkeit der Anregung durch peri- 
odische Krafte hingewiesen, die nicht | 
von Unwuchten, sondern vom Mit- 
schwingen des Tragers herriihren. 

Nimmt man, um sich auf den 
einfachsten Fall zu_beschranken, 
an, daB die Welle in oder mit ihren 
Lagern eine zusdtzliche, zu ihrer 
Achse normale Translationsschwin- 
gung mit der Kreisfrequenz w aus- 
fiihre, so wird ihre Bewegung relativ 
zum mitschwingenden System durch 
die Differentialgleichungen (3.12) be- 


: : Abb. 5. EinfluB der Kreiselwirkung auf die kritische Winkelgeschwindigkeit; 
glieder der Reihe nach durch [k? = (C— A)/m im Gleichlauf; k? = —(C + A)/m im Gegenlauf]. 
b-) 
oo fo) 
> (Acncosnot + B,,sinnot), > (Ay, cosnot + B,, sin n ot) 
n=1 n= 


sowie zwei Nullen zu ersetzen sind. Der Ubergang zur komplexen Schreibweise (3.13) ergibt 
statt (3.14) 


miter—Yy= 3S Cyexp(nios), 


(4.6) 
Ab—iCwmi—Fr+7v=0, | 
wobei die C, komplexe Konstanten sind. Mit dem Ansatz 
r= aexp (niat), v = bexp (ni ot) 
kommt jetzt statt (3.15) 
fr (@?) = (4 —nt mo?) 5 4+n(C—1n A) o = ele 0, le oes (4.7) 


+ 1 zunichst auf die Beriicksichti- 
gung der Grundschwingung beschrankt, erhalt man neben (4.1) die Bedingung (4.5). Diese unter- 
scheidet sich von (4.1) nur im Vorzeichen von C und liefert die beiden kritischen Winkelge- 
schwindigkeiten des Gegenlaufs, die wieder aus (4.3) und (4.4) und damit auch aus Abb. 5 er- 
halten werden, wenn man statt (4.2) jetzt 

—(C+A)=mFP (4.8) 
setzt. Mit n = +2,...gewinnt man sodann aus (4.7) unendlich viele weitere kritische Winkel- 
geschwindigkeiten, die den héheren Harmonischen der Stérung entsprechen und damit auch ohne 
Kreiselwirkung vorhanden sind. 

Das statische Verfahren, das als einzige und (bei Vernachlassigung der Dampfung) gleich- 
wertige kritische Zustande Gleich- und Gegenlauf liefert, gibt demnach die wirklichen Ver- 
haltnisse verzerrt wieder. Tatsichlich ist ohne du®ere Erschiitterungen nur der Gleichlauf 
kritisch, und zwar zufolge Resonanzgefahr mit den von (statischen oder dynamischen) Unwuchten 


und wenn man sich in der Entwicklung der Stérung mit n = + 


1 A, Stodola, Z. ges. Turbinenwes, 15 (1918), S. 253, 
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herrithrenden Trigheitskraften. Sind auBere Erschiitterungen vorhanden, so besteht auch Re- 
sonanzgefahr mit diesen, und damit werden neben dem Gleichlauf (und zwar auch ohne Kreisel- 
wirkung) unendlich viele weitere Drehzustande kritisch, unter denen (mit Riicksicht darauf, 
daB die Entwicklung der Stérung praktisch rasch konvergiert) der Gegenlauf am starksten in 
Erscheinung tritt. 

Da®B der Gegenlauf tatsachlich nicht nur ungefahrlich, sondern sogar nur selten nachzuweisen 
ist, wird iiblicherweise mit der inneren Dampfung der Welle erklart, die sich bei diesem Zustand 


dauernd deformiert. Ein weiterer Grund dafiir besteht aber zweifellos darin, daB die Erschit- — 


terungen der Welle praktisch meist gering sind. Beide Tatsachen (samt der raschen Konvergenz 
der Entwicklung) erklaren schlieBlich auch, daB die zu n= —+2,... gehdrenden kritischen 
Winkelgeschwindigkeiten in der Praxis nicht beobachtet werden. 


5. Der Einfluf& der Belastung. Ist die Scheibe nach Abb. 4 durch eine axiale Kraft D sowie ein 
Moment %& belastet, die Welle mithin auf Druck und Torsion beansprucht, so kann (2.6) nicht 
mehr durch (3.1) ersetzt werden. Vernachlassigt man, um den EinfluB von D und Y% auf die 
kritischen Drehzahlen zu erhalten, die Kreiselwirkung, so folgt aus der dritten und vierten Be- 
ziehung (2.6) 


Meas, Mi, Ta (5.1) | 


Im stationaren Betrieb ist ferner 
M, == We, (5.2) 
nach der letzten Gleichung (2.6) also wieder g = wt, und damit reduziert sich das System auf 
DV ee be OF my, -+ P= 0. (5.3) 


Auch hier kénnen die rechten Seiten mit Hilfe von EinfluBzahlen durch die Lagekoordinaten 
der Scheibe ausgedriickt werden. Man fihrt zu diesem Zweck mit o;;, (t, k = x, y) die in der 
i-Richtung gemessene Durchbiegung an der Stelle z, unter dem EinfluB einer hier in der 
k-Richtung wirkenden Einheitskraft, der Druckkraft D und des Momentes Y ein, wobei damit 
zurechnen ist, da unter der durch YW gegebenen Belastung Oxy AO, dy, AO und ayy Foxx 
ist. Mit diesen EinfluBzahlen gilt 


X= xx Pp ey Pe Y1 = My x di Recs Ce (5.4) 


und wenn man noch die Determinante 
N= Oly, gy One (5.5) 
einfiihrt, die+ wie (3.3) stets positiv ist, erhalt man hieraus durch Umkehrung 
1 il 
IPs =F (yy Mey Va) Py =F (—4y« X, + Onn 3) « (5.6) 


Durch Einsetzen von (5.6) in (5.3) gewinnt man die Differentialgleichungen 


| (5.7) 


my Wile als ae 


die als Ausdruck des Jmpulssatzes fiir das aus Scheibe und Welle bestehende System aufgefaBt | 


werden kénnen. 


Die EinfluBzahlen hingen von der Lagerung der Welle, bei kurzen Lagern, die eine Schief- 
stellung der elastischen Linie zulassen, aber auch davon ab, welcher Art die hier angreifenden 
Momente sind. Wirkt an einem kurzen Lager beispielsweise ein — durch ein Kreuzgelenk 
realisiertes — quasitangentiales Moment, so sind infolge der Drehung des Gelenkes die o;; sogar 
periodische Funktionen der Zeit. Um in dieser ersten Untersuchung zeitabhangige Einflu8- 
zahlen auszuschlieBen und zugleich die Fallunterscheidungen auf ein ertragliches MaB zu re- 


reduzieren, beschrinken wir uns im folgenden auf lange Lager, die fiir die elastische Linie als | 


Kinspannungen gelten kénnen. Unter dieser Voraussetzung sind nicht nur die a;, konstant, 
sondern samtliche an der Welle angreifenden Krifte bei quasi- wie semitangentialer Belastung 


' C. B. Biezeno u. R. Grammel, Technische Dynamik, S. 89, 
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der Scheibe konservativ. Es muB also 


aPy  OPs 
Ox, Oy, 
sein, und daraus folgt nach (5.6) 
Oyx = Oxy » (5.8) 


a. :h; die Richtigkeit des Maxwellschen Satzes auch in diesem Falle. Da ferner die potentielle 
Energie 


= 94 ry Or ae ey X11 aie eiVi) (5.9) 


des ganzen Systems fiir Belastungen, die unter der kritischen bleiben, positiv definit, das System 
also stabil ist, bleibt die allgemeinste Lésung von (5.7), wie man ibrigens leicht durch Integration 
bestatigt, beschrankt. 

La8t man jetzt wieder Unwuchten zu, wobei man sich bei Vernachlassigung der Kreisel- 


wirkung konsequenterweise auf statische zu beschranken hat, so kommt statt (5.7) das inhomo- 
gene System 


ayy Oxy 
tease aed = mew" cos wt, | 


(5.10) 
my, 79 ® 9g ets = wew* sin wt | 
Ayo a 1 ‘ 


Dieses besitzt ein partikulares Integral 
x, =a,coswt +b, sin wt, 


¥, = a,cos wt + by sin wt, 


wobei 
_ bew One ) _ pew? Syed 2 
a, ot 7 mo?) , by = F(a) \ A 1 = \\ 3 
ew a&xy 
b. ey f(a?) A 
und 


{(oe)= a —m o) e —m 0) = Say (5.11) 


ist. Man zieht hieraus wie schon in Abschnitt 3 den Schlu8, daB die kritischen Winkelgeschwin- 
digkeiten durch f=0 gegeben sind. 

Da das reduzierte System (5.7) im allgemeinen keine Lisung der Form x, = acos wt, 
y,=asin wt besitzt, kann die Gleichung f=0 hier nicht mehr als Bedingung fiir die Existenz 
einer nichttrivialen Gleichgewichtslage im mitrotierenden Koordinatensystem gedeutet werden. 
Damit ist beilaufig festgestellt, daB das Grammelsche Aquivalenzprinzip+ hier ebensowenig gilt 
wie etwa bei der unbelasteten Welle mit zwei verschiedenen Biegesteifigkeiten”. 

Ist das Moment Y semitangential, so kann es, wie ein Blick auf (1.2) zeigt, vektoriell 
angeschrieben werden. Da dann keine Achse vor der anderen ausgezeichnet ist, unterscheiden 
sich die Verschiebungsvektoren (a, , %yx) und (Ay, &yy), die von Einheitskraften in der x- bzw. 
y-Richtung erzeugt werden, nur durch eine Drehung um 7/2. Hieraus folgt aber ayy= ax sowie 


Oy, = —d,y und mit Riicksicht auf (5.8) 


Onn = yy =O, Oyx = Oxy =0, 1A ene (5.12) 


Der Punkt S auf der elastischen Linie (Abb. 4) verschiebt sich demnach unter einer beliebig 
gerichteten horizontalen Kinheitskraft in deren Richtung, und zwar um die Strecke x, die jetzt 
freilich nicht mehr mit der Einflu®zahl « der unbelasteten Welle iibereinstimmt, sondern mit 
Riicksicht auf die Belastung durch D und  gréfer ist. 

1 Vgl. C. B. Biezeno und R. Grammel, Technische Dynamik, 5. 782. 

2 Ebenda, S. 784- 
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Mit (5.11) und (5.12) geht die Bedingung f = 0 in 


(yd es eee (5.13) 


m 


QI 


iiber. 
Die lang gelagerte, mit einer einzigen Scheibe besetzte Welle besitzt also ohne 


Kreiselwirkung bei semitangentialer Belastung eine einzige kritische Winkel- 
geschwindigkeit @,. Diese ist kleiner als wy und wird mit (5.13) erhalten, wobei 
die EinfluBzahl unter Beriicksichtigung der Belastung D, % zu ermitteln ist. 

Im Falle eines quasitangentialen Momentes (1.1) ist das Problem nicht mehr rotations- 
symmetrisch. Die Beziehungen (5.12) gelten dann im allgemeinen nicht, und die Bedingung 
f =0 geht mit (5.11) in die Gleichung 


8(w*) = —— m (x2 + dyy) =: +m/4=0 (5.14) 


iiber. Diese hat denselben Aufbau wie (3.16) beim verlangerten Rotor, und die fiir diesen durch- 
gefiihrten Uberlegungen zeigen, da die beiden Wurzeln 


er Pea Ce er (5.15) 


) i. Z 


positiv sind. 

Die lang gelagerte, mit einer einzigen Scheibe besetzte Welle besitzt mithin 
ohne Kreiselwirkung bei quasitangentialer Belastung zwei kritische Winkel- 
geschwindigkeiten, die mit (5.15) erhalten werden. Sie kénnen freilich — min- 
destens inerster Naherung — zusammenfallen. 


6. Semitangentiale Belastung. Im semitangentialen Fall (mehrfach bis gleichmabig beauf- 
schlagte Scheibe) ist die kritische Winkelgeschwindigkeit durch (5.13) gegeben und ihre Be- 
stimmung damit auf die Ermittlung der Einflu8zahl & unter der Belastung D, ¥$ zuriickgefiihrt. 

Als Beispiel greifen wir die fliegende Welle (Abb. 1) heraus, die in 
Abb. 6 samt ihrer Belastung durch die Axialkraft D, das semitangen- 
tiale Moment (1.2) sowie eine horizontale Kraft §§ mit den Kom- 
ponenten P,, P, wiedergegeben ist. Die Differentialgleichungen 
ihrer elastischen Linie lauten 


EIx” = M,, Ely” =— M,, (6.1) 
wenn mit EJ die Biegesteifigkeit und mit M:, M, die Komponenten 
des Biegemomentes beziiglich der Hauptachsen (Abb. 3) des Schnittes 


P bezeichnet werden. Liest man diese Komponenten in Abb. 6 ab, 
so erhalt man 


tt , 1 / 
EIs" = —Wy +3 Wy, + Dy — 2) +P0—2),) 
(6.2) 
lad {2 1 ff 
Ely” = Ws’ —1 Wx, + Diy, —y) +P, 12), | 
Abb. 6. Fliegende Welle unter einer Oder, nach z abgeleitet, 
Druckkraft und einem semitangen- : 
ialen Momen eaufschlagung in Ee ing f = 
’ mehr ae Peay ; Eig a Wy a5 Dx ats FE, it) 0, | (6.3) 


Ely'” —Wx'' + Dy’ + P, =0. J 


Fiihrt man mit » (v,=2’, vy=y’, v;=1) wieder den tangentialen Einheitsvektor sowie die 
Abkiirzungen 


Was DE Px P 
EI=™” ET =» Eq = P» at =p, (6.4) 


ein, so gewinnt man durch komplexe Zusammenfassung 
Vt bya Px + 1 py = p (6.5) 
an Stelle von (6.3) eine einzige Differentialgleichung 


vw’ —iwv +bv+p=0, (6.6) 


XX. Band 1952. Ziegler: Kritische Drehzahlen unter Torsion und Druck. 387 


deren allgemeinste Lésung unter der Voraussetzung, daB b £ 0 ist, mit 


1 ee 
Ay,2= 5 (w +0), a= fut +46 (6.7) 


und den Integrationskonstanten A, B durch 


v = Aexp (iA,2) + Bexp (iA,z) —? (6.8) 
gegeben ist. 
Die Randbedingungen lauten 


0, =U, vi — fiw = 0. (6.9) 


Die erste davon driickt aus, da die elastische Linie am unteren Ende vertikal ist; die zweite 
folgt aus der Forderung, da neben (6.3) — insbesondere am oberen Ende — auch die Differential- 
gleichungen (6.2) erfiillt sein miissen. Die Bestimmung der Integrationskonstanten mit (6.9) 


fiihrt (6.8) in 


Pp exp (tA, 1) — w/o F exp (iA,1) + w/o ahh 
b = A EP ET Ge i ree GA) + SAG Oa) 


iiber. Die Verschiebung des oberen Endes wird hieraus mit 


= 


I 
r= % +1y, aH vdz (6.11) 


erhalten und hat, wie man nach einigen Zwischenrechnungen feststellt, die Form 


(6.12) 


Da in dieser Gleichung nur r, und p komplex sind, hat die Verschiebung erwartungsgem4B die 

Richtung der Kraft S$, und die gesuchte Einflu8zahl ist, wie ein Blick auf die beiden letzten 
Beziehungen (6.4) zeigt, 

_ wl 

sin — 

] w? + o? ol it 2 
Ed toh <2) 2") coro 
2 


Unter der Annahme, daf die Druckkraft D im Vergleich zur Knickkraft klein sei, und dah 
entsprechendes auch fiir das Moment W gelte, kann & in eine nach Potenzen von b und wv? fort- 
schreitende Reihe 


a= Bin (6.13) 


Ah Ey we rr FT (6.14) 
3 EI 5 80 ies Z 
entwickelt werden, und dafiir kann man unter Beriicksichtigung von (6.4) sowie der Tatsache, 
da® der Faktor vor der Klammer die EinfluBzahl « der unbelasteten Welle darstellt, auch 
a rAee AS tha LA 
eg aT (a) 
schreiben. Man hat also nach (3.17) und (5.13) fir das Verhaltnis der kritischen Winkel- 
geschwindigkeiten mit und ohne Belastung in erster Naherung 


(6.15) 


O,, Diz wi\? 
peel 6.1 
rican. ki ey Aes : (6.16) 
wobei 
1 | 
= TAL k,= [60 (6.17) 


ist. 
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Die Korrekturfaktoren k, und ky, die den Einflu8 der Druckkraft D und des semitangentialen 
Momentes W auf das Verhaltnis ©,/«) angeben, kénnen fiir andere Faille mit langen Lagern analog 
bestimmt werden und sind in Tabelle 1 fiir die drei haufigsten Falle zusammengestellt'. 


7. Quasitangentiale Belastung. Bei quasitangentialer Belastung (doppelt beaufschlagte 
Scheibe) sind die kritischen Winkelgeschwindigkeiten durch (5.15) gegeben und ihre Bestimmung 
damit auf die Ermittlung der EinfluBzahlen o,., yy, &vy Zuriickgefithrt. Fiir kleine Werte von 
D und W kann man sich auch hier auf eine erste Naherung beschrinken, d.h. auf eine ab- 
gebrochene Reihenentwicklung analog (6.16) ausgehen. Da diese fiir k, den bereits bekannten 
Wert liefern mu, darf man die Druckkraft von Anfang an weglassen. 

Im Falle der fliegenden, gema® Abb.7 belasteten Welle gelten 
nach (6.1) fiir die elastische Linie die Differentialgleichungen 


EIx«"” =—Wy’ a oP a(n z Jer) 
A aie Wx! —Wx,+P,(l—2).) 
Durch Ableiten nach z erhalt man hieraus wieder die Beziehungen 
(6.3), jetzt freilich mit D= 0, und mit den Abkiirzungen (6.4), (6.5) 


(7.1) 


zh VP folgt die komplexe Differentialgleichung 
le 
v’ —iwv’ + p=0 (7.2) 
mit der allgemeinsten Lésung 
4 = A exp (iwz) + B—iF az. (7.3) 


Die erste Randbedingung (6.9) gilt auch hier und liefert 
x v = A [exp (iw2) —1]—i~ a; (7.4) 


Abb. 7. Fliegende Welle unter einem 
quasitangentialen Moment 
(Beaufschlagung in zwei Punkten). 


die zweite ist jetzt durch die beiden Forderungen 

Vx, + wv,, = 0, Uy, = 0 (7.5) 
zu ersetzen, und diese machen die Riickkehr zur reellen Schreibweise nétig. Setzt man zu diesem 
Zweck A= A, +iA,, so zerfallt (7.4) in 


v, = A, (coswz—1)—A 


y SIN WZ Se ee 
w 


s (7.6) 
vy = A, sin wz + A, (cos wz — 1) =e Bes 
und die Integrationskonstanten A,, A, bestimmen sich aus (7.5) zu 
— 1 Py : Px 
ae Sale Ip.) sin wi +2], | 
} (Tea) 
Aaa ees | 
Ca ae 
Die Verschiebung des oberen Endes ist durch 
1 
; ) 
t= od = A,(#! — 1) ce ange V2 
1 / E aN Giga 7 + op 
‘ \ (7.8) 
si | he ee wl—] a pxP? 
Sil J Of & x Ai = A y ip l Day j 


gegeben. Fuhrt man die Integrationskonstanten (7.7) ein und setzt man gleichzeitig p, = 1/EI, 
Py = 9, so erhalt man nach (6.4) die Verschiebungen unter einer Einheitskraft in der x-Richtung, 


d.h. die EinfluBzahlen 
gen?! 1 2 l 2 
ae ayl(4s +2) +e wt <3 ( —1) ; 


1 l 1 1 P 
Oy ¢ = eves 
sf | w ge ws i wl 1 *| ed 


1 Die Werte fiir die beidseitig gelagerte Welle verdanke ich Herrn Chr. Wehrli, dipl. Math. ETH. 


(7.9) 


| 
| 


i 
J 
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und ganz entsprechend kommt mit pat: Py = 1/EI neben On y= Oy x 
Be aii! l 
Oyy = By piel. ; (7.10) 
Die Entwicklung von (7.9) und (7.10) ergibt 
ak 3 8 2 
Onn ae 50 we...) > ty eet eee, | > 
. (7.11) 
Ox y = lt a wel, . ; ‘ | 


Schreibt man den Randikanden von (5.15) mit Hilfe von (5.5) in der Form (6, , == yi) + 4 Oey ; 
so stellt man an Hand von (7.11) leicht fest, daB® sich die Wurzel in der ersten Naherung, auf die 
wir uns nach dem Verzicht auf die Mitnahme von D notwendig beschraénken miissen, auf 2 |ox | 
und der erste Klammerterm auf 


218 
ee (7.12) 
reduziert. Man hat also in dieser Naherung 
1 
op = M(% + Oey)» (7.13) 
zufolge (3.17) 
ss wee) 
wale (7.14) 
und nach (7.11) sowie (7.12) 
el eee] (7.15) 
Wo 16 , ‘ 


Nimmt man die Druckkraft schlieSlich wieder mit, so kommt gema (6.4) und (6.17) endgiiltig 


Or, DIP Wi 
wobei 
il 3 ‘ 
= =o | = Fs (7.17) 
ist. 


In Tabelle 1 ist neben k, und k, auch der Korrekturfaktor k, fiir das quasitangentiale Moment 
eingetragen?, 

8. SchluBbemerkungen. Durch Zusammenfassung von (6.16) und (7.16) erhalt man die bei 
langer Lagerung stets giiltige Beziehung 


Oy, DI Wil Wi\2 
1k p thay bs, (8.1) 


Mo 


deren Beiwerte k,, ky, k,in Tabelle 1 fiir drei einfache Falle zusammengestellt sind. 

Im Fall II ist im Gegensatz zu I und III die Scheibe selbst unbelastet, und da die Lager 
lang vorausgesetzt sind, braucht bei den Endmomenten nicht mehr zwischen dem semi- und 
dem quasitangentialen Kraftangriff unterschieden zu werden. Es ist hier ky = 0, und fiir ky, k, 
werden erwartungsgem48 die von R.Grammel? angegebenen Werte erhalten. 

Auch in den iibrigen Fallen ist bei semitangentialer Belastung (Beaufschlagung der 
Scheibe in mehr als zwei Punkten) k,=0; sowohl die Druckkraft D wie das Moment W haben 
(wie schon im Fall II) ein Absinken der kritischen Winkelgeschwindigkeit zur Folge. 

Bei quasitangentialer Belastung (Beaufschlagung in zwei Punkten) verhalt sich Fall III 
(in der beobachteten Naherung) genau so wie bei semitangentialer. Wie in Fall II, wird auch 
hier @, durch D und W gesenkt. Im Fall I dagegen ist k, +0 (so daB der dritte Term in (8.1) 


1 Auch hier verdanke ich die Zahlenwerte fiir die heidseitig gelagerte Welle Herrn Chr. Wehrli, dipl. 


Math. ETH. 
2 C. B. Biezeno und R. Grammel, a. a. O. S. 800. 


Ingenieur-Archiv 


390 Ziegler: Kritische Drehzahlen unter Torsion und Druck. 


unterdriickt werden mu). Hier wird «, durch D erniedrigt und durch W in zwei kritische 
Winkelgeschwindigkeiten aufgespalten, von denen die eine kleiner, die andere gréBer ist als 
diejenige unter Druck allein. . 

Aus (8.1) und den Zahlenwerten in Tabelle 1 geht hervor, dafB das quasitangentiale Moment, 
sofern es sich iberhaupt vom semitangentialen unterscheidet, die kritische Winkelgeschwin- 
digkeit eher starker modifiziert als dieses. Die Beriicksichtigung von W drangt sich daher vor 
allem bei doppelter Beaufschlagung und hier in erster Linie bei der fliegenden Welle auf. 


Tabelle 1. EinflufB von Druck und Torsion auf die kritische Winkelgeschwindigkeit 
; (ky, ke, kg Korrekturfaktoren in (8.1)]. 


Fal! I I i 


Lagerung D 


Belastung 


Auch in den Fallen I und III stimmen die Werte von k, mit den Grammelschen iiberein; 
diejenigen von k, unterscheiden sich aber davon, und ky stellt eine Korrektur dar, welche den 
mit W? proportionalen Term vidllig iiberschattet. 

Bei kurzen Lagern sind (in den Fallen II und III) andere Ergebnisse zu erwarten. Hier 
ist es fiir die Bestimmung der HinfluBzahlen wesentlich, welcher Art die am Ende der Welle 
angreifenden Momente sind. Vom praktischen Standpunkt dirfte vor allem das rotierende 
quasitangentiale Moment von Interesse sein, wie es in Kreuzgelenken auftritt. Wie schon in 
Abschnitt 5 bemerkt wurde, werden aber in diesem Falle die «;; im raumfesten Koordinaten- 
system periodische Funktionen der Zeit, und die Integration des Systems (5.10), das damit 
quasilinear wird, stellt ein véllig neues Problem dar. 


(Eingegangen am 20. Mai 1952.) 
Anschrift des Verfassers: Prof. Dr. H. Ziegler, Riischlikon bei Ziirich, Weiherweg 6. 
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Uber die Anwendung der Mellin-Transformation zur Lisung einer 


Aufgabe der Plattenbiegung. 
Von S. Woinowsky-Krieger. 


1. Einleitung. Zu den noch wenig geklarten Fragen der Theorie der Plattenbiegung gehért 
die Frage nach der Spannungsverteilung in der Ecke einer diinnen Platte bei einem beliebigen 
Offnungswinkel und bei beliebigen Grenzbedingungen lings der beiden anschlieBenden Platten- 
rander. Sind diese Rander unendlich ausgedehnt, so lat die Aufgabe auch eine strenge Lésung 
zu, die sich bequem mit Hilfe der Mellinschen Umkehrformel aufbauen laBt. 


Als Anwendungsbeispiel wahlen wir eine unendlich erstreckte keilférmige Kragplatte mit 
beliebigem Offnungswinkel. Da es sich vor allem darum handelt, die Brauchbarkeit der ge- 
wahlten Transformationsmethode zur Liésung von Plattenaufgaben der vorliegenden Art zu 
zeigen, so mégen im Folgenden nur einige spezielle Belastungsfille behandelt werden. Es ergibt 


sich dann von selbst der Weg zu einer Verallgemeinerung der Lisung fiir eine beliebig verteilte 
Belastung. 


2. Die durch eine Einzelkraft am Rande belastete Kragplatte. Die in Abb. 1 dargestellte, langs 
des Randes #=0 eingespannte Platte trage an der Stelle r—ry, ®—a ihres freien Randes eine 
Punktlast P. Die elastische Flache der Platte w(r, 0) geniigt dann iiberall, auBer in dem Punkte 
(ro, x), der Differentialgleichung 


engespanntrer Rand, 


AA w= 0. Gy 2 


Wahlt man ferner eine nichtnegative Zahl s als Parameter, 
so kann man die Lésung der Gleichung (1) in der Form 


W(s) =r *O(s) (2) 
mit 
-O(s) = A(s) sin sP + B(s) cos s# + C(s) sin (s ++ 2) 9} (3) Abb. 1. Die Kragplatte im 
+ D(s) cos (s + 2)0 | Koordinatensystem. 


anschreiben, worin A(s), B(s), C(s) und D(s) gewisse von s abhangige Koeffizienten sein sollen. 
Der Mellinschen Umkehrformel! entsprechend sei jetzt die Biegeflache der Platte in der Form 


Tee ee 
w(r, 0) = ae | r*@(s)ds (4) 
ausgedriickt, wobei der Integrationsweg in gewissen Grenzen willkiilich bleibt und die Wahl 
der reellen GréfSe o erst spater erfolgen soll. 


Nun sind zur Bestimmung der vier Beiwerte im Ausdruck (3) ebensoviele Randbedingungen 
erforderlich. Drei dieser Randbedingungen, die die Einspannung der Platte langs des einen 
Randes #=0 festlegen und die tangentialen Biegemomente lings des anderen Randes =a 
gleich null setzen, lassen sich sofort hinschreiben. Sofern vy die Querdehnungszahl des Platten- 
materials bezeichnet, erhalt man 


G40) ree | 


1 dW(s) bat 
ie ao Pa ' ; (5) 


f @W(s) , 1 W(s) , 1 ie aN, | 
= 10 


) 
or? | Tr Or 7 ane 


1 Siehe z. B. R. Courant und D. Hilbert, Methoden der mathematischen Physik Bd. 1, 5. 87, 2. Auflage. 
Berlin 1931. 
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Indem man die Relation (2) beachtet und dann die vorstehenden Gleichungen nach r aus- 
differentiiert, lassen sie sich folgendermafen umschreiben: 


lean | 


d@(s) 

ED leeward | (6) 
d*?@(s) 2 = 

aa + [vs? — (1 —1)]0(9)] == 0. | 


Zur Aufstellung der noch fehlenden vierten Randbedingung bendtigen wir eine geeignete 
Darstellung fiir die an der Stelle rr) der Geraden ?=«a angreifende Kinzelkraft. Unter 
Benutzung der Mellinschen Umkehrformel 1aBt sich eine willkiirliche Funktion p(r) in folgender 


Weise darstellen 
o+ ct co 


1 
1 DU ees | rds [et pedo. (7) 
o— ot 0 


Der gegebenen Lastverteilung entsprechend verschwinde die Funktion p(r) iiberall auBer 
in r=r, und geniige iiberdies der Bedingung 


[p@)ar Say (8) 
0 
Die Integralformel (7) liefert fiir eine derart gewahlte Funktion die Darstellung 
P o+ ct 
Ves 
PO) = a0 | (F) ds. (9) 


Der Gré8e o kann man dabei einen beliebigen Wert geben, da ja das Integral 


C= PO) dO = ET, a (10) 


18 


fir 0 <r, <oo absolut konvergiert. 

0 & 

Nun driicken sich die am Plattenrand «a wirkenden, in der Richtung nach unten als 
positiv angenommenen Scherkrafte einschlieBlich der Zusatzkrafte aus den Torsionsmomenten 


durch das Integral 


o+ ct 
hy DS 2—v ®W(s) , 2v—1 Ws) , 2(1—r) AW(s) 1 &W(s) 
diss Qect i r or? 0d 1 r or od | r3 Oo r 068 ds (11) 


aus, worin D die Biegesteifigkeit der Platte bezeichnet. Benutzt man wieder die Gleichung (2) 
zur Substitution, so ergibt sich 


o+ct 
D d* @ 9 = 
oe; | [aa tLe +A — 9) +0 +2] Hy Gt as, 
cakes idee ayes 


Setzt man jetzt 
3 . 
a= (2) pO. (13) 


so entspricht der Verlauf der Randscherkrafte immer noch der gegebenen Randbelastung P(r), 
die ja iiberall, auBer in r=ry, Null ist. Nach Einfithrung der Ausdriicke (12) und (9) in die 
Relation (13) ergibt sich nun die Randbedingung 


dB O(s iN r,® 2 
fo +18 40-9 tHe +2 Ol Pale (14) 


die im Zusammenhang mit dem System der Gleichungen (6) ausreicht, um die vier Koeffizienten 
A(s) bis D(s) der Funktion (3) zu bestimmen. Das Ergebnis dieser Rechnung lautet dann auf 
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Grund von (4) 

o+ ct 
Prr, [(-  EO,—F®, fete d 
4niD Tee aniceaTl os 


o—wi 


r 


sofern zur Abkiirzung 
O, = (s + 2) sins &— ssin (s +2) 8, 
O, = scos s } — scos (s +2)%, 
E = (1 —»)s cos sx —[s(1 —v) +4] cos(s +2) a, 
a (deep (sit 2) cin 6a = [31 py |e ein (s <2 2) ae its) 
Ne 41 +9 + (1 —?) cos(2s + 2) a] — 
— (1 — v)? [s(s +. 2) —(s +1)? cos 2a eos (2s +2)a] | 


gesetzt wird. Fiihrt man schlieBlich als neue Variable 


w= —(s-+1)i (17) 


ein und wahlt den Integrationsweg o=—1, so 1aB®t sich der Ausdruck (15) nach einigen Ver- 
einfachungen auf die reelle Form 


/ 


oe} Tr 0 rT 
pee Part Geos (win 7) + Hsin ne oF 18 
aD Nu (1+ u?) (18) 


0 
bringen, wobei 


G = [2 Go) aucosa + (1 —y)u Git ausin a] (Gin uf cos & — u Coj u# sin Y) + } 
+[(1 + )u Gin aucosa + (2 + u® — yu?) Go) cusin a] Gin ud sin 3, | 
H = (1 +) [(Ginaucos « —u €oj ausin a) Gin u sin & — (19) 
— Sin «usin a (Gin u# cos H — u Coj ud sin V)], | 
N=4+(1 +7)? + (1 — 7) (3 +7) Co) 20u + (1 —v)?u? (1 — cos 2a) 
sein soll. 

Man kann die Biegeflache (18) und ahnlich gebaute Liésungen fiir andere Randbedingungen 
sehr wohl auch unmittelbar, aus einem passend gewahlten reellen Ansatz, gewinnen. Hat dabei 
die keilférmige Platte nur ecingespannte oder gelenkig gelagerte Rander, so verschwindet in der 
Lésung, wie sich nachweisen laBt!, das Glied mit sin (uln r/ry) und der direkte Ansatz gestaltet 
sich dann verhaltnismaBig einfach. Bei anderen Randbedingungen sind die Sinus-Funktionen 


mitzuberiicksichtigen, und die Anzahl der Freiwerte im reellen Ansatz vergréfert sich dermaBen, 
daB die Transformationsmethode dann entschieden den Vorzug verdient. 


3. Durchbiegungen und Einspannmomente der Platte mit Einzellast. Wir wollen die Lésung (18) 
dazu benutzen, den Verlauf der Durchbiegungen w, langs des Randes #=«, sowie den Verlauf 
der tangentialen Biegemomente m,, langs des Randes #—0 zu verfolgen. Nimmt man noch, 
wie vielfach itiblich, »y 0 an, so lauten die betreffenden Ausdriicke 


Cc 


2 P rr, (Gin 2xu—usin 24) cos (win 7) en 
a aud Nu (1 +u2) ’ (20) 
0 
co r P r 
Pree (u mn) + Jsin (un =) ; F 
ics aaa ae’ | N( +2?) - ey 
0 
mit den Abkirzungen 
IT—uGinaucosa +(2 + uv?) oj ausina, | (22) 


J= Gin aucos a —u Codj ausina, 
wahrend N sich aus (19) ergibt, indem man dort y gleich Null setzt. 


1 §. Woinowsky-Krieger, The bending of a wedge-shaped plate, Journal of Applied Mechanics, paper 
No 52—A-6 (erscheint demniachst). 
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Die Integrale (20) und (21) wurden fiir die Offnungswinkel «= 2/4 und «=271/2 ausgewertet. 
Dies geschah nach der Simpsonschen Formel und zwar in Schritten von 4u=0,1 und in einem 
Intervall 0 <u <u, mit uy=3 bis uy=6 je nach der Konvergenz. Der Rest wurde dann nach 
entsprechender Vereinfachung des Integranden geschlossen ausgewertet. Aus den Abb. 2 und 3 
die die numerischen Ergebnisse wiedergeben, sieht man, da bereits bei einem Offnungswinkel 
von «=2/2 die Randmomente m, mit r— 0 gegen minus Unendlich gehen. Dasselbe gilt 
offensichtlich fiir die radialen Biegemomente langs des freien Randes bei ebenfalls r — 0, wahrend 
in der Nahe der Last diese Momente stets positiv werden und an der Stelle r= ry eine logarith- 
mische Singularitat! aufweisen miissen. ; 


Abb. 2. Einspannmomente m4, und Durchbiegungen Abb. 3. Einspannmomente me, und Durchbiegungen 
des freien Randes wy beim Offnungswinkel ~ = 45° des freien Randes wy beim Offnungswinkel « = 90° 
(Querdehnungszahl » = 0). (Querdehnungszahl y = 0). 


4, Anwendung der Residuen-Rechnung. Enthalt der Spannungszustand der Platte im Eck- 
punkt r=0 eine Singularitat, so erweist sich die Lésungsform (18) zur Darstellung der Platten- 
beanspruehungen in der nachsten Umgebung dieser Stelle als wenig geeignet; in der Tat, die 
Form des Ausdruckes (21) fiir die Randmomente gestattet nicht einmal, die Frage nach dem 
Vorhandensein einer solchen Singularitat direkt zu beantworten. 


Man kann sich indes die Tatsache zunutze machen, daf das Integral (18), dessen Integrand 
eine gerade Funktion von uw ist, eine bequeme Anwendung des Residuensatzes erméglicht. Dieser 
Vorgang setzt vor allem die Kenntnis der Nullstellen der Funktion N(u) im Nenner des Inte- 
granden voraus. Da die Wurzel der Gleichung N(u)=0 sich nur fiir «= 27 und a= explizite 
angeben lassen, so wird die folgende Untersuchung auf den Fall «=a beschrankt und den 
anderen Fall kénnte man ganz analog behandeln. 


Mit «=z geht die keilférmige Kragplatte in eine Halbebene iiber, deren Rand langs seiner 
einen Halfte (®=—0) eingespannt, langs der anderen Halfte (#7) frei ist und in r=ry, 0=2 
der Voraussetzung gemaB eine Einzelkraft P tragt (Abb. 4). Da in erster Linie der Bereich 
r<7rq interessiert, so schreiben wir die allgemeine Lisung (18) in die Form 


‘ Lg f) 
cy ia Hsin (win ™ 


aD Nu (1 + u?) 


du (23) 
0 
um. Mit «=z und den fritheren Bezeichnungen (19) erhalt man 


G = 2 Coj mu(u Coj ud sin b — Gin ud cos 0) — (1 +7)u Ginazu Gin ud sin 
H= — (1+) Ginzu GinuP sind, (24) 
N=4+(1 +»)? + (1—») (3 +7) Gof 22u. 


1 Siehe z. B. A. Nadai, Die elastischen Platten, S. 204. Berlin 1925. 
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Die Anwendung des Residuensatzes auf das Integral (23) fithrt nun unter Beachtung der 
vorliegenden Form des Integranden! auf den Summenausdruck 


_2Prr 
Reis a 


)) fir die elastische Flache der Platte. Gilt zuniichst r<rp, so sind unter R, die Residuen der 
Funktion G exp(ui In r/r))/Nu(1 + u*) und unter R, die Residuen der Funktion 


— Hexp(uiln r/r))/Nu(1 + u?) 


in samtlichen Polen des Integranden zu verstehen, deren imaginarer Teil positiv ist. Das Re- 
siduum der beiden Funktionen in u=1 ist null und alle anderen Polstellen sind in der Formel 


(xi R, +x2R,) (25) 


aad k ii 3 
ot ey (26) 
mit der reellen Zahl 
ee Ee 5+2y+2 
i = Wr oj Re yaaa (27) 


und m=1,3,5,... enthalten. Nach Ausrechnung der Werte R, und R, erhalt man den im 
Bereich r< ry giiltigen Ausdruck 
ko 


& 4 Per < ee weg Cee ee Ley 
7 DVG—s)(G-49) z > G oh rraruree some 2 


== Is 


pee ae 
toed fics a 4 we 2 een es f 


ko 
—(mutkaje 2 sin dcos ee + 


(28) 


ko 


m+1 
+ (ku —miye * sin sin MH] =e) 


wo! 3 


2 A? + pe? ia 


worin zur Abkirzung 
A= k3 +4k—3km’, | 
See 2D 2 
Lt m3 + 4m + 3k?m, (29) 


T9 


aia, 
if 


gesetzt ist. Die elastische Flache im Bereich r> 7 ergibt sich durch einfache Vertauschung der 
GréBen r und r, in den Formeln (28) und (29). Die Reihe (28) samt ihren partiellen Ableitungen 
_»zch r und ry konvergiert sehr gut gerade bei kleinem Argument r/ry und eignet sich daher vor- 
ziiglich zur Darstellung der Spannungskomponenten in der Nahe der Plattenecke. Man erhalt 
insbesondere fiir die Durchbiegung des freien Randes =z den Summenausdruck 


L co ol fiver * 
oo .—'— ae 5) 
eee teehee 2 2 (*) (30) 
7 #D(1—v)(3+7) 2? + pw? To 
ee ae 
und fiir die tangentialen Momente langs r—0 den Ausdruck 
4 Pry SP sey s kl 
m, = > [(m? + 2m — k?) A — 2 (m + 1)k ys] cos — + 
‘3 anrVa—)3+%), —, | | 2 
See — ‘ (31) 
AEs 2 air \e 
+ [(m? + 2m —B) +2 (m + 1)k Ajsin HE ON (2 | 


‘die beide fiir r<ry gelten. Der Verlauf dieser GréBen fiir y= 0, also fiir k= 1/z Ut Co] 5/3 = 0,34970 
ist in Abb. 4 zur Darstellung gebracht. 
1 Vel. etwa E. T. Whittaker und G. N. Watson, A course of modern analysis, $. 115, 4. Auflage. Cam- 


bridge 1927. 
28 
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Wie der Ausdruck (31) lehrt, wird das Einspannmoment an der Plattenecke von der Ordnung 
(r,/r)/2 mit negativem Vorzeichen unendlich. Daf das absolut genommene Einspannmoment 
mit zunehmendem 1/ry) sehr rasch abnimmt und dafs die einklemmende Wirkung des Randes 
sich praktisch auf einen (relativ zu rp) sehr kleinen Bereich an der Ecke erstreckt, ist wohl auf 
die besondere Stellung der Einzellast — in der Fortsetzung des eingespannten Randes — zuriick- 
zufithren. Man kann schlieBlich unschwer nachweisen, da® auch die Auflagerkrafte mit r/T)—>0 
unbegrenzt anwachsen und zwar in der Ordnung von (7)/r)°””. 

5. Kragplatte bei beliebig verteilter Belastung. Bisher wurde nur eine am Plattenrand 
angebrachte Kinzelkraft in Betracht gezogen, die auch der Lésung (18) zugrunde liegt. Bei 
einem beliebigen Gesetz der Lastverteilung 
lings des Randes ergibt sich die zugehorige 
Biegefliche durch einfache Quadratur des 


0 


Abb, 4. Einspannmomente ms, und Durchbiegungen des freien Randes Abb. 5. Kragplatte mit konzentrischer 
wy beim Offnungswinkel ~ =180° (Querdehnungszahl » = 0). Linienlast. 


Ausdruckes (18) in der Verdnderlichen ry. Da die Gré8e ry nur in Potenzen (mit komplexem 
Exponent) auftritt, ist diese Quadratur bei jeder praktisch méglichen Lastverteilung in ge- 
schlossener Form ausfiihrbar. 

Etwas komplizierter wird die Frage im Fall einer im Inneren des Plattengebietes, etwa in 
r=T, 0=09,)( <a) angreifenden Punktlast. Macht man dann fiir w(r, #) in den beiden Gebieten 
0<d<d, und 3) <3 <a zwei verschiedene Ansatze und zwar jedesmal von der Form (4) in Ver- 
bindung mit (3), so stehen zur Bestimmung von insgesamt acht unbekannten Koeffizienten, 
etwa A,(s), B(,s),..... ALS) 5) Big(S)oneterotare folgende Bedingungen zur Verfiigung: zwei geo- 
metrische Bedingungen langs des eingespannten Randes ?=0; zwei Bedingungen fiir das Ver- 
schwinden der Schnittkrafte langs des freien Randes #=o; drei geometrische Kontinuitats- 
bedingungen langs des Halbmessers = @%, und zwar fiir w, dw/d00 und d?w/d0*; schlieBlich 
eine Diskontinuitatsbedingung, die die Differenz der Scherkrafte — D(1/r) 0(4w)/00 zu beiden 
Seiten von #= %) und die Lastfunktion (9) in der Form der Relation (13) gleichsetzt. Die so 
erhaltene EinfluBfunktion w(r, 0, ro, 09) wiirde es dann erlauben, die Wirkung einer in beliebiger 
Weise verteilten Belastung zu erfassen. 


6. Kragplatte mit konzentrischer Linienbelastung. Bei gewissen einfacheren Lastanordnungen 
erweist es sich als zweckmafig, die Konstruktion einer EinfluBfunktion zu umgehen und die Lésung 
unmittelbar auf der Mellinschen Umkehrformel aufzubauen. Die Kragplatte trage beispiels- 
weise eine gleichférmig iiber den Kreisbogen vom Halbmesser ry verteilte Linienlast von der 
Intensitat q (Abb. 5). Die elastische Flache der Platte hat in diesem Falle der Differential- 
gleichung 


o+ cot —(s-+ 4) 
—_ a a 
NN lore, | T, ds. (32) 
zu geniigen, deren rechter Teil aus (9) hervorgeht und zwar durch Multiplikation mit = To ; 
lg 


Man iiberzeugt sich leicht, daB dieser rechte Teil die durch D dividierte vorgeschriebene Last- 
verteilung als Funktion von r darstellt. Die Lésung der Gleichung (32) kann nun in der Form 


w= Wy + Wy (33) 
angesetzt werden. Hierin ist 
, a+ colt 
RL apy pe 1 
VOT SD | & s? (s +2)? (34) 
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eine Partikularlésung von (32), wahrend w, als Lésung der homogenen Gleichung AAw = 0 
in der Form (4) angenommen werden darf. Fiir die Gesamtlisung w,-+w) miissen die Durch- 
biegung und die Neigung der elastischen Flache langs )=0 wee inttden desgleichen die 
tangentialen Biegemomente und die Auflagerkrafte lings }— «; diese vier Better t scleeiien 
sich in Ubereinstimmung mit (6) und (12) folgendermafen hin: 


ee 3 
O19), 0 ae ies cee 
d@(s) ted 
baw gui 
ee ne a 1 aatt teeny (35) 
| di | ‘ eG) am +r Ds (s+ 2)? —=0, 
O(s 2 
Spe te ta—y e+ ye +O! mo, 


10 
Sind die Koeffizienten A(s) bis D(s) der Funktion O(s) aus den Gleichungen (35) bestimmt 
so bleibt es nur noch iibrig, die Gesamtlésung (33) durch die frithere Substitution u= —(s ean 
in die reelle Form itiberzufiihren. Das Endergebnis sei nur fiir den Sonderfall y—0 angegeben. 
Die beiden Bestandteile von (33) sind dann 


oO 7 2 
ait cos (u In 4 a A (14m cy ; (r <1) | 
0 2D Que ON) os ( se (36) 
4 4D (14+m*), eat 
grr “ K cos (win + Lsin (un) 
__ 4To r 
ae hed | Nu @ + u?)? = du (37) 


mit 

K = (u? sin 2a — 3u Gin 2au — 4 Coj au sin a) (u Coj ud? sin  — Gin ud cos 7) + 
+ (3u? Coj 2au + u? cos 2a — 4 Co) aucosa + 4u Gin au sin a) Cin ud sin & — 
— Nu Goj ud cos 3. 

L = 2 (u ©oj ausin « — Gin au cos a) (u of ud sin 9 — Gin u# cos 9) + 
+ 2[(1 — u?) Sin ausin a + 2u(Coj aucosa —1)] Ginud sin #, 

N=5+3 @oj 2au + u? (1 — cos 2a). 

Das Integral in (37) konvergiert etwas besser als das analoge Integral in (18) und 1aBt sich 
in ahnlicher Weise numerisch auswerten. Man erhalt beispielsweise an der Stelle r=ry, 0=a 
und fiir die Offnungswinkel «=—2/4, 2/2 und a die Durchbiegungen k Qr,?/D mit k=0,0305; 
0,0951 und 0,123, sofern Q= qr,« die gesamte, iiber die Bogenlange ry verteilte Last bezeichnet. 
Verteilt sich schlieBlich die Last gleichmaBig auf eine durch die beiden Radien 0=0, 0D=« 
und durch zwei Kreisbégen mit dem Mittelpunkt in r=0 begrenzte Flache eines Sektors, so 
ergibt sich die zugehérige elastische Flache aus der Lésung (36) bis (37) wiederum durch einfache 
Quadraturen. 

7. SchluBbemerkung. Wie die durchgefithrte Untersuchung zeigt, lassen sich durch die An- 
wendung der Mellinschen Umkehrformel bequem Lésungen konstruieren, die beliebige Last- 
bedingungen und auch beliebige Randbedingungen langs der beiden Grenzradien der Platte zu 
erfiillen vermégen. In einer etwas abgednderten Form eignen sich diese Lésungen auch dazu, 
nihere Auskunft iiber die etwa vorhandene Spannungskonzentration in der Plattenecke selbst 
zu geben. Die angewandte Methode ist sicherlich nicht hinreichend, um zu gleicher Zeit und in 
strenger Weise Randbedingungen an etwaigen weiteren Begrenzungslinien des Plattengebietes zu 
befriedigen; es ist aber wohl denkbar, ein derart erweitertes Problem durch Kombination einiger 
elastischen Flachen von dem hier behandelten Typus mit einer beliebigen Annaherung zu lésen. 


(Eingegangen am 29, Mai 1952.) 
Anschrift des Verfassers: Prof. Dr. S. Woinowsky-Krieger, Quebec (Canada). 200, av. Marguerite 
Bourgeoys, Apt. 8. 
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Ein Verfahren zur Beurteilung des Giitegrades von Mischungen. 
Von K. Stange. 


1. Aufgabenstellung. Grundbegriffe und Voraussetzungen. Es werden zwei ,,kérnige“ Stoffe 
(P) und (Q) miteinander ,,gemischt‘+, und zwar so, da das Gesamtvolumen V der Mischung 
die absoluten Anteile Vp bzw. Vg oder die relativen Anteile 


Vp ioe (1.1 
Ww = iP bzw. fn = Q ) 
beider Stoffe enthalt. Die Werte P und Q sind bekannt. Wegen Vp + Vo = V ist 
Oi Ae (1.2) 


Durch eine Stichprobe S, die man der Mischung entnimmt, und die praktisch erheblich 
kleiner als V sein soll, will man beurteilen, ob an der Entnahmestelle gut durchmischt ist. 
Es ist einleuchtend, daB die GréBe einer solchen Stichprobe nicht ganz beliebig sein darf. Das 
zeigt die Betrachtung der beiden Grenzfalle: 

a) Die Stichprobe S entspricht dem Gesamtvolumen V. Dann ist kein Urteil itber dieMischung 
moglich, denn man findet bei der Analyse nur das triviale Ergebnis, dafs die Anteile von (P) 
und (Q) genau den Werten P und Q entsprechen. 

b) Die Stichprobe hat nur die GréBe eines einzelnen Kornes oder Elements der gemischten 
Stoffe. Dieses Teilchen gehért dann entweder zu (P) oder zu (Q); und es ist auch dann kein 
Urteil iiber die Mischung méglich. 

Die Frage nach der notwendigen Gré8e einer Stichprobe erscheint danach sinnvoll. Umgekehrt 
kénnen wir unter der Voraussetzung guter Mischung danach fragen, wie groB das kleinste Vo- 
lumen v= v ist, welches innerhalb vorgeschriebener Genauigkeitsgrenzen in bezug auf die Be- 
standteile (P) und (Q) als homogen angesehen werden darf. Es soll versucht werden, diese Fragen 
im folgenden zu beantworten!. 

Wir legen der Untersuchung ein einfaches mathematisches Modell zu Grunde, das die folgenden 
sechs Voraussetzungen erfiillen soll: 

I. Die Mischung besteht nur aus zwei Stoffen (P) und (Q). 

Il. Die ,, Korner“ seien kugelférmig gestaltet, so daB eine méglichst liickenlose dichte Packung 
im Raume entsteht. 

III. Die Einzelelemente (Kérner; Kugeln) seien fiir jeden einzelnen Bestandteil (P) bzw. (Q) 
der Mischung gleich gro. Die verschiedenen Bestandteile werden durch unterschiedliche Farben 
gekennzeichnet; sie sollen entweder schwarz = (P) oder weiB = (Q) sein. 

Das Volumen eines einzelnen Teilchens sei ep bzw. ¢9; ihre Anzahl im Gesamtvolumen sei 
Np baw. No. Dann ist 


Np ep + No eg = V. (153) 
Bezeichnen wir die Gesamtzahl der Teilchen im Volumen V mit N, so ist 
Np No= N: (1.4) 


Die bereits in (1.1) erklarten Quotienten Vp/V = Pund Vo/V = Q bezeichnen wir als Volumen- 


1 Durch Herrn Dr. Weidenhammer, Clausthal, wurde ich darauf hingewiesen, daB die Frage, wie man 
Mischungen von kérnigen Stoffen beurteilt, schon mehrfach behandelt worden ist, u. a. von B. Baule 
und A. Benedetti- Pichler [Z. analyt. Chem. 74 (1928), S. 442], U. Graf und H. J. Henning [Z. Erzberghb. 
Metallhiittenwesen 5 (1952), S. 127], A. B. Manning [Fuel XX VIII—3 (1948), S. 49] und J. Visman 
[De monsterneming van heterogene binomiale korrelmengsels, in het bijzonder steenkool, 94S. Gro- 
ningen 1947]. In diesen Arbeiten steht jedoch eine andere Fragestellung im Vordergrund: Eine Mischung 
besteht aus zwei kérnigen Stoffen(P) und (Q), die beide einen interessierenden dritten Stoff (X)im Ge- 
wichtsverhaltnis Xp [%] und Xg [%] enthalten. Gefragt wird, wie groB eine der ,,Gesamtmischung“ ent- 
nommene Stichprobe S mindestens sein mu, damit man aus der Analyse der Probe ein mit vorgeschriebe- 
ner Wahrscheinlichkeit zutreffendes Urteil dariiber abgeben kann, wieviel Gewichtsprozent an Stoff (X) 
die ,,Gesamtmischung*‘ enthalt. — Im folgenden wird eine ganz andere Frage erértert. Thre Bearbeitung 
wurde vom Leiter des Instituts fiir Chemische Technik der Technischen Hochschule Karlsruhe, Herrn 
Prof. Dr. F. A. Henglein, angeregt. 
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haufigkeit, die Quotienten 


Np N, 
7 Pp und y=4 (1.5) 
als Teilchenhaufigkeit der Mischungsanteile (P) und (Q). Selbstverstandlich ist wegen (1.4) 


auch 


p+, (1.6) 


Der Zusammenhang zwischen den Teilchenhaufigkeiten (p; gq) und den Volumenhiufigkeiten 


(P;Q) ist durch die folgenden Bezichungen gegeben, die man leicht aus den vorausgehenden 
Grundgleichungen herleitet : 


as Plep Qleg 
= her,” T= Pep -+ Oleg’ 7) 
oder umgekehrt 
ae PEp a q£o 
peeesGes 98” Q= pESeGes a5" (1.8) 


IV. Zunachst sollen die Teilchen (P) und (Q) gleich groB sein, also ep = ég. GemaB (1.8) 
und (1.6) wird dann 


pL und 708 (1.9) 


so dafs wir Volumenhaufigkeiten und Teilchenhiaufigkeiten vorerst nicht zu unterscheiden 
brauchen, was die Erérterungen ganz wesentlich vereinfacht. Spiter werden wir uns von dieser 
Einschrankung freimachen. 


Selbstverstandlich stellen I. bis IV. schwerwiegende einschrankende Bedingungen dar, 
denen in Wirklichkeit kaum eine Mischung entspricht, und es ist nicht leicht zu iibersehen, ob 
eine soweit vereinfachte Modellvorstellung iiberhaupt fiir die Erfassung wirklicher Verhdltnisse 
niitzlich ist. Schieben wir diesen Einwand vorlaufig einmal zuriick und stellen wir uns auf den 
auch sonst zweckmaBigen Standpunkt, dafi man mit einfachsten Fallen beginnen und darauf 
weiter aufbauen mu, dann kénnen wir uns die Entstehung der Mischung nach Abb. 1 folgender- 
mafen vorstellen: Je ein Behalter ist mit schwarzen (P) und wei®en (Q) Kugeln gefiillt, ein 
dritter, in dem die Mischung entstehen soll, ist leer. Jetzt lassen wir einen Zufallsvorgang 
(Wirfelspiel, Roulette, oder dergl.) ablaufen, bei dem zwei Ereignisse (P) und (Q), die wir den 
Mischungsanteilen zuordnen, mit den Wahrscheinlichkeiten p und q sich einstellen. Jedesmal, 
wenn das Ereignis (P) eintritt, legen wir ein Element e¢p in das dritte Gefé%. Entsprechend 
verfahren wir mit ¢9, wenn der ZufallsprozeB das Ergebnis (Q) liefert. 

Auf diese Weise entsteht in dem dritten Gefa8 eine Zufallsmischung der Stoffe (P) und (Q) 
mit den Teilchenhdufigkeiten p und g. Die mathematischen Gesetzmafigkeiten, denen der 
ZufallsprozeB bei wachsender Anzahl] der beteiligten Elemente folgt, sind bekannt. Wir brauchen 
nur die Ergebnisse der mathematischen Statistik nutzbar zu machen. Zu dem Zweck fassen 
wir m zusammen liegende Teilchen zu einem Raumteil S, der Stichprobe, zusammen. Wir wollen 
die Raumteile S so wahlen, da®B sie alle gleiche Gestalt haben und selbst wieder eine Aufteilung 
des Gesamtraumes ergeben. (Das ist eine an sich unwesentliche Einschrankung, die wir nur 
zur besseren Veranschaulichung einfiihren). Wir setzen weiter voraus, daB 

V. die Anlagerung benachbarter Teilchen unabhangig voneinander vor sich geht, was gleiches 
oder wenigstens nahezu gleich 's spezifisches Gewicht der Bestandteile erfordert. Diese Forderung 
schlieBt bei wirklichen Mischungen auch Faille aus, in denen etwa Teilchen aus Kupfer und 
magnetisiertem Stahl miteinander gemischt werden sollen. Dann ist die Wahrscheinlichkeit 
AW(n,,n,;n), daB ein aus n Teilchen bestehender Raumteil S(n) genau n, schwarze und ny, 
weiBe Kugeln enthalt, bekanntlich 

n! 


AW (n,, Ny3 1) = i el pda mit n,-+ny=n wd, O0srngs nm, J (110) 
2 hs ig 


Betrachten wir nun die Gesamtheit der Raumteile S unter der weiteren Voraussetzung, 


VI. daB S <V oder n <N ist, daB also die Zahl der Stichproben, die aus dem Gesamt- 
volumen V entnommen werden kénnen, sehr grof ist, oder daB die Teilchenzahl n der Stichprobe 
klein im Vergleich zur Gesamtzahl N bleibt, so verteilen sich die Stichproben in der Mischung 
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nach der Binomialverteilung 


n 


n! iE Bs 5 
(p+q)"= > ai ce ee mit n,-+ny =n. (1.11) 
n, = 0 A y 

Damit ist die Aufgabe grundsatzlich schon gelést, aber die Lésung ist insbesondere fiir groBe n 
noch nicht praktisch verwendbar. Um sie brauchbar zu machen, bendtigen wir ein paar be- 


kannte Tatsachen der mathematischen Statistik. 


2. Die Wahrscheinlichkeitsverteilung lings der Stichprobengeraden. Durch (1.10) bzw. (1.11) 
ist in der (n,, my)-Ebene eine Wahrscheinlichkeitsverteilung festgelegt. Es handelt sich um eine 
entartete ebene Verteilung, bei der nur gewissen Gitterpunkten (n,, n,) der Geraden n, + ny=n 
eine von Null verschiedene Wahrscheinlichkeit zu- 
geordnet ist, und zwar den Gitterpunkten des ersten 
Quadranten (Abb. 2). Tragen wir iiber jedem Punkt 
(n,, my) der Geraden n, +-ny—=n die ihm ent- 


sprechende Wahrscheinlichkeit AW(n,, ny; n) auf 


iS iY is Z A 
TAD 
Oma, 


Xe) 
RY 
XK 
@\ 
\ 
IN 


CX WW 
NOOO 
NO \S\ 
\G OO! 
NOOO 
COO\ 
CWEX 


eV 
WX a 
© 
@ 
. 
KX 


IN 


WYUmUYS a--0 for 
G GeeZgZ 0 2 % o 8 = Te 
bX 1) GZ 0 b2 OY 06 08 40 
Abb. 1. Zur Entstehung einer Zufallsmischung. Abb. 2. Die ,,ebene‘* Wahrscheinlichkeitsverteilung 4 W (ny, nys n) 


Eine Stichprobe S(n) umfaBt n Teilchen. fir p = 0,4; q= 0,6 und n= 10. 

und legen die entstehende Verteilung in die (n,, ny)-Ebene um, so finden wir die in Abb. 2 ge- 
gebene Darstellung. ZweckmaBig fiithren wir hier noch anstelle der absoluten Haufigkeiten (n,, ny) 
fiir das Vorkommen schwarzer und weiBer Kugeln die relativen auf die Stiickzahl n der Stich- 
probe bezogenen Werte 


mn 
nN, und oh 
mn 


sy) mit *-+-y= 1 (2.1) 
ein. 

Die in Abb. 2 durch einen Kamm dargestellte diskrete binomiale Verteilung darf nun fir 
geniigend groBe n in guter Naherung durch eine Normalverteilung mit dem Mittelwert m(p; q) 
und der auf den Mittelwert bezogenen Streumatrix 


p(1—p) Pq 
Oxx Oxy eas a —— = 
C= = (2.2) 
a a Axoap q(1—q) 
ae XY, n nm 


ersetzt werden!. DaB es sich hier um eine entartete ebene Verteilung handelt, bei der die 


1 R. von Mises, Wahrscheinlichkeitsrechnung, S. 418. New York 1945. — H. Cramér, Mathematical 
Methods of Statistics, S. 318. Princeton 1946. 
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Gesamtwahrscheinlichkeit in Wirklichkeit langs einer Geraden verteilt ist, kommt dadurch zum 
Ausdruck, daB die Determinante der Streumatrix verschwindet, 


lela —4 
oS | Ey on Panes 04 (2.3) 
P q 
was man wegen p -+-q=1 sofort sieht. 

Wesentlich ist, daB die einzelnen Glieder der Streumatrix mit wachsender Teilchenzahl n 
der Stichprobe gegen 0 streben. Die Verteilung zieht sich also mit wachsendem n mehr und 
mehr auf eine kleine Umgebung des Punktes m(p; q) zusammen. 

Fiir die weitere Untersuchung haben wir nun zwei Méglichkeiten. Entweder betrachten wir 
die Verteilung als eine zweidimensionale iiber der (x, y)-Ebene liegende ausgeartete Verteilung. 


Dann bleibt der Vorteil der Symmetrie in den Veranderlichen x und y erhalten. Oder wir pro- 
jizieren die Verteilung 


AW(x, y) = AW(x, 1 — x) = Funkt. (x) = AW(1 — y, y) = Funkt.(y) 


auf die x- oder auf die y-Achse. Dann erhalten wir dort eine eindimensionale Verteilung, haben 
aber entweder x oder y einseitig zum Nachteil der anderen Veranderlichen bevorzugt. Im Grunde 
genommen ist ein Weg so bequem wie der andere. Von Bedeutung werden diese Unterschiede 
erst bei drei oder mehr Komponenten in 


der Mischung. (nr 2 of 
8 plr-py” g 
iY 
WA \x-7| Ann} 7 
6 
i) 
4 
3 
2 4 a 
of =- 
asl oh pot abn wae 
— 
100 98 6 4 Se 10 Fy? 
eee eee eS | a Lee ————l a Se ee 
0 a, 4 6 8 0 
Abb. 3. Zur Berechnung der Wahrscheinlichkeit Abb. 4. Zusammenhang zwischen Wahrscheinlichkeit W, Irrtumswahr- 
W {|x— p| < 4p; n}. scheinlichkeit 1— W = J und Teilchenzahl n der Stichprobe S(n). 


Wir berechnen nun die Wahrscheinlichkeit W, daB bei gegebenen Werten fiir die Teilchen- 
haufigkeiten (p,q) und gegebener Teilchenzahl n der Stichprobe die beobachteten relativen 
Haufigkeiten (x, y) in den in Abb. 3 geschrafften Bereich 


|s—p|<Ap bzw. |y—q| <4q= Ap (2.4) 


hineinfallen. Diese Wahrscheinlichkeit hangt von p, Ap und n ab und ist mit praktisch aus- 
reichender Genauigkeit durch die Flache unter einer normierten Normalverteilung gegeben, also 


durch 


| 1 
Wi\\x—p| < Ap; n Tee 


22 


[- 2 at, (2.5) 


2 


wobei die Integrationsgrenzen -+ @ aus der Beziehung 
@=o(p. Ap, n) = = yn (2.6) 


berechnet werden miissen. Abb. 4 stellt den durch (2.5) und (2.6) vermittelten Zusammenhang 
zwischen den Veranderlichen p, Ap, n und W dar. Zu drei gegebenen Gréfen laBt sich die vierte 
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einfach berechnen. Zu p=0,4, Ap= 0,04 und W=0,95 findet man z. B. mit 9? = 3,84 die not- 
endi i nzahl n= n ~ 580. ' 
Berne mit Abb.5 folgendermafen. Setzen wir einmal voraus, daf das Volumen V mit 
den Stoffen (P) und (Q) im Verhaltnis p:q= 0,40: 0,60 gut durchmischt ist und entnehmen wir 
einer sehr groB gedachten Grundgesamtheit V nicht nur eine, sondern zahlreiche Stichproben 
S,(n), S,(n), S,(n), .. . von gleicher Teilchenzahl n= 580, so werden nach der eben durchgefithrten 
Rechnung fast alle diese Proben. (namlich 95%) einen relativen Anteil x im Bereich I oder 
0,36 <x <0,44 aufweisen. Nur 5% der Proben besitzen davon abweichende Werte x;; mit 
|x17 — p| > Ap = 0,04 im Bereich II. Entsprechendes gilt sinngemaf fiir andere Zahlenwerte. 


3. Der statistische Schluf. Jetzt sei um- 
gekehrt eine Mischung V vorgelegt, aus der 
wir eine einzige Stichprobe S,; mit n= 580 
Teilchen entnehmen. Diese Stichprobe mag 
bei der Analyse fiir (P) und (Q) die Teilchen- 
haufigkeiten (x,, y,) ergeben. Wir erwarten 
(x;, ¥;) im Bereich I. Liegt der gefundene 
Stichprobenpunkt (x,, y,) jedoch entgegen 
unserer Erwartung in dem nicht geschrafften 
Bereich II, so stehen grundsatzlich drei Még- 
lichkeiten offen, aus diesem Tatbestand eine 
Folgerung zu ziehen: 


Abb. 5. Der Schwankungsbereich I fiir die Teilchenhaufigkeiten x 
der Stichprobe bei gegebener Irrtumswahrscheinlichkeit I. 


(A) 
Die Voraussetzungen itber die 
Mischungsanteile treffen nicht 
zu, d.h. das Mischungsver- 
haltnis nach Teilchen ist nicht 
p/q, sondern irrtiimlich ein 
anderes. 


(B) 
Die Voraussetzungen sind 
richtig, d.h. das Mischungs- 
verhaltnis im Gesamtvolumen 
ist p/q. Es ist gut gemischt. 
Aber trotz guter Mischung ist 
eine stark abweichende Stich- 


(C) 
Die Voraussetzungen sind 
richtig, d.h. das Mischungs- 
verhaltnis im Gesamtvolumen 
ist p/q. Es ist schlecht ge- 


mischt. 


probe herausgegriffen worden, 
d.h. es ist ein sehr unwahr- 
scheinliches Ereignis ein- 
getreten. 


Die Méglichkeit (A) schlieBen wir selbstverstandlich aus, obwohl sie in der Praxis vorkommen 
kann, da kein Mensch fehlerfrei arbeitet. Wir nehmen also an, daB® bei der Herstellung der 
Mischung kein Versehen vorgekommen ist. 

Die Méglichkeit (B) schlieBen wir bei statistischen SchluBfolgerungen verabredungsgema4B 
aus. Wir stellen uns also auf den Standpunkt, daf ein Ereignis, dessen Wahrscheinlichkeit sehr 
klein ist, bei einmaliger Realisierung seiner Vorbedingungen, also bei einmaliger Beobachtung, 
nicht eintritt!. 

Dann bleibt also, wenn der Stichprobenpunkt (x, y) im Bereich II liegt, nur die unter (C) 
getroffene Entscheidung ,,es ist schlecht gemischt“ iibrig. 

Die statistische Hypothese H, die wir im vorliegenden Falle mit Hilfe der Stichprobe S 
testen wollen, heifBt: ,,Es ist an der Stelle, an der S entnommen wird, gut durchmischt*‘. 

Liegt der Stichprobenpunkt (x,, y;) in dem 2-Bereich I, so urteilen wir: Die Hypothese H 
ist mit dem Befund des Versuchs vertraglich. H wird angenommen. 

Liegt jedoch (x, y,) auBerhalb des geschrafften Bereichs I im Bereich II, so lehnen wir die 
Hypothese H als unvertraglich mit dem versuchsmabigen Ergebnis ab. Das Urteil lautet dann: 
»,Hs ist an der Stelle S nicht gut gemischt*. 

Da auch bei guter Mischung immerhin noch 5% aller Stichproben im Bereich IT liegen, so 
werden wir bei der eben erlauterten Schlu®weise auf lange Sicht in 5% aller Falle ein Fehlurteil 
abgeben, némlich immer dann, wenn trotz ,,guter Mischung*, also bei zutreffender Hypothese H, 
ein Wertepaar (%,, y,) im Bereich II erscheint, der nur die geringe Wahrscheinlichkeit 0,05 fiir 
sich hat. Die Wahrscheinlichkeit J=1—W+=I—0,95=0,05 heif®t in diesem Zusammenhang 
auch die Irrtumswahrscheinlichkeit. Es ist also wesentlich, dafs man bei allen statistischen 


1 Man vel. dazu die Erérterungen bei B. L. van der Waerden, Studium Generale JOS tes (oy, 
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Schliissen eine solche Irrtumswahrscheinlichkeit zulassen muB. Erscheint J=0,05 als ein zu 
groBer Wert, so steht nichts im Wege, eine kleinere Irrtumswahrscheinlichkeit, etwa J—0,01 
zu wahlen und alle Urteile auf diesen Wert aufzubauen. Die Wahl einer Irrtumswahrschein- 
lichkeit ist kein statistisches Problem. Man ist frei in der Entscheidung, welche Wahrschein- 
lichkeit man als so klein betrachten will, da®8 man mit dem Eintreten des Ereignisses bei ein- 
maligem Versuch nicht rechnet. Es handelt sich um eine Willensentscheidung, die von Fall 
zu Fall verschieden ausfallen wird, und die man in erster Linie nach den Folgen ausrichtet, die 
eine irrtiimliche Entscheidung nach sich ziehen kann}. 

Sind wir damit zufrieden, wenn die beobachtete Teilchenhiufigkeit x der Stichprobe innerhalb 
eines gegebenen Schwankungsbereichs 4p um den Sollwert p der Grundgesamtheit liegt, lassen 
wir also eine Schwankung der x-Werte im Bereich |x—p|<Ap zu, und schreiben wir eine be- 
stimmte Irrtumswahrscheinlichkeit J vor, so entnimmt man aus Abb. 4 die Teilchenzahl n der 
Stichprobe, die zur Beurteilung der Mischung an der Entnahmestelle erforderlich ist. Fir 
J=0,05 und eine zuliassige Schwankung vom Betrage Ap—0,04 um den Sollwert p=0,40 
muf die Stichprobe mindestens aus 7 = 580 Teilchen bestehen (wie wir vorhin schon ausgerechnet 
haben). Umfassen die Stichproben mehr Teilchen (n>7),) so sind diel Schwankungen der 
Haufigkeit x geringer, fiir n<n sind sie gréBer als Ap. 


Hiutigheit/oni : Hitigkeitf oni 


Op : om 
ép 


Abb. 6. Die Haufigkeitskurven fiir die TeilchengréBen. 


4, Losung fiir ungleiche TeilchengréBe ep == &g. Alle bisherigen Betrachtungen galten unter 
der Voraussetzung IV, d.h. wir haben fir die ,,Elementarvolumen“ ep bzw. eg der Hinzel- 
teilchen gleiche Werte ep=¢9 angenommen. Wir wollen uns im folgenden von dieser Be- 
schrankung freimachen. In Wirklichkeit wird nicht nur ep + €g sein, sondern die TeilchengréBe 
wird von Korn zu Korn schwanken. Eine vor der Mischung vorgenommene feinere Analyse der 
TeilchengréRe ¢ und ihrer Haufigkeiten fp und fg mag fiir (P) und (Q) zwei Verteilungskurven 
fe(e) und fo(e) nach Abb. 6 geben. Wir berechnen fiir beide Haufigkeitskurven die Mittelwerte ep 
und ¢g durch 


co oO 


Ep = it é fp(e) de und, -€9 == ‘h €fo(e) de (4.1) 
0 0 
und die auf diese Mittelwerte bezogenen Streuungen op und o@ durch 
ob= f (c—er)*fole)de und of = | (¢ —e0)*fole)de- (4.2) 
0 0 


Wir wollen voraussetzen, daB die StreumaBe op und og ,,klein‘« gegen die entsprechenden 
Mittelwerte ep und ég sind, 

Op = Ep und 0Q << EQ » (4.3) 
und daB ep und eg wenigstens von gleicher Gré®enordnung sind (daf man also, anschaulich 
gesprochen, nicht Meh] und Pflastersteine zusammen mischt.) In diesem Falle wird man anstelle 
der von Korn zu Korn veranderlichen Teilchengré8e die Mittelwerte ep und eg zugrundelegen 
dinfen. Diese Naherung ist in einem bemerkenswert weiten Bereich brauchbar, wie eine weitere 


Untersuchung gezeigt hat. 
Wir haben nun die Ergebnisse, die wir fiir die Teilchenhaufigkeiten (x; y) gefunden haben, 


auf die Volumenhaufigkeiten 


xX. cae os und Ye — = (4.4) 


v v 


1 Hl. Miinzner, Mitteilungsbl. math. Statistik 1 (1949), S. 5, 69 und 185. 
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zu iibertragen. Dabei ist v, das von den n, Teilchen (P) eingenommene Volumen. Entsprechend 
ist vy erklart, und vist das Gesamtvolumen der Stichprobe. 
Es ist ahnlich wie in (1.8) und (1.7) 


wear ye 4.5 
sears Boo i Bed Gs 
und a “i 
es Sp ale by se ey 4.6 
Se GPE EP iy epae NI eD a (2-6) 


Abb. 7 zeigt den Zusammenhang zwischen den Volumenhdufigkeiten (P, Q) und den Teilchen- 
haufigkeiten (p,q) im Gesamtvolumen fiir verschiedene Verhiltnisse eép/eg der Elementar- 
volumen. Einer Mischung mit den vorgeschriebenen Volumenhdufigkeiten (P; Q) und zulassigen 
Toleranzbereichen fiir die Stichproben (AP; 4Q) mit AP + AQ=0 entspricht ein Wertepaar 
(p;q) der Teilchenhaufigkeiten mit zulassigen Toleranzbereichen (4p; 4q). Zur Berechnung 
dieser Bereiche fiihren wir zweckmaBig die HilfsgréBe ¢, das mittlere Elementarvolumen in der 


Mischung, durch 


Vag yursepe LNG 
j= . ER eee apa Ta 4.7 
0 ‘ = €=pép + qég bzw . ne (4.7) 
f H | | ein. Dann findet man aus (1.8) oder epp= éP 
OZ} 08 durch eine kurze Rechnung leicht 
2 
= P 
Ap ipa A 
a4} 06 es 
bgfeg-1 und qr, (4.8) 
[ 2 mit Ap +Aq= 
G6} 04 7h f 
i Fiir die relativen Anderungen folgt daraus 
10 Ap 26 AF Aq__& AQ 
G8 02 pe ips und q c8 Q . (4.9) 
i 50 | Zu W baw. 1—W=J und den eben ge- 
pale ’ fundenen Werten p und Ap bestimmt man 
i oak anh ee ee irs lee ae Ae g nach den Ergebnissen des Abschnitts2 die 
i a a OM Ue Z notwendige Teilchenzahl n=n. Das dazu- 


Abb. 7. Der Zusammenhang zwischen den Volumenhaufigkeiten gehorige Volumen v=v der Stichprobe ist 
(P, Q) und den Teilchenhaufigkeiten (p, g) in Abhangigkeit e 
vom Volumenverhiltnis €plég: ann 


v= (np)ep +(ng)eg=ne. (4.10) 
Damit ist die Gréfe des innerhalb gegebener Toleranzbereiche (AP, AQ) als ,,homogen“ 
anzusehenden Volumens v bestimmt, wenn die Volumenhaufigkeiten (P, Q) gegeben sind. 


5. Proben S(n) gleicher Teilchenzahl und Proben S(v) gleichen Volumens, Mit einem Einwand 
miissen wir uns noch auseinandersetzen. Die in den vorausgehenden Abschnitten entwickelte 
Mischungstheorie verlangte die Entnahme von Proben S(n) gleicher Teilchenzahl n=n. Das ist 
aber nicht so einfach zu verwirklichen. Wir kénnen leicht Proben gleichen Volumens, etwa 1 cm’, 
1dm°,.... entnehmen. Wie sollen wir aber Proben von n—5600 Teilchen oder dergleichen 
erfassen ? Wenn es auch keine uniiberwindlichen Schwierigkeiten macht, Verfahren zur Ent- 
nahme solcher Stichproben zu ersinnen (etwa fiir Laborbetrieb), so wird man praktisch doch nur 
mit Proben S(v) gleichen Volumens v=% arbeiten kénnen. Wenn jedoch die Elementarvolumen 
ép und €g sich voneinander unterscheiden, so sind Proben S(n) fester Teilchenzahl n—n nicht 
zugleich Proben gleichen Volumens, sondern mit schwankender Zusammensetzung der Probe 
mufs bei fester Tcilchenzahl n= 1n das Volumen schwanken, wie es Abb. 8 anschaulich zum 
Ausdruck bringt. Aus Griinden der einfachen und raschen Handhabung wird man nun stets 
Proben S(v) gleichen Volumens v=% wahlen. Dann mu man Schwankungen in der Teilchenzahl 
in Kauf nehmen, und die Frage ist, ob dadurch die Ergebnisse unserer Theorie wesentlich be- 
einfluBt werden, d. h. ob die nur fiir Proben S(n) hergeleiteten GesetzmaBigkeiten, insbesondere 


die Toleranzbereiche, auch fiir Proben S(v) gelten, Wir wollen im folgenden zeigen, daB das inner- 
halb gewisser Grenzen der Fall ist, 
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Zu dem Zweck berechnen wir in Abb. 9 bei gegebenen Werten (p, q) der Teilchenhaufigkeiten, 
zugelassenen Toleranzen (Ap, Ag) und vorgeschriebener Irrtumswahrscheinlichkeit J—1—W 
die Grenzkurven K’ und K’’, welche die Bereiche fiir die Teilchenzahlen n, und ny bei gewahltem n 
abgrenzen. 


Gema8 (2.4) und (2.6) sind diese Grenzlinien in Parameterdarstellung durch 


mn —=np+oe(J)/pqyn und ny=ng+ o(J)/pq yn (5.1) 
gegeben, wobei n der Parameter ist. 
Ist nun n die zu den festen Werten Thy 
(p, Ap, J) nach Abschnitt 2 gehérende 7 7 
notwendige Teilchenzahl der Probe und \ > 1=200 
=ne ihr Volumen, so hat eine Linie tt 
konstanten Volumens in Abb.9 die Gestalt 
NyEp + Ny €g= v=nNe. (5.2) 720 1 NK (nz,ny)~ 
Entnehmen wir nun zahlreiche Proben 
gleichen Volumens v, so werden diese bei ’ 
guter Mischung ohne Riicksicht auf den yy NM, ; Ks 
Wert von n jedenfalls in dem durch K’ Bee 
und K’’ abgegrenzten Konfidenzbereich 5 INS 
: : ee Ve i 
liegen, der sich nach auSen hin fiir wach- yon 
sende n mit //n erweitert. 80 y BN 702110 
7 
g 
V=0 Ir Ep> Eg 
Ue 
%9 7 7 20 40 60 80 700 
Abb. 8. Zur Veranschaulichung der Volumenschwankung ; Tipr 
Av bzw. Ae in Stichproben S(n) gleicher Teilchenzahl Abb. 9. Zur Beriicksichtigung der verdnderlichen Teilchenzahl bei 
m— 7) tur P — 0540; 10'=— 0,60; Ep/eg = 2; AP = 0,04. Stichproben S(v) gleichen Volumens v=7. 
Im duBersten Falle enthalten diese Proben demnach 
n= n,+n, = n+ An’ bzw. n’ =n, tny=n+An" (5.3) 
Teilchen. Die Wertepaare (nz, n,) und (n{, n,) miissen sowohl die Gleichung (5.2) der Geraden 


v= vals auch die Gleichungen (5.1) der Grenzlinien K’ und K”’ erfiillen. Aus beiden Bedingungen 
findet man nach kurzer Rechnung zur Bestimmung von An’ die Beziehung 


—e — elie 
An’ +9 2— 2 ypqn yi+e =0. (5.4) 


é 


Fiir An’’ hat man nur ¢ durch — ¢ zu ersetzen. 
Die letzte Gleichung wird man (nicht als quadratische Gleichung fiir An’ sondern) iterativ 


durch 


ee ee aed 
A Net oe yt = (5.5) 


é n 
lésen (s= Schrittzahl). Man findet in erster Naherung 


ee et (5.6) 


nm nm é 


n 


406 Stange: Ein Verfahren zur Beurteilung des Giitegrades von Mischungen. Ingenieur-Archiv 


Ein zweiter Schritt ist im allgemeinen nicht notwendig, (wenn man die sonstigen in der Theorie 
enthaltenen idealisierten Voraussetzungen in Betracht zieht). 
Berechnet man nun aus n, und n’ die zugeordnete Teilchenhaufigkeit x’, so erhalt man 


(as n’p telpa)n’ oder mit (5.3) nies preje GAPE : : 
n n n a a An’ 
n 
Mit (5.6) wird daraus 
ae Q pq (1 use? 0 Q Vpq a7 
ed Die renee (5-7) 


Bezeichnet man die Unterschiede zwischen den ,,iuBersten‘ Stichprobenwerten (x’; x’’) und 
dem Sollwert p mit (4x’; 4x’’), so erhalt man aus (5.7) mit (2.6) 
ép— 

2€ 
Fiir Ax’’ hat man @ durch — 9 oder Ap durch — Ap zu ersetzen. 

Das Ergebnis dieser Rechnung kénnen wir folgendermafen zusammenfassen: Entnimmt man 
der Grundgesamtheit Stichproben S(v) gleichen Volumens ¥, so schwankt deren Teilchenzahl n 
zwischen n’ und n’’, Der damit verbundene Ein- 
fluB auf die Teilchenhaufigkeit x ist nach (5.8) 
klein von zweiter Ordnung in (Ap). Solange 
also der als Faktor bei (Ap)? stehende Quotient 
(ep —€9)/2€ keine extrem grofen Werte annimmt, 
geben Stichproben gleichen Volumens nahezu die 

fir §(n)- gleichen Ergebnisse wie Stichproben gleicher 
Abb. 10, Die Schwankung der Volumenhiufigkeit X fiir Proben Leilchenzahl. Rechnen wir (5.8) schlieBlich auf 
S(n) gleicher Teilchenzahl m und ProbenS(v) gleichen Volumensv- die Volumenhaufigkeiten (X, Y) der Proben S(v) 
um, so finden wir, da diese nicht mehr genau 

in dem zugelassenen Toleranzbereich (AP, 4Q), sondern in einem Bereich 


x —p= Ax’ ~Ap+ Ef Ap le. (5.8) 


AP? 


=—AP” <AX= X — P<AP’ (5.9) 
liegen, der nach rechts durch 
ee Ep &o Pa) Ep—€g € 
AP’ =~ = Ap (1 aS Ap) ao Ale (1 1 waar AP] (5.10) 
und nach links durch 
AP" ~ AP (1 a2 Chane AP) (5.11) 
2 Ep£g 


begrenzt wird. Dieser Bereich ist nach Abb. 10 gegeniiber dem urspriinglichen + AP fiir ep > €9 
nach rechts und fiir ep < ég nach links verschoben. Diese Verschiebung ist jedoch nur von 
der GréBenordnung (AP)?. 


Mit diesen Ergebnissen haben wir uns von den eingangs aufgezahlten einschrankenden Voraus- 
setzungen II,, III. und IV. befreit. 


6. Mischung aus drei Stoffen; die Wahrscheinlichkeitsverteilung. Bisher haben wir uns auf 
den einfachen Fall einer Mischung zweier Stoffe (P) und (Q) beschraénkt. Jetzt nehmen wir einen 
dritten Mischungsanteil (R) hinzu. Wir lassen also auch noch die Voraussetzung I. fallen. 
Die vorgeschriebenen Volumenhaufigkeiten der drei Stoffe im Gesamtvolumen seien P, Q und. 
Rmit P+Q+R=1, die zulassigen Toleranzen entsprechend AP, AQ und AR mit AP +AQ 
+AR=0. Die gegebenen Volumenhaufigkeiten P, Q, R, die mittleren Elementarvolumen €ép, 


UG ee i zugeordneten Teilchenhaufigkeiten p, q, r sind jetzt entsprechend zu (1.7) 


ép p=eP, 0 q=é&Q, er tr=éR (6.1) 
miteinander verkniipft, wie man durch eine der friiheren ahnliche Rechnung nachweist. Dabei 
ist 4ahnlich wie in (4.7) die HilfsgréBe @ jetzt durch 


é= pep +qeg+rer oder ky soe oi R 
é€ 


(6.2) 


XX. Band 1952. Stange: Ein Verfahren zur Beurteilung des Giitegrades von Mischungen. 407 


erklart. Diese Gleichungen lassen sich leicht auf mehr als drei Mischungsanteile verallgemeinern. 


Es tritt dann jeweils ein neuer gleich gebauter Summand hinzu. Die Reihenfolge der Stoffe (P), 
(Q) und (R) denken wir so gewahlt, da8 


PROS T (6.3) 
ausfallt, was keine Einschrinkung der Allgemeinheit unserer Untersuchung bedeutet. 


Kine aus n Teilchen bestehende Probe, die fiir die Anteile (P), (Q) und (R) genau die Teilchen- 
zahlen n,, ny und n, besitzt, hat bekanntlich die Wahrscheinlichkeit 


! 

AW(n,, n,, n,3n) = aba Big PO ee ot, T= mt, (6.4) 
Durch (6.4) ist eine dreidimensionale Wahrscheinlichkeitsverteilung im (n,, ny, 1,)-Raum gegeben, 
die jedoch wegen der Zusatzbedingung in (6.4) zu einer ebenen Verteilung entartet, und die allen 
im ersten Oktanten gelegenen Gitterpunkten der Ebene 
n, +ny +n,=7n eine von 0 verschiedene Wahrschein- 
lichkeit zuordnet, wie es Abb. 11 anschaulich zeigt. 


7. Die Ersatzverteilung fiir gro®e Teilchenzahlen. 
Von praktischer Bedeutung ist nun die Tatsache, daB 
auch die hier auftretende trinomische Verteilung, die 
durch Entwicklung des Ausdrucks (p-+q-+r)" ent- 
steht, fir geniigend grofe Teilchenzahl n der Stich- 
probe wieder durch eine Normalverteilung ersetzt 
werden darf. ZweckmaBig fihren wir ebenso wie 


friiher die Teilchenhaufigkeiten der Probe durch 


ny Wy 
? i= ahes a ri (7.1) 
mit 


xtytz=] (7.2) 
nm B Abb. 11. Die ,,raumliche‘‘ Wahrscheinlichkeitsvertcilung 
ein. Dann entsteht fiir geniigend groke Werte von n 4 W(n,, ny, n,3 n) fiir p = 0,25 q=0,35 r=0,5 und n=5, 


tiber der Ebene (7.2) des (x, y, 2)-Raumes eine Normal- 
verteilung mit dem Mittelwert m(p,q,r) und der auf diesen Mittelwert bezogenen sym- 
metrischen Streumatrix 


2 


Cee Tay Ou Ox Ge GyiOng Os Or, 
G(m) =| Oye Oy <0, = (eyorx0ys Oy Oy Oz Oy: 
Ge Gey On OLGx Orel 20, Oy Orr Oz 
(7.3) 
ae (lp) at! = us 
=—| — 9p q(1—4q) = OF = 16", 
1 ONE Sp 7 Lee) 


wobei G* eine Matrix ist, deren Elemente nur von den Teilchenhaufigkeiten p, q, r abhangen. 
Auch hier zeigt man leicht, daB die Determinante der Matrix verschwindet, und dai die Matrix 
selbst den Rang 2 hat, denn die in der linken oberen Ecke stehende zweireihige Determinante 
hat wegen p+q-—+r=1 den Wert pqr= 0. 

Zur weiteren Untersuchung stehen nun wieder die bereits im Abschnitt 2 genannten zwei 
Wege offen: 

a) Entweder projizieren wir alle ,,Gitterpunkte G(x, y,z) der Ebene (7.2) in die ihnen 
entsprechenden Gitterpunkte G’ (x, y) der waagerechten (x, y)-Ebene und ordnen den Punkten G' 
die durch (6.4) gegebene Wahrscheinlichkeit AW zu. Dann geht aber die Symmetrie in den 
Verdnderlichen p, q,1r verloren. ZweckmafBiger ist deshalb die zweite Méglichkeit. 

b) Wir bestimmen die der Streumatrix ©* zugeordneten Hauptachsen durch Lésung der in A 
kubischen Hauptachsengleichung (© = Einheitsmatrix) 


Io* —AG| = 0. (7.4) 


Da die Determinante |6*| der Matrix ©*-verschwindet, so sieht man, daB A=A,=0 eine Lisung 
der Hauptachsengleichung ist. Zur Berechnung der beiden andern Lésungen findet man aus (7.4) 
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nach laingerer Rechnung die quadratische Gleichung 


A? —2(pq + qr +rp)A + 3pqr=0 (7.5) 
mit den Lésungen 
Ae = (pa tar +rp) + V(pq+ar+rp)—3pqr. (7.6) 
Die dem Wert A,= 0 zugeordneten Richtungswerte sind 
uf 
0 PO 0 yee (7.7) 


d. h. es sind die Richtungswerte des Normalenvektors der Stichprobenebene (7.2). Fiir die den 
beiden anderen Lésungen /, und A, zugeordneten Richtungswerte findet man 
a? = f[p(1 — p)A —2pqr], 
B= flq(l —q)A—2pqr], | , (7.8) 
Y=f[r(l—r)A—2paqrl), 
mit 


j22= L (7.9) 


Ley Reise La ye 


Zur Berechnung von (4,, B,,y;) hat man A =A, einzusetzen, und fiir (a5, By, Ya) ist A=A, zu wahlen. 
Da 0, (64, By, 771) und A9(M, By, v2) zwei aufeinander senkrecht stehende, in der Stichprobenebene 
(7.2) gelegene Einheitsvektoren sind, so hat man die Vorzeichen von a, B und y so zu bestimmen, 
daB «+ B+y=0 wird. Ohne die Allgemeingiiltigkeit unserer Untersuchung einzuschranken, 
diirfen wir y>0 wahlen. Dann ist a, in der Stichprobenebene eindeutig festgelegt. GemaB 
Abb. 12 gilt fiir den Winkel y zwischen dem Einheitsvektor e(— y2/2, ob y2/2, 0), der sowohl 
in der (x, y)-Ebene als auch in der Stichprobenebene liegt, und a, die Beziehung 


cos = + ¥2(B, — 1) - (7.10) 


8. Der Zufallsbereich E (J, n) in der Stichprobenebene. 
Die in der Stichprobenebene auf die Achsen a, und a, bezogene Streumatrix wird 


0 ee al 
te ee ’ (8.1) 
BNO eA5 0. as 
Thr ist die ebene Wahrscheinlichkeitsdichte mit dem Differential 
n u \2 v \2 
~|{—=) + (= 
ee : lz) (iz) | a(*) d Ge (8.2) 
4 Vay} \Wae 
zugeordnet. Die Ellipse, welche die Hauptstreuungsmafe 0, = )/A,/n und o, = ya,/n zu Halb- 


achsen hat, wollen wir die Hauptstreuungsellipse H nennen. Man bestitigt leicht, daB das 
Produkt der Hauptstreuungen den Wert 


O71. Op V3pqr/n 
hat, daB also die Flache 0,0, der Ellipse H mit wachsender Teilchenzahl n der Probe wie I/n 
gegen 0 strebt. Grenzt man also um m einen kleinen festen Bereich B ab, der m als inneren 
Punkt enthalt, so fallt mit wachsendem n ein immer gréBerer Anteil von Stichproben in diesen 
Bereich hinein. 
Durch die Substitution 

u v 

atl é und aad (8.3) 
wird die Hauptstreuungsellipse H der (u, v)-Ebene in den Einheitskreis £2 +1? = 1 der (&,7)-Ebene 
iiberfiihrt. Die in dieser Ebene einem beliebigen Kreis K vom Halbmesser 9 zugeordnete Wahr- 
scheinlichkeit ist (wie man leicht nachrechnet) 

1 


—— 


W (BK) = le a. (8.4) 
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_ Daraus erhalt man fiir die im Abschnitt 3 erklarte Irrtumswahrscheinlichkeit den Wert 
ee 51 Qo? 
1—_W=J=e (8.5) 


| Gibt man nun einen festen Wert W der Wahrscheinlichkeit bzw. 1 — W— J der Irvrtums- 
_ wahrscheinlichkeit vor, so 1aBt sich aus (8.5) der zugeordnete Halbmesser o des Kreises K zu 


oJ) = —2in J 


(8.6) 


_ berechnen (Abb. 13). Je kleiner J vorgeschrieben wird, um so gré®er mu 9 gewahlt werden. 


| 
30 
Stichprobenebene 
G5 
W:19? 
40, — 
J10" 0) 95 700 
70 5 oO 
ye" Abb. 13. Der Zusammenhang zwischen der Irrtums- 
E wahrscheinlichkeit J und dem Halbmesser Q eines 
Abb. 12. Die Stichprobenebene mit den Hauptrichtungen qa, und a). Kreises, der (W-100)% der Beobachtungen enthalt. 


Dem Kreis K(g) der (&, 1)-Ebene entspricht in der Stichprobenebene (7.2) eine Ellipse E mit 


dem Mittelpunkt m/(p, q, r) und den Halbachsen 


1 tb ae 
eal und = Vag 


(8.7) 


in Richtung der Einheitsvektoren a, und a,. Im Innern dieser Ellipse E liegen bei der zugrunde- 
gelegten Irrtumswahrscheinlichkeit J bei ,,guter Durchmischung* der Stoffe W-10? (%) aller 
Stichprobenpunkte (x, y,z). Der so festgelegte und durch eine Linie gleicher Wahrschein- 
lichkeitsdichte begrenzte Bereich E ist bei festgehaltener Irrtumswahrscheinlichkeit J = konst. 
der kleinste in der Stichprobenebene mégliche Bereich. Er umfaBt diejenigen Stichproben- 
punkte (x, y, z), welche die gréBte Wahrscheinlichkeitsdichte besitzen. Allen Punkten auBerhalb 


von E ist eine kleinere Wahrscheinlichkeitsdichte zugeordnet. 


Kehren wir noch einmal zu den Gleichungen (8.7) zuriick. Der erste Faktor 1/yn ist nur von 
der Teilchenzahl n der Stichprobe abhangig. Der zweite Faktor hangt nur von den Teilchen- 
haufigkeiten (p, q, r) ab und bleibt bei zyklischer Vertauschung von p, q, r ungedndert, da auch 
die Lésungen A, und A, der quadratischen Gleichung (7.5) dabei fest bleiben. Der dritte Faktor @ 
wird nach (8.6) allein durch die zugrundegelegte Irrtumswahrscheinlichkeit J = 1— W bestimmt. 


9. Die praktische Form der Lésung. Mit Hilfe der vorausgehenden Ergebnisse lassen sich die 
Bereiche E in der Stichprobenebene ein fiir allemal festlegen. Man setzt in (8.7) willkilich 
n=1 und 9=1 und kann dann die fiir diesen Sonderfall entstehende Ellipse E* mit den Halb- 


achsen 


= VA und OL Vrs 


(9.1) 


in den Richtungen a, und @, einzeichnen. Die einem belichigen Wertepaar (n;Q) zugeordnete 


Ellipse E findet man dann einfach durch eine Mafstabsinderung. 
In diesem Zusammenhang sind zwei Bemerkungen niitzlich. 
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a) Da es auf die Reihenfolge der Mischungskomponenten nicht ankommt, so sind die 3!=6 


Permutationen der Elemente p, q, r einander gleichwertig. 


Es gehéren also in dem Merkmal- 


dreieck der Ebene x-++-y +-z=1 gemaf Abb. 14 je 6 Punkte zusammen, welche durch Spiegelung 


a Q 
i 


nN 
We 


/) INAS 
VAN AV: 
CIA 


ok 


ADK 
IBID 


(Ze 


an den drei Mittelloten p= q, gq=r, 
r= p des Dreiecks auseinander her- 
vorgehen. Man beherrscht also die 
ganze Verteilung, wenn man sie 
in einem der Teildreiecke, etwa in 
dem geschrafften, kennt, was fiir 
die numerische Auswertung der 
Formeln von Vorteil ist. 

b) Dividiert man _ schlieBlich 
noch alle Abmessungen in der Stich- 


probenebene durch /2, so gelangt 
man zu der bekannten Darstellung 
durch Dreieckskoordinaten in Netz- 
papier, wie es bei der Analyse von 
Dreistoffgemischen Verwendung 
findet. In dem Falle ist die Seiten- 
lange des Stichprobendreiecks zu 1 
normiert. 


Die folgende Zahlentafel 1 und 


‘Abb. 15 enthalten schlieBlich fiir 


einige Werte (p, q, r) die Lange und 
Richtung der Hauptachsen!. Will 
man die zu bestimmten Werten der 


Abb. 15. Richtung und Lange der Hauptachsen der Streuungsellipsen im Stichprobendreieck der Ebene x + Yes 
(die Seitenlange des Dreiecks ist auf den Wert 1 normiert). 


Irrtumswahrscheinlichkeit J und der Teilchenzahl n gehérende Zufalls-Ellipse E(J; n) in 
Dreieckspapier einzeichnen, so hat man wegen der eben genannten Normierung die Werte 


* Fiir die gewissenhafte Durchfiihrung der umfangreichen numerischen Rechnungen habe ich meiner 


Mitarbeiterin Fraulein Dipl.-Math. E, Steinmetz zu danken. 
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o* der Zahlentafel abweichend von (8.7) mit dem Faktor o(J)/y2n zu multiplizieren. Die mit 
Hilfe des Randmafistabes aus Abb.15 abgelesenen Lingen der Halbachsen entsprechen, wenn 
man sie mit 10 multipliziert, den Werten o* der Tabelle 1. 


Tabeile 1, Mi 
Pp q r oy =a, ves = Va, 
0,1 0,1 0.8 0.490 0,216 
1 22 7 0,578 0,354 
1 3 6 0.638 0,264 
1 4 5 0,667 0,367 
2, 2 6 0,600 0,447 
2s 3 5 0,623 0,482 
2 4 4. 0,633 0,490 
0,3 0,3 0.4 0,600 0,548 


Die Ermittlung der notwen- 
digen Teilchenzahl n—7ni fiir die 
Stichprobe S(n) geht nun fol- 
gendermafen vor sich. Zu den 
vorgeschriebenen Volumenhaufig- 
keiten P,Q, R und den zugelas- 
senen Toleranzbereichen AP, AQ, 
AR berechnet man mit (6.1) und 
(6.2) die Teilchenhaufigkeiten p, 
q, r und die ihnen zugeordneten 
Toleranzbereiche Ap, Aq, Ar. Der 
Toleranzbereich sei in der (X,Y, 
Z)-Ebene X+Y+Z=1 durch 
ein Parallelogramm ABCD ge- 
geben, welches den Mittelpunkt 
M(P,Q, R) besitzt, und dessen 
Seiten durch die beiden (unab- 
hangig voneinander wahlbaren) 
Fehlergrenzen 4P und AQ fest- 
gelegt sind. Die Seiten des Par- 
allelogramms laufen zu zwei 
Grundrichtungen parallel (Abb. 
16). Die vier Eckpunkte bildet 
man am einfachsten mit Hilfe 


der Gleichungen 


pe NPP AG. 4 (9:2) 


1 
Ep 


mit 


— de = (2 4 ah a 


Ep 70) ER 


auf die (x,y,z)-Ebene x+y +2 
=1 ab. Entsprechende Formeln 
gelten fiir Aq und Ar. Bei dieser Abb. 16. Zur Berecbnung der Toleranzbereiche abed baw. a’b’c’d’ im Stichproben- 
Umrechnung geht das Parallelo- dreieck der (¥,y, «)-Ebene. 
gramm ABCD wieder in ein 
Parallelogramm abcd mit dem Mittelpunkt m/(p,q,r) der Stichprobenebene itber. Seine Seiten 
laufen aber nicht mehr zu den Grundrichtungen parallel. 

Nun zeichnet man in dem Merkmaldreieck um den Punkt m(p,q,r) fiir eine ganz beliebige 
Teilchenzahl n eine Streuungsellipse E mit den Halbachsen oj o/f2n und of o/f2n. Dieser 
Ellipse E sind Toleranzbereiche Ax, Ay, Az zugeordnet, die man nach Abb. 17 anschaulich findet, 


indem man an E die zu den Linien p= konst., g=konst., r= konst. parallelen Tangenten zeichnet. 
29 
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FormelmaBig sind sie durch 
Ax="- Ippo)» dy= = Val —gelJ) z= yr —r old) (9.3) 


gegeben. Man wahlt nun n=7 so groB, daB die der Teilchenzahl n zugeordnete Ellipse E ganz 
in dem durch abcd gegebenen zulassigen Toleranzbereich liegt. Das ist stets méglich, da sich 
die Hauptachsen der Ellipse bei wachsender Teilchenzahl mit I/yn zusammenziehen. 

Ein bequemes (hinreichendes, aber durchaus nicht notwendiges) Verfahren, um E in das 
Parallelogramm abcd hineinzuzwingen, geht folgenden Weg. Man bestimmt nach Abb. 16 
das gréBte, ganz im Innern von abcd gelegene Parallelogramm a’b’c'd’, dessen Seiten zu zwei 
Grundrichtungen parallel laufen. Diesem Bereich sind die Abweichungen Ap’ < Ap, Aq’ < Aq, 
|Ar’|<|Ar| zugeordnet. Halt man 
nun die Bedingungen 


[As] <Apee Aye Age 
|Az| < |Ar’| 


ein, so liegt die Ellipse FE sicher im 
Bereich abcd. Die dazu notwendige 
Teilchenzahl n=n findet man ge- 
ma (9.3) aus den drei Ungleichungen 


=e 20) aN (9.4) 


indem man daraus den gréBten Wert 
fiir n auswahlt. Das Volumen v der 
Probe ist dann v=v=nNneé. 

Im folgenden ist der Gang der 
Rechnung mit den wichtigsten Zah- 


lenwerten fiir ein Beispiel zusammen- 

gestellt, wobei wir des Zusammen- 

hanges wegen einige Ergebnisse des 

aie LW) £  niachsten Abschnitts vorwegnehmen. 

a7 02 3 O# : ‘ ; 

eAy - Die Rechnung und die Abbildun- 

Abb. 17. Elliptischer Zufallsbereich E, (Jn) fir P:Q:R = 0,2: 0,3: 0,5; gen 18 und 19 zeigen, da man auf 

AP = 0,015; AQ = 0,023; |AR| = 0,038; J = 0,01 und ny = 200. Mit wache dem beschriebenen Wege die zu- 
sendem n zieht sich die Ellipse ahnlich zusammen, bis sie im Parallelo- e 

gramm abed bzw. a’b’c’d’ liegt. gelassenen Toleranzbereiche AP, 


AQ und |A R| fiir die Anteile (P), (Q) 
und (R) gar nicht ausnutzt, da erstens das Parallelogramm a’b’c'd’ ganz im Bereich abcd liegt, 
und da zweitens von den drei nach (9.4) méglichen Anzahlen n,, n, und n, die gréBte aus- 
gewahlt werden mu&. Die Bedingungen |AX| < AP, |AY| < AQ und |4Z| < |AR| werden da- 
mit gewissermafen alle ,,iibererfiillt*, wie Abb. 19 zeigt. 

Nachdem n bekannt ist, erscheint es deshalb zweckmabig, zu n und (p;q;r) durch zuriick- 
laufende Rechnung die durch % gewahrleisteten Toleranzen (Ap, dq, Ar) und (4 P, AQ, AR) 
zu ermitteln. Man findet dann die in den Zeilen 5 und 6 des Beispiels zusammengestellten Werte, 
wobei wegen n= it, hier natirlich Ap = Ap’ sein muB. Da8 die Summen Ap + Aq — |Ar| und 
AP +AQ—|AR| nicht verschwinden, liegt daran, da® im elliptischen Zufallsbereich die 
auBersten Werte Ap und Ag bzw. AP und AQ gar nicht gemeinsam auftreten. Infolgedessen 
ist sowohl |Ar| < Ap + Aq, als auch |AR| < AP + AQ. 


10. Stichproben S(v) gleichen Volumens. Die vorausgehenden Betrachtungen gelten fiir Stich- 
proben S(n) gleicher Teilchenzah]l. Wir wollen ahnlich wie im Abschnitt 5 abschatzen, wie gro8B 
der Unterschied An=n—vn der Teilchenzahl ausfallt, wenn wir Stichproben S(v) gleichen 
Volumens v= v entnehmen. 


413 


Stange: Ein Verfahren zur Beurteilung des Giitegrades von Mischungen. 


XX. Band 1952. 


‘uIa1EqN usMMOYTOA qostyyeid ospe meus eyorerag s1q ‘3][93S03 
-1ep Weyqieqn qoejyuyez puis uasunyqowaqy wazjperqosas sq ‘jaidsraqualqez 
urs my Z‘A ‘XK Yoorwspueqoidygonsg wr susunjoa ueyores (2)g¢ naeqosdyos 
pun [yezueqpey ryore[s (u)g ueqoidqgong my oqorereqsy[eynz 1G “61 “GqV 
Ove 
20-0 


520 


Z0 Hf 0 
T 


‘ala1aqN wsMIMOY]OA yosmyvid osje usMUMIYs eYoIIIVg aIq “31]93803 
-1ep IEYIeGN yYyoeyuyez puts useZunyormqy ueisyeryoses eid ‘jerdstequo]qez 
ula Ing 2% 6K Se yoorsrpueqoid gong Wi susmIN[OA Ueyqo1s[s (a)y usqoidqong 
pun [qezueyopeay, rzeyore]s (u)g ueqoidqong imy eyerereqsy[eynz Aq “gl ‘qqV 


Lp 


£20 0=2 wiypuldp_ 220 420 b 
T T 


7) 


a 


29* 


414 Stange: Ein Verfahren zur Beurteilung des Giitegrades von Mischungen. Ingenieur-Archiv 


Tabelle 2. (Z bedeutet Zeichnung.) 


| Beispiel (P) (0) (R) RGer 


1 Gegeben sind die Volumenhaufigkeiten P=0,20 Q=0,30 R=0550 
mit den zulassigen Toleranzbereichen AP=15-10 | -4O=23-10-* ~|, |AR| = 38-10 
und den Elementarvolumen Ep =3 ég=2 Ep=1 

2 Das mittlere Volumen wird € =1,395 (6.2) 

3 Dazu gehoren die Teilchenhaufigkeiten p=9,093 q=0,209 r=0,698 (6.1) 
mit dem zulassigen Toleranzparallelo- 
gramm abcd Ap 9,81 0= 8a Agi —22.5- 0m Arma, Ll mo ime oee) 
Diesem ist das Parallelogramm a’ b’c’ d’ 
mit den zulassigen Toleranzbereichen Ap’ ==6,8-107* | Aq’ =16,5°10~ | \Art 2375310 Z 
zugeordnet. 

4 Fir W= 0,99 erhaélt man die not- Sa 16,8 - 10° n, = 5,6 -10° n, = 3,6-10° 
wendige Teilchenzahl n und das Vo- (9.4) 
lumen v der Stichprobe n=16,8-103 » = 23,4-103 : 


5 Die durch n gewahrleisteten Toleranz- Ap =6.8-10-8 Aq=9,6-10-3 | Ar=10,7-10-3 | (9.4) 
bereiche sind in der (x, y, z)-Ebene 


6 | baw. in der (X, Y, Z)-Ebene AP=13,1-10-3 | AQ=12,5-10-3 | AR=12,9-10-3]  Z 


7 Mit ¢,.,, wird die gréBte Schwankung & = 0.659 (10.11) 
der Teilchenzah] (4n/n) in Proben ig, 


max 


S(v) gleichen Volumens v = v (An/M) ax = 141-107 ae 
8 Dadurch andern sich die Toleranz- Ties aaa cee 
bereiche fiir die aus Proben S(v) be- de Roy: adh He NGO 
rechneten Teilchenhaufigkeiten Ap 055% Aq 05% Ar = 0,075 | Ce 
(x,y,z) um héchstens 
9 Die entsprechenden Anderungen der OD Wty; Tip 
Toleranzbereiche fiir die Volumen- A*P =0,4% AQ ==, 3°/4 A*R =0,6%, Yi 
haufigkeiten in der (X, Y, Z)-Ebene AP AQ AR 
sind 


Der Konfidenzbereich fiir die Teilchenzahlen n,, n, und n, bei festen Werten p, q, rin Ab- 


hangigkeit von nist ein réhrenférmiges Gebilde. Seine Mittellinie ist die Gerade der Erwartungs- 
werte 


aot (6 Ny =1q, =i) (10.1) 
Die Ebene n, +-ny --n,=n schneidet die Rohre in einer Ellipse, deren Halbachsen 
a, = in Vay 0 und i — jn Vas 0 (10.2) 


parallel zu den aus (7.8) berechneten Kinheitsvektoren @, (a, B,, 1) und (a, By, yy) liegen. Die 
mit wachsendem n in der Schar paralleler Ebenen entstehende Ellipsenschar ist bei gegebenen 
Werten (p,q, r) demnach fest im Raume orientiert. 


Bezeichnen wir mit §{ den Ortsvektor vom Ursprung zu einem Punkt (n,, ny, n,) der Roéhren- 


flache und mit c(p, q, r) einen Vektor mit den Komponenten p, q und r, so ist die Gleichung der 
Rohrenflache 


R= ne + yn e(oF cos yp a,+o% sin y a,) , (10.3) 


wobei n und p die Flachenparameter darstellen. Haben wir nun die zu den vorgeschriebenen 
Volumenhaufigkeiten (P,Q, R), den zugelassenen Toleranzen (AP, AQ, AR) und der Irrtums- 
wahrscheinlichkeit J gehérende notwendige Teilchenzahl n der Probe (nach Abschnitt 9) be- 
rechnet, so lautet die Ebene, welche die Proben S(v) gleichen Volumens v= 7% enthiilt, 


N,Ep + Nyég + N,Eg = V = NE. (10.4) 
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Eine kurze Rechnung zeigt, daB ein Wertetripel (n,, n,, n.), welches der Rohrenflache (10.3) 
und der Ebene (10.4) gleichzeitig angehort, der Gleichung 


. né+oyn (Aep+ Beg+ Cer) =e (10.5) 
geniigt, wobei 
A(y) = of a, cos y + oF a, sin p, 
B(y) = of B, cos p + oF By sin p, (10.6) 
C(y) = of y, cos p + 02 Vp sin | 


eS ok spater bemerken wir, daB die gré®ten Werte von A, B und C gema8 (7.6) und (7.8) 
ure 


Ave a yp(l =P) ? Bixaz ad Va —q) > Cnae — yr(l —r) (10.7) 
gegeben sind. Mit dn—n—v7n erhalt man aus (10.5) als Bestimmungsgleichung fiir An 
Aép+ Bengt Ce =~ i Ale 
A P Q Ree / An 
n+ : eyn 1+ > =0. (10.8) 
In erster Naherung ist 
An Atp+ Beg+Cep 1 
sat . o(J) ae (10.9) 


Die relative Anderung der Teilchenzah] beim Ubergang von Proben S(n) fester Teilchenzahl 


n= n zu Proben S(v) festen Volumens v= 1 ist von der Gréfenordnung Lyn. 
Die gré8te mégliche Anderung An,,,, erhalt man, wenn man in (10.9) den Ausdruck 


e(y) = Aep+ Beg+ Cer 
gemaf (10.6) als Funktion von yp betrachtet. Aus 
€(y) = of (0; ep + By eg +71 Er) cos p + oF (A. Ep + Bo eo+ 72 €r) Sin y (10.10) 
findet man mit (7.6) und (7.8) nach langerer Rechnung 
Emax = Ep P(1 — p) + €8 q(1 — q) ter (1 —r) —2(ep eo pa tegen gr +érep rp). (10.11) 


Damit gilt schlieBlich fiir die gesuchte relative Schwankung der Teilchenzahlen 


™ éE 


Mit Hilfe dieses Ergebnisses sind wir nun imstande, den Einflu8 der veranderlichen Teilchen- 
zahl bei der Berechnung der Teilchenhaufigkeiten abzuschatzen. Bei Proben S(n) gleicher 
Teilchenzahl n=n gilt in der Ebene n=7n auf dem Rande des Konfidenzbereichs fiir die 
Teilchenhaufigkeiten, z. B. fiir x, 
me "Pt Vn eAly) 


i oe 


n 1m 


oder fiir den Unterschied gegeniiber dem Ex wartungswert p, 
1 
ae (4x), == eA(y) eae 
\ ri 


Fiir Proben S(v) gleichen Volumens v= 1% hat man in der Ebene v= v auf dem Rande des Kon- 
fidenzbereichs entsprechend 


(x —p), (10.13) 


nm, _np+\n @A(y) 


) 


a 


also mit n=n-+An, 


(Ax) ~ (4%), =(1 peak (10.14) 
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Verwendet man also an Stelle der eigentlich notwendigen Proben gleicher Teilchenzah] solche 
gleichen Volumens, so wird der gréBte mégliche Unterschied fiir die aus diesen Proben berech- 
neten Teilchenhaufigkeiten mit (10.7), (10.12) und (10.13) dem Betrage nach 

| 1 oe pera s 2 10.15) 
ae ae yp(l P) iS DS Pot ( : ) 
Entsprechende Gleichungen gelten fiir Ay und Az. Man hat nur p durch q bzw. r zu ersetzen. 
Der Ubergang zu Proben S(v) ist demnach erlaubt. Sie liefern in erster Naherung die gleichen 
Konfidenzbereiche fir die Teilchenhaufigkeiten x, y, z wie Proben gleicher Teilchenzahl. 

Es ist niitzlich, aus (10.15) den Sonderfall von nur zwei Komponenten (P) und (Q) in der 
Mischung herzuleiten. Aus (10.11) folgt mit r=0 zunachst 


— (Ax) 


Emax = (Ep — €0)*PQ- ‘ 
Aus (10.15) erhalt man damit 


| (43), See 


ew =n 


Ep—eE 1. ep 
a eel pgs eee oe ol (py. 


Das stimmt mit (5.8) tiberein, wenn man beriicksichtigt, daB dort x’ — p= (Ax), _— darstellt 


und da8 Ap damals durch (2.6) mit 9 und ” verkniipft war. Damit ist der Zusammenhang mit 
den fritheren Ergebnissen hergestellt. 

11. SchluBbemerkungen; Zusammenfassung. Die vorausgehenden Untersuchungen wurden 
mit Hilfe der Begriffe ,, Volumenhaufigkeit P,Q, R und ,,Teilchenhaufigkeit“ p, q, r entwickelt. 
Der Gedankengang bestand im wesentlichen aus drei Schritten. Die Volumenhaufigkeiten P,Q, R 
wurden in Teilchenhaufigkeiten p, q, r umgerechnet. Fiir p, q, r bzw. die Einzelteilchen wird bei 
der Entstehung der Mischung die Giiltigkeit statistischer GesetzmaBigkeiten vorausgesetzt. Die 
Ergebnisse gelten dann zunachst nur fiir Stichproben S(n) gleicher Teilchenzahl. SchlieBlich 
wird gezeigt, daB die Ergebnisse der Teilchen-Theorie cum grano salis auch fiir Stichproben S(v) 
gleichen Volumens Giiltigkeit besitzen. Es wurde versucht, die wichtigsten Ergebnisse in einer 
Form darzustellen, die wohl ohne wesentliche Anderungen auch praktischen Bediirfnissen gerecht 
wird, wenn sich herausstellt, daB die Theorie zur Beschreibung und Beurteilung wirklicher Misch- 
vorgange brauchbar ist. Ein ganz anderes, hier nicht beriihrtes Problem besteht darin, wie man 
aus einer Stichprobe mit den Teilchenhaufigkeiten x, y, z bei unbekannten Anteilen P, Q, R der 
Gesamtmischung einen RiickschluB von x, y, 2 iiber X, Y, Z auf diese unbekannten Haufigkeiten 
machen kann. 

Selbstverstandlich bleiben alle Gedankenginge vollig ungedndert, wenn man an Stelle von 
P,Q, R die relativen ,,Gewichtshaufigkeiten“ PO R, mit iP +Q, +R,=1 zugrundelegt. 
Die Elementarvolumen ép, ég, €g sind dann durch die Elementargewichte lp, Mg, 4x der Einzel- 
teilchen zu ersetzen. Es bietet keine Schwierigkeiten, diese Gewichtshaufigkeiten mit den 
Teilchenhaufigkeiten und den Volumenhaufigkeiten zu verkniipfen. In manchen Fallen scheint 
eine Kennzeichnung der Mischung durch Gewichtsanteile sogar sinnvoller zu sein. 

Fiir kleine Werte der Teilchenzahl n verlieren die Ergebnisse ihren Sinn, weil der Ersatz der 
diskreten Multinomialverteilung durch eine Normalverteilung unzulassig wird, und weil sich 
die statistischen GesetzmaBigkeiten erst bei groBen Teilchenzahlen auswirken kénnen. Diese 
Gefahr der zu geringen Anzahl wird aber bei praktisch vorkommenden Mischungen wohl kaum 
bestehen. 

Der EinfluB des zwischen den Einzelteilchen bestehenden leeren Volumens ist nicht beriick- 
sichtigt worden. Da dieser Einflu8 nicht nur von den Korngré8en ep, €g, €r, sondern haupt- 
sichlich von der Gestalt der Kérner abhangt, wird man zu Korrekturwerten Zuflucht nehmen 
miissen, die man wohl nur durch einen Versuch im konkreten Falle bestimmen kann. 

Es ist schlieBlich selbstverstandlich, daB die Uberlegungen auf Vierstoffgemische und noch 
weiter verallgemeinert werden kénnen. Gedanklich ist dann im Vergleich zu der vorausgehenden 
Entwicklung nur wenig Mehrarbeit zu leisten, wahrend der mathematische Aufwand._natiirlich 
rasch anwachst und man fir fiinf und mehr Mischungskomponenten nicht mehr die Méglichkeit 
einer anschaulichen Darstellung hat. Der Fall von vier Komponenten P, Q, R, T liegt in dieser 
Hinsicht gerade an der Grenze. Dadurch, da® man den vierdimensionalen (X, Y, Z, T)-Raum 
(bzw. den x, y,2,t-Raum) durch die ,,Ebenes X+Y+1+Z2+7T=—1 (bzw. x +y¥ +2-+t=1) 
schneidet, entsteht bei geeigneter orthogonaler Substitution und passender Normierung ein 
dreidimensionales regelmaBiges Tetraeder der Seitenlange 1. Dieses Stichprobentetraeder ent- 
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spricht dem Stichprobendreieck aus Abbildung 14. Die Punkte M (bzw. m) im Innern dieses 
Tetraeders bilden die Gesamtheit der méglichen Haufigkeitskombinationen X, Y, Z, T’ bzw. 


x, y,%,t. Mit Hilfe von Tetraederkoordinaten ]48t sich dann fiir vier Mischungskomponenten 


eine 4hnliche Darstellung erzielen wie bei einer Dreikomponentenmischung im Dreiecknetz der 
Abbildung 14. 


Wenn die Mischung auf dem in Abbildung 1 dargestellten Wege entsteht, gibt die 
Theorie ein gutes Abbild der Wirklichkeit, wie das folgende Beispiel zeigt. Es wurde eine aus 
N=4-10! Teilchen bestehende Mischung aus zwei Komponenten hergestellt. Als Zufallsvorgang 


% ASummenhdutigheit 
99,398 


999 


i 9 Vo-Bereich cal | 
Zahl der Teilchen (P) in der Mischung 
ae ee Me 
sa 7 Ti om) 35 60 


Abb. 20. Verteilung von 400 Stichproben aus einer Zufallsmischung mit N = 4-104 Teilchen im Wahrscheinlichkeitsnetz 
nach Rechnung (Gerade) und Versuch (Eiuzelpunkte) fir p = q = 0,5 und n = 100. 


dienten die Ergebnisse des Roulettespieles der Spielbank in Baden-Baden, indem man die 
Teilchenhaufigkeiten p und q mit den Wahrscheinlichkeiten fiir ,,rot** und schwarz gleichsetzte. 
Die Mischung wurde in 400 Stichproben S(n) von je n = 100 Teilchen zerlegt. Die Stichproben 
wurden dann auf ihre Zusammensetzung untersucht. Das Exgebnis dieses Versuchs im Vergleich 
zu dem nach der Theorie zu erwartenden Verlauf zeigt Abbildung 20. Natiirlich ist die gute Uber- 
einstimmung noch kein Beweis dafiir, dafB das mathematische Modell auch die wesentlichen Ziige 
einer Mischung wiedergibt, die auf ganz andere Weise in einem Mischgerat erzeugt wurde. Die 
Nachpriifung, wieweit die entwickelte Theorie bei praktisch auftretenden Mischungsvorgangen 
ein den Tatsachen entsprechendes Bild liefert, soll weiteren Untersuchungen vorbehalten bleiben. 


(Eingegangen am 1]. Juni 1952.) 


Anschrift des Verfassers: Dr. K. Stange, Karlsruhe, Technische Hochschule, Institut fir Mathematik 
und ihre techn. Anwendungen. 
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Experimentelle Untersuchungen am schiebenden Fliigel* 
Von W. Jacobs. 


1, Einleitung. Bei allen bisherigen Messungen an Pfeilfliigeln itherdecken sich grundsatzlich 
drei verschiedene Einfliisse, namlich 
1. der EinfluB des schiebenden Fliigels, 
2. die besonderen Verhaltnisse in der Fliigelmitte, wo zwei Fliigel mit verschiedenem Schiebe- 
winkel zusammentreffen, 
3. die Umstrémung der Fliigelenden. 
Um iiber Pfeilfliigelfragen weiter Klarheit zu erhalten, erscheint es notwendig, diese Einfliisse 
getrennt zu behandeln. Die folgenden Versuche haben den Zweck, lediglich den EinfluB des 
schiebenden Fliigels naher zu untersuchen. 
Um die beiden anderen Einfliisse auszuschalten, wurden Fliigel, mit den Schiebewinkeln 
B = 0°, 30° und 45° mit Endscheiben versehen (Abb. 1, Anordnung I). Diese Messungen bilden 
somit einen Beitrag zum schiebenden Fliigel unendlicher Spannweite (,,quasiebenes*‘ Problem). 
Zu Vergleichszwecken wurden die 
| q ee Fligel 6 =0° und 45° auch ohne 
Endscheiben vermessen. Diese An- 
g ordnung kann dazu dienen, bei einem 
schiebenden Fligel den Einflu8 der 
Randumstrémung (gestrichelte Fla- 
chen der Anordnung II in Abb. 1) zu 
untersuchen. 

Bei den theoretischen Uberlegun- 
gen zum Pfeilfliigelproblem spielt das 
cos f-Gesetz fiir den Auftrieb eine 

Abb. 1. Erlauterungsskizze: Anordnung I: Schiebender Fliigel zwischen End- prrelnee ea ee Pte Were a 
scheiben; Anordnung II; Schicbender gerader Fliigel; Anordnung ITI: Pfeilfligel- unter, daB beim Ubergang vom nicht 
schiebenden Fliigel auf einen Fligel 
mit dem Schiebewinkel f der Auftrieb mit cos f kleiner wird, wenn dabei das Fliigelprofil und 
der Anstellwinkel in der Anstrémungsrichtung konstant bleiben. Dieses cos B-Gesetz soll durch 
unsere Messungen experimentell gepriift werden. Es wird sich ergeben, daB es bis zu grofen 
Schiebewinkeln sehr gut erfiillt ist. 
Diese Messungen bilden einen weiteren Beitrag zu den systematischen Pfeilfliigelmessungen 
des Aerodynamischen Institutes der T. H. Braunschweig 1 2. 


2. Bezeichnungen. Im folgenden sind die Bezeichnungen nach dem Normblatt DIN L 100 


benutzt worden, namlich 

XY, ® Achsen des Koordinatensystems (der Koordinatenursprung liegt im 1/4-Punkt der 
Profilsehne im Mittelschnitt des Fligels, Abb.1) 
Fligeltiefe in Anstrémungsrichtung, 


l 

b Fliigelspannweite senkrecht zur Anstrémungsrichtung, 
F Fligelflache, 
A 
a 


iy) 


= b?/F Seitenverhiltnis, 
Anstellwinkel des Fliigels, 


* Die hier mitgeteilten sehr einfachen Messungen, die im Jahre 1944 im fritheren Aerodynamischen 
Institut der Technischen Hochschule Braunschweig (Leitung: Prof. Dr. H. Schlichting) duseetal ee Beagle a! 
haben eine wichtige grundsatzliche Bedeutung fiir den Pfeilfliigel. Trotz der in der Zwischenzeit stark 
ausgedehnten Forschung iiber Pfeilfliigel sind systematische Messungen der hier vorliegenden Art noch 
nicht bekannt geworden, so da eine Veroffentlichung unserer Ergebnisse auch heute noch angebracht er- 
eee Dies ergab sich insbesondere aus einer Diskussion auf der Tagung des VDI-Fachausschusses fiir 

tromungsforschung in Berlin am 7. 2. 1952. Die Bearbeitung des Manuskriptes fiir die Veréffentlichung 
Bowie die theoretischen Uberlegungen des Abschnittes 5 stammen von E. Truckenbrodt 
: W. Jacobs, Ing.-Arch, 18 (1950), S. 344; 19 (1951), S. 83 und S, 400. 
H. Trienes und E. Truckenbrodt, Ing.-Arch. 20 (1952), S. 26, 
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a” unkorrigierter geometrischer Anstellwinkel, 
p Schiebewinkel des Fliigels, 
Gq + v? Staudruck der Anstrémung, 
DaiPs gemessener Druck am Profil, 
A : ; 
— an Auftriebsbeiwert des ganzen Fliigels, 
0 
W : : 
Pa aa Widerstandsbeiwert des ganzen Fliigels, 
. 
On = qF 1 Momentenbeiwert des ganzen Fliigels, schwanzlastig positiv (Momentenbezugsachse 
0 
die y-Achse), 
Ca srtl értlicher Auf- 
triebsbeiwert, 
Cosril értlicher Wider- 
standsbeiwert. 


3. Modellbeschreibung. 
Fiir die durchgefiihrten Un- 
tersuchungen sind drei Flii- 
gel konstanter Tiefe mit den 
Schiebewinkeln 6 = 0°, 30° 
und 45° verwendet worden 
(Abb.2). Die Spannweite der 
Fliigel betragt b = 0,60 m, 
die Tiefe 1 = 0,20 m. Daraus 
ergibt sich ein Seitenver- 
haltnis von A =3. Als Profil 
wurde NACA 0015 gewahlt. 
Langs der Spannweite wa- 
ren 9 DruckmeSschnitte mit 
je 18 DruckmeBstellen vor- 
handen, deren Lagen aus 
Abb.2 zu ersehen sind. Die 
Fliigel waren mit Endschei- 
ben von rechteckigem Um- 
riB (0,35 x 0,60 m) versehen 

und aus PrefSschichtholz 
 hergestellt. 


4, Umfang und Ergeb- 
nisse der Messungen. Die 
Messungen wurden noch vor 
Kriegsende in dem fritheren 
Aerodynamischen Institut 
der Technischen Hochschule 
zu Braunschweig durch- 
gefihbrt. Der Strahldurch- 
messer des Kanales war 
1,2m. Die Geschwindigkeit 
betrug v— 40 m/sec, was 
einer Reynoldsschen Zahl 
Re = vl/y = (5,6 + 10° . ent- 
spricht. Es wurden folgen- 
de Messungen ausgefiihrt: 

a) Dreikomponen- 
tenmessungen (C,, Cy», €m) 


Endsche/be 
ty 
| Momentenbezugsachse 
Flage/ 1 
=80° 
Momentenbezugsachse 
Mittelschnitt 
2 

Fliigel Ws es 
p=? ®& Z 


aa 


0 20 “| 
Protil NACA 0075 


Abb. 2, Ubersichtsskizze der verwendeten Fliigel und Lage der Mefschnitte und 
MeSanhokrungen. MaBe in mm. 


fiir alle drei Fliigel mit Endscheiben und fiir die Fliigel 6 = 0° und B = 45° ohne Endscheiben. 
Die Polaren, Auftricbs- und Momentenkurven sind in Abb.3 und 4 dargestellt. 
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1a 


b) Druckverteilungsmessungen langs Spannweite bei 6 Anstellwinkeln (a’ = 0, 3, 6, 9, 
12, 15°). Auf die Wiedergabe simtlicher Druckverteilungen an den einzelnen Fligelschnitten 
wird verzichtet. In Abb.5 sind die Verteilungen im Mittelschnitt (Nr. V) fiir die Fliigel mit End- 
scheiben bei drei verschiedenen Anstellwinkeln zusammengestellt. Die aus den Druckverteilungen 
ermittelten ¢,¢-Verteilungen sind in Abb.7 und 8 fiir die Fliigel mit und ohne Endscheiben 
dargestellt. rw 

c) Impulsmessungen an verschiedenen Stellen langs Spannweite fiir verschiedene Anstell- 
winkel zur Bestimmung der drtlichen Profilwiderstandsbeiwerte (Abb. 9). snl 

d) Messungen der Grenzschicktdicke an der Hinterkante des Fligels 6 = 45° (mit 


Endscheiben) bei verschiedenen y-Werten (Abb.10). 
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Cn» VonEndtscheiben| | G2 Go a tl 47 Ge 3 Gy 
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Abb, 3. Polaren, Auftriebs- und Momentenkurven fiir Abb. 4: Polaren, Auftriebs- und Momentenkurven fiir die 
die Fligel B = 0°, 30°, 45° mit Endscheiben. Fliigel B = 0°, 45° ohne Endscheiben. 


Aus den Messungen kinnen folgende Ergebnisse entnommen werden: 


1. Auftrieb. Die Auftriebskrafte des schiebenden Fliigels mit Endscheiben sind hei 
gleichem geometrischen Anstellwinkel in Ubereinstimmung mit den theoretischen Uberlegungen 
des Abschnittes 5 proportional dem Cosinus des Schiebewinkels: 


€q = (€q)p= 9-008 p'. 
Dieser Zusammenhang ist fiir den Gesamtauftrieb aus den Auftragungen von c, iiber a zu ent- 
nehmen (Abb. 3), in welche die aus den Werten des geraden Fliigels 6 = 0 umgerechneten Werte 
gestrichelt eingezeichnet sind. 

Dasselbe Gesetz findet auch seine Bestitigung fiir den drtlichen Auftrieb des Fliigels langs 
Spannweite. Diese Werte wurden gewonnen aus den Druckverteilungsmessungen an den einzelnen 
Profilschnitten. Wie die Auftragungen der ¢,5,-Werte langs Spannweite zeigen (Abb.7), ist fiir 
die Fliigel mit Endscheiben bei kleinen Anstellwinkeln die Auftriebsverteilung lings Spannweite 
mit recht guter Genauigkeit konstant. Bei gréBeren Anstellwinkeln treten jedoch infolge von 
Grenzschichteinfliissen Abreiferscheinungen an dem rickliegenden Ende des Fliigels auf. Die 
gestrichelten Kurven in Abb.7 sind aus den entsprechenden gestrichelten Mittelwerten des ge- 
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raden Fliigels durch Multiplikation mit dem Faktor cos f erhalten worden. Man sieht, daB die 
Proportionalitat mit cos § im gesunden Anstellwinkelbereich sehr gut erfiillt ist. 

Weiter wurde untersucht, ob die Druckverteilung am Profil, also der Druck in jedem einzelnen 
Punkt, ebenfalls dem obigen cos f-Gesetz unterliegt. Da die Druckverteilungen bei anliegender 


[ero 
Jo 


40 


Abb. 5. Druckverteilungen tiber Fliigeltiefe fiir 
den Mittelschnitt (Nr. V) der Fliigel 8B = 0°, 30°, 
45° mit Endscheiben. a’ = 3°, 6°, 9°. 


Strémung an den einzelnen Schnitten gut iibereinstimmen, haben wir diese Kontrolle nur fiir 
den Mittelschnitt durchgefihrt, fiir den in Abb.5 die Druckverteilungen aufgetragen sind. Bei 
Giiltigkeit der Proportionalitat des értlichen Druckes mit cos, miissen die Verteilungen fiir 
alle Fliigel zusammenfallen, wenn wir (p— po)/cos f itber der Tiefe auftragen. Wie Abb. 6 zeigt, 
liegen diese umgerechneten Werte tatsiachlich fiir alle drei Fliigel auf einer einzigen Kurve. 
Hiermit ist also bewiesen, da® die gesamte Druckverteilung am Filiigel sich mit dem Schiebe- 


winkel nach dem cos /-Gesetz andert. 
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Der Gesamtauftrieb des schiebenden Fliigels ohne Endscheiben (f =45°) folgt ebenfalls 
niherungsweise dem cos f-Gesetz (vgl. Abb.4). Beziiglich der Auftriebsverteilung langs Spannweite 
(Abb.8) miiBte eine genaue Auftriebsverteilungsrechnung durchgefiihrt werden, um einen Vergleich 
zwischen dem nicht schiebenden und dem schiebenden Fliigel durchfithren zu kénnen (vg). Abschn.9). 


pejto 7 VARIG pe 
7 cosh 10 Jo cos 
7 


Abb. 6. Druckverteilungen iiber Fliigeltiefe fiir 
den Mittelschnitt (Nr. V), umgerechnet auf den 
geraden Fligel mit dem Faktor 1/cos. Fligel 
B = 0°, 30°, 45° mit Endschciben. a’ = 3°, 6°, 9°. 


2. Neutralpunktlage. Der Neutralpunkt bleibt bei den Fliigeln mit Endscheiben 
fiir samtliche Schiebewinkel annahernd im 1/4-Punkt des Mittelschnittes erhalten, wie aus Abb. 3 
hervorgeht. Es ergibt sich fir alle Fliigel dey/de, = 0,01, bezogen auf den 1/4-Punkt des Mittel- 
schnittes; der Neutralpunkt liegt also auf 24°% der Fligeltiefe hinter der Fliigelnase. Die Ergeb- 


nisse der Fliigel ohne Endscheiben sind fir den schiebenden Fligel leider verloren ge- 
gangen. 
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3. Profilwiderstand lings Spannweite. Einen besonderen EinflufB hat der Schiebe- 
winkel auf die Grenzschicht und damit auf die Verteilung des Profilwiderstandes langs Spann- 
weite des Fliigels. Die Grenzschichtbeeinflussung rithrt daher, daB beim schiebenden Fliigel in- 
folge der in Strémungsrichtung vorhandenen Verschiebung der Profile ein Druckgradient in der 
Querrichtung entsteht, dem vor allem das abgebremste Grenzschichtmaterial infolge seiner ge- 
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: indigkeit und der daraus resultierenden geringen Trigheitskrafte zu folgen ver- 
Sect cchiober ded Fliigel tritt auf der Saugseite des Fliigel nach dem zuriickliegenden ae 
ende eine Druckabnahme auf, so daB das Abwandern des Grenzschichtmaterials in prs Bees 1- 
tung erfolgt und daher eine Anhadufung im riickhegenden duBeren Fliigelteil erzeugt ne ; oes 
bedeutet erhéhten Profilwiderstand und erhéhte Abreifgefahr in diesem Gebiee Deane os ad 
stellwinkeln. Der EinfluB ist umso starker, je gréBer der Schiebewinkel des Fligels ee ie r- 
scheinung geht deutlich aus unseren Messungen hervor (vgl. Abb.3 und 4). Der Fliigel mit gréStem 


424, Jacobs: Experimentelle Untersuchungen am schiebenden Fliigel. Ingenieur-Archiv 


Schiebewinkel hat im oberen c,-Bereich die gréBten Widerstande, und seine Polare knickt dort 
zuerst ab. Der eingetragene Widerstand stellt die Summe von Fliigelwiderstand und Widerstand 
der Endscheiben und der Aufhangung dar. 

AufschluB itber die Profilwiderstande lings Spannweite geben die bei mehreren Anstellwinkeln 
in verschiedenen Schnitten hinter den Fliigeln durchgefithrten Impulsmessungen. Die hieraus er- 
mittelten Profilwiderstandsbeiwerte sind fir die drei Fliigel in Abb. 9 aufgetragen. Die Messungen 
wurden ausgefithrt im Abstand einer halben Fliigeltiefe von der Hinterkante. Die Auswertung 
erfolgte nach der Formel yon B. M. Jones. Die Auftragungen zeigen beim geraden Fliigel wie 
zu erwarten, konstante ¢,,-Werte langs der Spannweite. Bei den schiebenden Fliigeln da- 
gegen tritt mit wachsendem Anstellwinkel ein immer stérkeres Anwachsen des értlichen ¢,- 
Wertes nach dem zuriickliegenden Ende hin auf. Diese Zunahme erfolgt fast linear bis kurz vor 
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Abb. 8. cq ,4-Verteilungen lings Spannweite fiir die Fliigel 8 = 0° und f = 45° ohne Endscheiben. 
ot I (i Uy ee SG: 


die Stelle, wo zuerst AbreiBen der Strémung eintritt. Besonders bemerkenswert ist, daf vor 
der Stelle des eigentlichen AbreiBens der Strémung eine Abnahme des Ortlichen Profilwider- 
standes zu beobachten ist. Dies rithrt offenbar daher, daB der Druckgradient in Querrichtung, 
der im vorderen Teil des Profiles in der Gegend des Druckminimums am starksten ist, an 
der Grenze zwischen abgerissenem und anliegendem Gebiet infolge des AbreiBens iiber ein 
langeres Stiick der Fliigeltiefe wirksam ist und daher von dem noch anliegenden Grenzschicht- 
material des Nachbarschnittes besonders viel abgesaugt wird. Infolgedessen sinkt der Profil- 
widerstand an der betreffenden Stelle. 


Die Zunahme der Grenzschichtdicke am Fliigel 6 = 45° ist in Abb.10 fiir den Anstellwinkel 
a’ = 14° dargestellt. Es ist hier das Ergebnis einer Gesamtdruckmessung dicht hinter der Hinter- 
kante des Fliigels aufgetragen. Die Messungen wurden an drei verschiedenen Schnitten am zu- 
riickliegenden Teil des Fliigels durchgefithrt. Die Koordinate z ist positiv nach unten gemessen. 
Die Zunahme der Grenzschichtdicke erfolgt ausschlieBlich auf der Saugseite, wahrend die Druck- 
seite nur geringe Anderungen aufweist. Hieraus geht hervor, daS der Einflu® der ,,Grenzschicht- 
abwanderung“ sich nur auf der Saugseite abspielt. 
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5. Theoretische Uberlegungen. Bei der Diskussion der Ergebnisse der Auftriebsmessungen fiir 
die Fliigel mit Endscheiben im vorigen Abschnitt wurde gefunden, daB sich die Anitiel tase 
bei gleichem in Anstrémungsrichtung gemessenen Anstellwinkel um den Faktor cos f Seo tlie 
den Werten des nicht schiebenden Fliigels (6 =0) verkleinern. Dieses cos f-Gesetz laBt aio fol- 
genderma8en leicht theoretisch begriinden: 

Wir nehmen eine zwischen Endscheiben befindliche ebene Platte an, die zur Anstrémungs- 
richtung um den Schiebewinkel f schraggestellt ist, und die mit der Anstrémungsrichtung Hie 
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Abb. 9. Ortlicher Profilwiderstand lings Spannweite bei verschiedenen Anstellwinkeln 
fiir die Fligel 8 = 0°, 30°, 45° mit Endscheiben. 


Anstellwinkel « bildet. Da durch die Endscheiben cine seitliche Umstrémung der Fliigelenden 
verhindert wird, kann geschlossen werden, daB der Druck auf Linien gleicher prozentualer Riick- 
lage langs Spannweite konstant ist. Es liegt also ein sogenanntes ,,quasiebenes“ Problem vor. 
Denken wir uns jetzt den Fliigel durch einen Potentialwirbel ersetzt, so kann dieser gemaB dem 
eben Gesagten nur in Richtung des schiebenden Fliigels liegen. Wird nun dieser Wirbel mit 
der Anstrémungsgeschwindigkeit v angeblasen, so ist nach Kutta-Joukowsky jedoch nur die 
senkrecht zur Wirbellinie herrschende Geschwindigkeitskomponente, in unserem Fall v* = vcos f, 
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fir die Auftriebserzeugung wirksam. Desgleichen ist auch der Anstellwinkel «* senkrecht zur 
Wirbellinie zu messen. Fiir kleine Anstellwinkel gilt a* = a/cos f. 

Die gesamte Auftriebskraft der unter dem Anstellwinkel «* stehenden ebenen Platte, die mit 
der Geschwindigkeit v* angestrémt wird, ergibt sich zu 


A= 2 ma on ee 
Beachten wir, da®B F* = F die Flache des schiebenden Fliigels ist, dann erhalten wir unter Ein- 


setzen der Beziehungen fiir v* und «* den Auftriebsbeiwert zu 
A 


iz ¢, = =2n0cosp. 
ah | ae 
i | B=45° ( 
-80 | Wird akonstant gehalten, dann 
N ie ergibt sich daraus das gesuchte 
71 ieee) if cos B-Gesetz 
‘Av teen ; C5 = C56 cos fb 
= 4Z 
4 bzw. fiir den Auftriebsanstieg 
2 50. | de, de, 
5 -(fae 
8 yp : , 
| Werden jetzt die Endscheiben 
oe | entfernt, so wird aus dem 
\ | IN P, eben behandelten _ ,,quasi- 
200m AS ebenen*’ Problem ein raéum- 
a AQ liches Problem. Fiir die Auf- 
~ om = triebsverteilung langs Spann- 
10 : weite ist jetzt eine genauere 
Se aw 7- fe Rechnung ; ERE a 
Hierzu eignet sich z. B. das 
Verfahren von J. Weissinger?. 
& 
8 te 6. Zusammenfassung. Es 
§ af wird itiber Kraft-, Druckver- 
| teilungs- und Impulsverlust- 
messungen an_ schiebenden 


30 : ; ; ; 
nea : Fligeln, die sich zwischen 
Abb. 10. Grenzschichtdicke an der Hinterkante des Fliigels 8 = 45° mit Endscheiben : : : 
fiir verschiedene Schnitte lings Spannweite. Endscheiben befinden, berich- 


tet. Beim schiebenden Trag- 
fliigel unendlicher Spannweite andern sich Gesamtauftrieb, értlicher Auftrieb und Druck- 
verteilung nach dem cos B-Gesetz (f/f = Schiebewinkel). Der Neutralpunkt liegt fiir samtliche 
Schiebewinkel annihernd im 1/4-Punkt des Fliigelmittelschnittes. Impulsmessungen zeigen, daB 
der Profilwiderstand langs Spannweite bei héheren Anstellwinkeln infolge der ,,Grenzschicht- 
abwanderung* eine betrachtliche Zunahme nach dem zuriickliegenden Fliigelende hin erfahrt. 
Dieser Effekt wachst mit dem Schiebewinkel des Fliigels staik an. 


(Eingegangen am 2. Mai 1952.) 
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